SINGULAR  LINE  THEORY  AND  CONTROL  SYSTEMS 


Constantinos  Spyros  A,  Cariniotakis 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

SINGULAR 

LINE  THEORY  AND  CONTROL 
by 

SYSTEMS 

Constantinos  Spyros  A.  Cariniotakis 

March,  1978 

Thesis  Adv 

isor:           George  J 

.  Thaler 

Approved  for  public  release;  distribution 
unlimited 


T 183199 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Data  Bntarad) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  OOVT  ACCESSION  NO 


1.     RECIPIENT'S  CATALOG  NUMBER 


4.     TITLE  (and  Subllll,) 

Singular  Line  Theory  and  Control  Systems 


5.     TYPE  OF  REPORT  •  PERIOD  COVERED 

Electrical  Engineer 
March.  1978 


•  .  PERFORMING  ORG.  REPORT  NUMBER 


7.   AUTMORf.) 

Constantinos  Spyros  A.  Cariniotakis 


•   CONTRACT  OR  GRANT  NUMBERS 


»   PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


10.    PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


II.     CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


12.     REPORT  DATE 

March   1978 


IS.  NUMBER  OF  PAGES 


296 


U.    MONITORING  AGENCY  NAME  a  AODRESSfi/  dlllarant  from  Controlling  Olllca) 

Naval  Postgraduate  School 
Monterey,  California  93940 


IS.    SECURITY  CLASS,  (ol  (hi.  riporl) 


UNCLASSIFIED 


ISa.     OCCLASSIFI  CATION/ DOWN  GRADING 
SCHEDULE 


16.     DISTRIBUTION  STATEMENT  (ol  (Ma  Raport) 


Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (ol  tha  abattact  anlarad  In  Block  20,  II  dlllmronl  ham  Raport) 


IS.    SUPPLEMENTARY  NOTES 


IS.    KEY  WORDS  (Contlnum  on  tavaraa  alda  II  nacaaaary  ltd  Idantltr  «r    block  numbar) 


Parameter  plane 
Parameter  space 
Singular  system 


Singular  line 
Singular  point 
Singular  root 


20.    ABSTRACT  (Contlnua  on  rovaraa  aldm  II  nacaaaary  and  idantltr  by  aloe*  maaaar) 

The  theory  of  singular  lines  on  the  Parameter  Plane  and  in 
Parameter  Space  are  derived,  and  applied  to  control  system  design. 

A  method  for  design  of  compensation  of  linear  feedback  control 
systems  using  singular  line  theory  is  presented, and  application  of 
this  design  method  for  self-adaptive  control  systems  is  considered. 

General  design  steps  and  required  procedures  for  using  this 
design  method  are  summarized  and  examples  are  presented. 


DD  ,!i;:M7l  1473 

(Page  1) 


EDITION  OP   I  NOV  St  IS  OBSOLETE 

S/N   0102-014-6601   | 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dmlm  Bntarad) 


Approved  for  public  release;  distribution  unlimited 


SINGULAR  LINE  THEORY  AND  CONTROL  SYSTEMS 


by 


Constantinos  Spyros  A.  Cariniotakis 
Lieutenant  Commander,  Greek  Navy 
B.S.E.E.,  United  States  Naval  Postgraduate  School,  1977 
M.S.E.E.,  United  States  Naval  Postgraduate  School,  1977 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

ELECTRICAL  ENGINEER 


from  the 
NAVAL  POSTGRADUATE  SCHOOL 
March  19  7j^^. 


.\ 


ABSTRACT 


The  theory  of  singular  lines  on  the  Parameter  Plane 
and  in  Parameter  Space  are  derived,  and  applied  to  control 
system  design. 

A  method  for  design  of  compensation  of  linear  feed- 
back control  systems  using  singular  line  theory  is 
presented,  and  application  of  this  design  method  for 
self-adaptive  control  systems  is  considered. 

General  design  steps  and  required  procedures  for  using 
this  design  method  are  summarized  and  examples  are  pres- 
ented. 
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INTRODUCTION 


Prior  to  1959  analysis  and  design  of  control  systems 
was  primarily  performed  by  classical  methods.   These 
methods  are  adequate  when  there  is  a  single  variable 
parameter  and  a  single  loop  system. 

In  modern  control  engineering,  where  multiloop  systems 
with  characteristic  equations  with  two  or  more  variable 
parameters  usually  result,  the  above  mentioned  design 
methods  become  in  general  very  complex.   The  designer  has 
to  complete  several  designs  in  a  trial  and  error  process 
and  to  choose  the  most  suitable  solution  and  to  make 
usually  some  arbitrary  choices  based  on  his  experience. 

Methods  for  studying  the  effect  of  variations  in 
parameters  of  a  multivariable,  multiloop  control  system 
have  been  developed  by  Mitrovic  [1]  ,  Siljak  [2]  ,  Thaler 
[3]  ,  and  other  investigators.   All  these  methods  are  based 
on  the  same  concept  and  their  dependence  on  a  computer  is 
characteristic. 

Mitrovic "s  [l]  method  which  was  introduced  in  1958, 
allowed  study  of  the  effect  of  variations  of  two  variable 
parameters  on  the  location  of  the  characteristic  roots  of 
a  system.   This  method  consists  of  specifying  as  variables 
the  two  lowest  order  coefficients  of  the  characteristic 
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equation.   Using  then  the  characteristic  equation  as  a 
mapping  function,  constant  zeta,  omega  and  sigma  curves 
are  transformed  from  the  s-plane  into  the  variable  co- 
efficient plane.   This  method  permits  the  designer  to 
adjust  these  variable  coefficients  so  that  the  roots  of  the 
characteristic  equation  may  be  set  at  desired  locations. 

In  1964,  Siljak  [2]  extended  Mitrovic's  [1]  ,  co- 
efficient plane  method  into  the  Parameter  Plane  method  in 
which  the  two  variable  parameters  may  appear  linearly  in 
any  or  all  coefficients  of  the  characteristic  equation. 
The  parameter  plane  method  provides  the  designer  with  a 
procedure  for  factoring  a  characteristic  equation  and 
displaying  the  results  in  a  parameter  plane  diagram,  and 
allows  him  to  obtain  information  about  system  stability 
and  how  it  is  affected  with  parameter  variations. 

The  parameter  plane  method  was  further  extended  in 
1965  [4]  to  a  general  case  in  which  the  two  variable 
parameters  may  appear  nonlinearly  in  the  coefficients  of 
the  characteristic  equation.   Another  extension  of  the 
Parameter  Plane  method  to  include  more  than  two  variable 
parameters  has  been  defined  as  the  Parameter  Space  method 
[5]  . 

When  in  general  only  two  parameters  (or  coefficients) 
are  involved,  for  any  complex  value  of  s  or  a  pair  of  (c,w  ) 
values  two  simultaneous  equations  in  two  unknowns  are  ob- 
tained by  setting  the  real  and  imaginary  part  of  the 
characteristic  equation  equal  to  zero. 
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In  1967,  Bowie  [6]  while  working  with  the  character- 
istic equation  of  an  inertially  stabilized  vehicle  on  the 
parameter  plane,  found  that  a  complete  set  of  roots  for 
the  characteristic  equation  could  not  be  determined  in 
the  real  root  area  by  using  the  existing  parameter  plane 
theory.   He  further  investigated  this  problem  and  found 
that  specific  pairs  of  C,un  values,  make  the  above  men- 
tioned system  of  equations  singular,  which  thus  corresponds 
to  a  single  straight  line  in  the  parameter  plane,  called 
a  "singular  line". 

The  singular  line  added  a  new  dimension  to  the  para- 
meter plane  because  it  provides  the  designer  a  way  to  hold 
a  root  fixed  while  varying  others,  by  varying  the  two 
parameter  values.   Another  important  characteristic  of 
singular  lines  is  also  the  fact  that  they  can  be  drawn  in 
a  parameter  plane  diagram  by  hand,  without  needing  a 
computer,  as  in  the  case  of  drawing  constant?  orcon  curves. 

Although  an  analytic  study  of  the  two  previous  men- 
tioned simultaneous  equations,  which  are  obtained  by 
setting  the  real  and  imaginary  part  of  the  characteristic 
equation  equal  to  zero,  makes  obvious  the  existence  of 
singular  lines  under  certain  conditions,  they  were  not 
found  earlier  due  to  the  standard  computer  programming 
practices  and  limitations. 

Singular  line  theory  is  based  on  the  special  case  of 
linearly  dependent  parameter  plane  solution  equations. 

From  the  standpoint  of  the  design  engineer  it  was 
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desirable  that  this  theory  be  investigated  in  order  to 
establish  a  standard  general  set  of  rules  which  would 
provide  a  method  to  obtain  systems  with  singular  lines. 

In  this  thesis  this  task  is  accomplished  and  possible 
applications  of  the  singular  line  theory  in  control  system 
design  are  investigated. 
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II.    SINGULAR  LINE  THEORY 


A.    SINGULAR  CONDITIONS 

Consider  the  characteristic  equation: 

n 
F(s)  =  Ls     a<  s<  =o  (2-1) 

K=0 

where  a  ^    {k=0,    1,    2,...,n}    ,    is   a   function  of   two  variable 
parameters  a  and  3  ,    i.  e.  ,   a  jc=g}c    (a,  3)  ,    and   s    is   a  point   on 
the   s   plane  which  can  be   expressed   in   rectangular   coordi- 
nates  as: 


s=  -  ^iJ^Vl-?2  (2-2) 

If  for  a  given  pair  of  ?s,  o)ns  values  {where  |.5S|   <1  and 
uns  >®  ^  »  equation  (2-1)  has  an  infinite  number  of  a  and  3  real  valued 
pairs  which  satisfy  it,   then  for  this  equation  singular  conditions 
exist.     In  this  case  the  point  of  the  s  plane  which  is  defined  by  the 
pair  £    ,  oons  is  a  singular  point,  and  the  characteristic  equation  is 
singular. 

The  locus  of  all  points   (  a,   3)  on  the  Parameter  Plane  which 
satisfy  equation  (2-1)   in  the  case  of  singular  conditions  is  called 
the  singular  line.     Each  such  line  corresponds  to  a  specific  singular 
point  or  in  other  words  to  a  specific  pair  of  c    ,  cons  values. 
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B.   SINGULAR  LINE  ANALYSIS  IN  THE  LINEAR  CASE 

A  simple  case,  which  was  analyzed  and  which  appears  in 
many  practical  problems  is  the  linear  case.   In  this  case 
the  two  variable  parameters  a  and  g  appear  in  the  coefficients 
of  the  characteristic  polynomial  in  the  general  form: 

V=  3Ka+cK  B+dK  (2-3) 

wherebK,cK    and  dK  are   real   constants. 

1.    Analysis  by  Using  Parameter  Plane  Equations 

Consider  the  characteristic  equation  (2-1),  where 
aK    is  defined  by  equation  (2-3).   By  introducing  equation 
(2-2),  into  equation  (2-1),  and  equating  real  and  imaginary 
parts  to  zero,  it  can  be  shown  [6]  that: 

V^  3   -   -  D,  (2.4) 

B2  a  +  C2  6    -    -   D2 

where 
n 


Bi-£  (-DKbK   %<   UK.1  CO         B,-r     (-l)\u)nKdK   UK(?) 

K  =  0  K^O 

Ci-£  (-1)%   a)nK  Uk-i  CO         C2=^     C-DKcKwn<uK    CO         (2-5) 

K  =  0  K  =  0 

d,-.£  l-l)\  %K  u,,.,  CO      D,-f;    (-i)%;:B[  t« 


K  =  0  K=0 
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The  functions  U,  (5)  and  uY_i(?)  are  Chebyshev 
functions  of  the  second  kind  for  which  U  (C)-0,  U,(?)=l, 
and  in  general : 

UkCC>  2cUk_1(C)  -  Uk_2C?)  (2-6) 

The  system  (2-4)  can  be  written  in  the  Matrix  form 


as 


or 


Bl   Cl 

B2   C2 


A  x  =  B 


D. 


-D- 


(2-7) 


When  the  solution  of  the  system  (2-7)  is  considered 
there  are  the  following  cases: 

a.   No  Solution 

When  Rank  (A)  =  1  and  Rank  (A,  B)  =  2,  which 
implies  that: 


i  =  c_i  +  D_i 


B2    C2 


D, 


(2-8) 


then  the  system  has  no  solution, 
b.   Unique  Solution 

When  Rank  (A)  =  2  and  Rank  (A,  B)  =  2 ,  which 
implies  that: 


Bl  .  Cl 


D, 


B2    C2      U2 
then  the  system  has  a  unique  solution. 


C2-9) 
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c.   An  Infinite  Number  of  Solutions 

When  Rank  (A)  =  1  and  Rank  CA,  B)  ■  1,  which 
implies  that 

Bl    Cl     Dl 

-i-  =  —  =   —  C2-10) 

B2    C2      D2 

then  the  system  has  an  infinite  number  of  solutions.   There- 
fore in  this  case,  since  there  are  an  infinite  number  of 
pairs  of  a, 6  values  which  satisfy  the  system  (2-4)  and  in 
consequence  the  characteristic  equation  (2-1),  it  implies 
that  for  this  equation  singular  conditions  exist. 

2 .    Analysis  by  Using  Complex  Variable  Theory 

The  characteristic  equation  (2-1)  is  considered 
again.   This  equation  can  be  rearranged  in  the  form: 

f^s)  a+  f2(s)  3  +  f3(s)  =  0  (2-11) 

where  f-,(s),  f2(s)  and  f,(s)  are  polynomials  of  s.   Intro- 
ducing equation  (2-2)  into  equation  (2-11),  these  poly- 
nomials will  be  in  general  a  complex  number  with  some 
magnitude  and  an  argument.   Thus  they  can  be  written  as: 

fx(s)  =  r1   ejei 

f2(s)  =  r2  ej92 

f3(s)  =  r3  e^3 

Introducing  the  above  equations  into  equation 
C2-11)  yields: 

r1   e^ia  +      r2    eje  2  3+      r3   eje  3   =    0  (2-12) 
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Using   the   Euler's    formula   eJ      =   cos9    +   jsine    , 
equation    (2-12)    becomes;  i 

r,{ COS0-.+J sine,  }0i+r2{cose2+j sin02}3+r3{ cos 9,+j sine,  }  -  0 

or 

{r,  cosS  ^a+r-cose^g+r^cose^^j^r,  sine  1a+r2sine23+r,sine^}  =  0 

Equating   real    and   imaginary  parts    to    zero   yields 
the    system: 

r, sin6ia+r2sin92  g=-r,sine-7 

(2-13) 
r,cos6  a+r2cos62  ft—r^cose^ 

By   applying    the    same   reasoning   as    in  paragraph   one, 
the    system    (2-13)    will    have: 

a.    No   Solution 

If  and  only   if: 

r,sin6,       r2sine2       r,sin6.. 


or 


r,cose-i       r2cose2       r^cose? 

tan8,    =   tane2   f   tane, 

b.  Unique  Solution 

If  and  only  if: 
tane1  f   tan62  f   tane.. 

c.  An  Infinite  Number  of  Solutions 
If  and  only  if: 

tan91  =  tane2  =  tane^  (2-14) 
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Equation  (2-14)  implies  that  for  the  character- 
istic equation  (.2-1)  which  has  been  rearranged  in  the 
form  (2-11)  singular  conditions  exist,  if  and  only  if: 

91  *  e3  +  kin 

e2  -  e3  +  k2n 

or 


Acs)  =A(s)  *  k2n 

where  k,  and  k2  are  integers. 


(2-15) 


3.    Singular  Conditions  for  the  Linear  Case 

Summarizing  previously  derived  results,  the  con- 
clusion is  that,  in  the  linear  case  where  the  character- 
istic equation  is  written  in  the  general  form: 


n 


F(s)  =  T*  a,.  sk   0 


k=0 
where  a.,     =  g,  a+  cl-B+  ^v»  singular  conditions  exist  if  and 
only  if: 

B2    C2    D2 

When  the  characteristic  equation  is  considered 
in  the  general  form: 

f,(s)a  +  f9(s)  3  +  f_(s)  =  0 
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then  singular  conditions  exist,  if  and  only  if; 

/f^s)  t^ Cs)  +  krn 
/f2(s)  =/3Cs)  +  k2n 

4.    Equation  of  a  Singular  Line 

When  for  a  given  pair  oft,  u   values  singular 
conditions  exist  for  the  characteristic  equation: 

n 

F(s)  =E  ak sk  =  ° 

k=0 


or 

f1(s)a  +  f2(s)B  +  f3(s)  =  0 

then,  according  to  the  previous  analysis  the  two  para- 

metersot  and3  will  be  related  by  one  of  the  following 
equations,  respectively: 

D-.    B,       D^  B0 

9>   =  .  _L  .  _L  a=.  _i  .  _i  a 

Cl    Cl      C2  C2 


(2-16) 


or 


rTsin9_   r^sine,     rTcos0T   r-,cos9i 
B  =  -  _J £.  .  J: ia  =  -  _J £.  -  _i £  a   (2-17) 

r2sin62   r2sin9  2    r2cos62   r2cos02 


Each  of  these  equations,  is  the  equation  of  the 
singular  line  corresponding  to  the  singular  pair  (5  ,  (Dnc.) 
which  was  considered.   This  singular  line  is  a  straight 
line  on  the  parameter  plane  which  is  tangent  to  the  constant 
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5   curve  at  the  singular  point  [6].   It  was  found  in  this 
thesis,  that  this  is  also  tangent  to  the  constant  us   curve 
at  the  same  point. 

5 .    Methods  to  Test  if  a  Point  s  is  Singular 

When  the  characteristic  equation  is  given  and  it 
is  desired  to  determine  if  a  given  point  C  C>  <*>_]  ,  (|c|  <1 
and  ui  >  0)  ,  is  singular  the  following  methods  can  be  used. 

a.  Algebraic  Methods 

By  using  the  C  ,  w  values  of  the  point  under 
consideration,  the  coefficients  of  equations  (2-5)  can  be 
evaluated.   Then  the  given  point  will  be  singular  if  and 
only  if  equations  (2-10)  are  satisfied. 

A  second  algebraic  method  consists  of  the 
evaluation  of  the  arguments  appearing  in  equation  (2-15)  , 
after  writing  the  characteristic  equation  in  the  form  of 
equation  (2-11).   Then  the  given  point  is  singular  if  and 
only  if  equations  (2-15)  are  satisfied. 

b.  Graphical  Method 

By  this  method  the  characteristic  equation 
must  be  written  in  the  form  of  equation  (2-11)  and  the  roots 
of  the  polynomials  f,(s),  f2Cs)  and  f?(s)  must  be  evaluated. 
After  these  roots  have  been  located  on  the  s-plane,  with 
the  use  of  a  spirule  the  arguments  of  the  previous  men- 
tioned polynomials  for  the  point  under  consideration  are 
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determined.   Then  this  point  will  be  singular  if  and  only 
if,  equations  C2 -15)  are  satisfied. 

c.   Use  of  a  Computer 

There  are  in  general  many  ways  by  which  a 
computer  can  be  programmed  to  answer  this  problem.   One 
such  program  is  the  << Singular  Point  Program  >>  which  was 
developed  with  this  thesis  and  which  appears  in  Appendix 
H.   In  order  for  this  program  to  be  used  in  the  case  under 
consideration,  the  corresponding  value  of  z,   of  the  point 
which  is  considered  is  used  as  an  input  data.   Then  at  the 
computer  print  output  the  corresponding  value  of  w   is 
checked  if  it  is  such  that  it  makes  the  system  (2-4)  have 
an  infinite  number  of  solutions.   When  this  is  the  case, 
the  point  (  C,  w )  which  was  considered  is  a  singular  point, 
otherwise  it  is  not  singular. 

6.   Methods  to  Test  if  a  Characteristic  Equation  is 
Singular 

In  the  previous  paragraph  a  characteristic  equation 
was  considered  to  be  given  and  it  was  desired  to  determine 
if  a  point  of  the  s-plane  was  a  singular  point  for  this 
equation.   In  this  paragraph  the  problem  is  reversed,  i.e., 
a  characteristic  equation  is  considered  to  be  given  and  it 
is  desired  to  determine  if  there  is  at  least  one  point  of 
the  s-plane  which  is  singular.   For  this  problem  the  fol- 
lowing methods  can  be  applied: 
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a.   Analytic  Method 

By  this  method  the  characteristic  equation  is 
written  in  the  form  of  equation  (2-11).   Then,  the  argu- 
ments of  f, (s) ,  f2Cs)  and  f^Cs)  are  evaluated  as  functions 

of  S  and  u    i.e., 

n      ' 


/f^s)  =  §1c  c,^n) 


'f2Cs)  =  g2(  c,o)n)  (2-18) 


/f3Cs}  =  *3C  ?'wn) 

Introducing  equations  (2-18)  into  the  system 
(2-15)  yields: 

g1(C,wn)  =  g3(c,wn)  +  kxn 

(2-19) 
§2U,%)  -  g3(?,wn)  +  k2n 

When  the  system  (2-19)  is  solved  for  the  un- 
knowns? and  a)  ,  any  acceptable  solution,  i.e.,  any  solution 

where  the  values  of?  and  oj   are  real  and  also   C  <   1 

n 

and  u  >  0,  will  define  a  pair  of  C,  w  values,  which  cor- 

n    '  r  '      n 

respond  to  a  singular  point  of  the  characteristic  equation. 
If  at  least  one  such  point  is  found  then  the  character- 
istic equation  is  singular. 

b.   Graphical  Method 

By  this  method  the  characteristic  equation 
which  is  considered  must  be  rearranged  again  in  the  form 
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of  equation  (2-11).   Then  the  roots  of  the  polynomials 
f,(s),  f?(s)  and  f,(s)  are  evaluated  and  they  are  located 
on  the  corresponding  points  of  the  s-plane.   By  the  use  of 
a  spirule  a  search  on  the  s-plane  must  be  performed  for 
points  which  satisfy  equations  (2-15)  .   This  search  can  be 
performed  either  by  moving  on  constant  Z,   lines  or  constant 
co   cycles.    This  method  will  in  general  be  useful  when 
the  polynomials  of  s  which  are  involved  in  equation  (2-11) 
are  of  low  order  and  it  is  desired  to  find  out  if  the  given 
characteristic  equation  is  singular  for  a  restricted  area 
of  the  s-plane  which  is  of  interest. 

c.   Use  of  a  Computer 

The  Singular  Point  Program  which  was  mentioned 
before  and  which  was  developed  for  this  purpose  corresponds 
to  one  among  the  many  ways  by  which  a  computer  can  be  pro- 
grammed to  answer  this  problem. 
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III.    SINGULAR  LINE  THEORY  AND  CONTROL  SYSTEMS 

A.    DEFINITIONS 

For  a  control  system  with  two  parameters  a  and  g 
appearing  in  its  characteristic  equation  the  following 
definitions  come  as  a  consequence  of  the  singular  line 
theory  which  was  previously  developed. 

1.  Singular  System 

A  control  system  is  singular,  if  and  only  if, 
its  characteristic  equation  is  singular. 

2 .  Singular  Characteristic  Equation 

The  characteristic  equation  of  a  control  system 
is  singular,  if  and  only  if,  the  system  of  the  two 
simultaneous  equations  in  two  unknowns,  a  and  3,  which  is 


obtained  by  setting  s  =  -  c;w  ±  jto _\/l-  S   into  the  character- 
istic equation  and  requiring  the  summation  of  reals  and  the 
summation  of  imaginaries  to  go  to  zero  is  singular,  for  at 
least  one  pair  of  r,  ,  w   values,  where   |c|  <1  and  ajn  >0  . 
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3.  Singular  Root 

A  root  of  the  characteristic  equation  of  a 
control  system  is  singular  if  it  makes  the  characteristic 
equation  singular. 

4.  Singular  Point 

Each  singular  root  corresponds  to  a  specific 
point  of  the  s  plane.   Such  a  point  is  called  a  singular 
point. 

B.    SINGULAR  LINE  THEORY  IN  CONTROL  SYSTEM  DESIGN 

In  general  the  analysis  and  design  of  any  control 
system  is  based  on  an  accepted  model  for  the  system  under 
consideration.   In  the  case  of  a  linear  control  system  of 
any  order  the  no  zero,  two  complex  conjugate  poles  system 
model  is  used  as  a  standard  one.   The  accuracy  of  this 
approximation  is  tightly  related  to  the  dominant  mode 
concept. 

At  this  point  the  use  of  the  singular  line  theory 
seems  to  be  of  utmost  importance  since,  if  a  system  is 
singular  the  complex  pair  which  provides  the  singularity 
remains  fixed  and  the  locations  of  the  other  poles  may  be 
adjusted  to  some  extent  because  an  infinite  number  of  a 
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and  3  pairs  is  available  for  their  adjustment.   In  addition 
some  other  performance  criteria  may  be  satisfied  by  an 
appropriate  parameter  variation. 

All  these  cases  are  analyzed  and  further  investigated 
for  other  probable  applications  of  the  singular  line  theory 
in  control  system  design. 

C.    SINGULAR  LINE  THEORY  AND  CONTROL  SYSTEMS 

1 .   Basic  Requirements 

From  all  that   has  been  stated  in  section  II,  it 
seems  that  in  order  for  the  singular  line  theory  to  be 
applied  on  control  systems  design  and  compensation  the 
following  basic  requirements  must  be  satisfied,  i.e.,  the 
system  under  consideration  must  be  singular,  a  singular 
root  of  the  system  characteristic  equation  must  meet 
problem  specifications  when  it  is  considered  the  dominant 
root,  this  root  must  be  made  dominant  by  adjusting  the 
parameter  ot  and  3  values,  and  finally  the  two  parameters 
a  and  3  must  be  adjustable. 

.  Therefore,  the  initial  problem  in  general  will  be  to 
select  the  structure  of  the  system  or  the  required 
compensator  for  a  given  system,  in  order  for  the  above 
basic  requirements  to  be  met. 

In  the  following  paragraph  a  preliminary  analysis  of 
the  form  of  the  characteristic  equation  of  a  control  system 
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was  performed,  as  the  first  step  in  developing  the  re- 
quired structure  of  the  control  system  to  obtain  the 
desired  singular  conditions. 


2.   The  Characteristic  Equation  of  a  Control  System 
in  the  Linear  Case 


In  general  the  characteristic  equation  of  a  con- 
trol system  in  the  linear  case  will  be  in  the  form: 

F(s)  -£  ak  sk  =  0  (3.1} 

k=0 

where  a,     =  b,  ct+c,  3  +d,  and  b,  ,  c,  and  d,  are  real 
constants.   Equation  (3-1)  can  also  be  rearranged  in  the 
following  general  form: 

f1(s)a  +f2(s) 3  +f3(s)  =  0  (3-2) 

where  f ,  (s) ,  f2(s)  an^  f-(s)  are  polynomials  of  s  of  any 
order.   Assuming  that  f-,  (s)  and  f2(s)  are  not  identical  to 
zero,  and  that  f,(s)  =  0  or  f?(s)  f    0,  yields  the  following 
two  cases: 

a.   f3(s)  f    0 

By  assuming  that  f3(s)  f    0,  dividing  equation 
(3-2)  by  f3(s)  yields: 

fT(s)     f2(s) 

-i a  +  — 3+1  =  0  (3-3) 

f3(s)      f3(s) 


or 
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F-.Cs)   a  +  F2(s)  B  +  l  m    0  (3-4) 

where  F,  (s)  and  F2(s)  are  either  a  ratio  of  two  polynomials 
of  s  or  a  polynomial  of  s. 

When  the  characteristic  equation  of  a  control 
system  is  considered  to  be  in  the  form  of  equation  (3-4), 
then  the  necessary  and  sufficient  conditions  for  this 
equation  to  be  singular,  according  to  equations  (2-15), 
are: 

/frjCs)  =  kxn 

(3-5) 

where  k,  and  k2  are  integers. 

It  must  be  noticed  that  according  to  the  process 
which  was  followed  in  the  derivation  of  equations  (2-15), 
they  are  valid  if  either  f-,(s),  f2(s)  and  f?(s)  are  con- 
sidered polynomials  of  s  or  ratios  of  two  polynomials  of  s. 
Based  on  this  remark  and  the  fact  that  the  argument  of 
unity  is  zero  or  an  even  multiple  of  n,  equations  (3-5) 
were  derived  from  (2-15). 

Considering  now  the  characteristic  equation  of  a 
control  system  in  the  form  of  equation  (3-4)  the  following 
cases  may  appear: 

(1)   Case  I:   F-L(s)  f    F2(s).   In  this  case  a 

point  s  of  the  s-plane  is  singular  for  the  charactereristic 
equation  (3-4),  if  and  only  if,  it  does  not  lie  on  the  real 
axis  and  satisfies  equations  (3-5). 
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In  order  to  determine  if  a  given  point  in  such  a 
case  is  singular  or  in  order  to  find  the  singular  points  of 
the  characteristic  equation  (3-4),  if  any,  one  of  the 
methods  which  were  described  in  the  previous  section  can 
be  used.   In  addition  to  these  methods  the  following  method 
which  is  related  to  the  Root  locus  concept  can  be  applied. 

The  characteristic  equation  of  the  control  system 
is: 

F1(.s)  a  +  F2(s)  g  +  i  =  0  (3-6) 

Suppose  that  the  Root  loci  of  the  following  equations  are 
drawn  on  the  same  drawing: 
A  F^s)  =  -1 

B  F2(s)  =  -1 

C  F1Cs)  =  -1 

D  F2(s)  =  -1 

where  A  and  B  are  Root  locus  variable  parameters  which  take 
values  from  zero  to  plus  infinity  and  C  and  D  are  Root 
locus  variable  parameters  which  take  values  from  zero  to 
minus  infinity. 

Based  on  the  angle  relation  on  which  each  Root 
locus  is  drawn,  all  the  points  of  each  of  the  above 
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Root  loci  will  satisfy  the  following  equations  respec- 
tively: 

/V-S)  =  C2n  ±  l)n 

/f2C.s)  =  (2n  ±  l)n 

/F1Cs)  =  2nn 

/F2(s)  =  2nn 

where  n  is  an  integer.   Therefore,  the  singular  points  of 
the  characteristic  equation  (3-6),  i.e.,  the  points  of  the 
s-plane  which  satisfy  equations  (3-5)  simultaneously,  will 
be  the  points  where  either  the  first  or  the  third  Root 
locus  intersect  with  the  second  or  the  fourth  Root  locus. 

From  the  above  described  method  it  can  be  con- 
cluded that  in  the  case  under  consideration  either  there 
are  no  singular  points  or  there  is  a  finite  number  of  such 
points.   Therefore,  when  in  general  F,  (s)  f    F?'-s-'  t*ie 
characteristic  equation  of  the  control  system  will  not 
always  be  singular  and  if  it  is  singular  the  probability  of 
finding  a  singular  point  which  meets  problem  specifications 
as  the  operating  point  of  the  system  will  be  very  low. 

(2)   Case  II:   F,(s)  =  F2(s).   In  this  case 
the  characteristic  equation  (3-4)  of  the  control  system 
becomes : 

Fx(s)  (a  +  B)  +1  =  0  (3-7) 
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Thus,  the  singular  conditions  given  by  (3-5)  become: 

/F1Cs)  -  kxn  (3-8) 

As  it  will  be  shown  in  the  following  case, 
when  there  is  only  one  angle  relation  to  be  satisfied, 
like  that  which  is  given  by  equation  C3-8) ,  then  there 
will  be  in  general  an  infinite  number  of  singular  points. 
Considering  equations  (2-5)  and  (2-16),  the  equation  of 
the  singular  line  of  any  singular  point  will  be  in  the 
general  form: 

a  +  6  =  c 
where  C  is  a  real  constant.   This  implies  that  although 
in  a  singular  case  there  will  be  an  infinite  number  of 
pairs  of  a  and  £  values  which  satisfy  equation  (3-7),  the 
roots  of  this  equation  will  be  fixed  for  any  of  these 
pairs.   For  this  reason  this  case  can  be  considered  as  a 
trivial  one  under  singular  considerations  and  not  useful 
for  applications. 

(3)   Case  III:   F2(s)  =  F1n(s).   When  F2(s)  = 

F,  (s) ,  where  n  is  an  integer,  then  the  characteristic 
equation  (3-4)  of  the  control  system  becomes: 

Fx(s)  a+F1n(s)  g  +  1  =  0  (3-9) 

and  the  singular  conditions  given  by  (.3-5)  become: 

/FjOO  =  k1n  (3-10) 
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i.e.,  in  this  case  a  point  s  of  the  s-plane  is  singular, 
if  and  only  if  it  does  not  lie  on  the  real  axis  and  satis- 
fies equation  (3-10) . 

In  order  to  determine  if  a  point  of  the 
s-plane  is  singular  for  the  characteristic  equation  (3-9) 
or  in  order  to  find  the  singular  points  of  this  equation 
one  of  the  methods  which  were  described  in  the  previous 
section  can  be  used.   In  addition  to  these,  the  Root  locus 
method  which  was  introduced  in  the  first  case  can  be  used 
as  follows: 

Suppose  that  the  Root  loci  of  the  following 
equations  are  drawn: 

A  F1(s)  =  -1 

B  FjCs)  =  -1 

where  A  and  B  are  Root  locus  variable  parameters  which 
take  values  from  zero  to  plus  infinity  and  from  zero  to 
minus  infinity,  respectively. 

Based  on  the  angle  relation  on  which  each 
Root  locus  is  drawn,  all  the  points  of  each  Root  locus 
will  satisfy  the  following  angle  relations  respectively: 
/FjCs)  =  (2k  +  l)n 

/FjCs)  =   2kn 

where  k  is  an  integer.   Therefore,  all  the  points  of  these 
two  Root  loci  will  also  satisfy  equation  (3-10) .   After 
these  remarks  it  can  be  concluded  that  in  the  case  under 
consideration  there  will  be,  in  general,  an  infinite 
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number  of  singular  points.   Also  the  probability  of  finding 
a  singular  point  which  meets  problem  specifications  as  the 
operating  point  of  the  system  will  be  relatively  high. 

In  the  simplest  case  where  n  =  2,  i.e., 

2 
F?(s)  =  F,  (s)  ,  the  characteristic  equation  [3-9)  becomes: 

F1(s)  a+F12(s)  3+1  =  0  C3-11) 

and  the  necessary  and  sufficient  conditions  for  this  equa- 
tion to  be  singular  will  be  given  again  by  equation: 

/Fj(s)  =  kn  (3-12) 

where  k  is  an  integer. 

(4)  Case  IV:   F1(s)  •   F2n(s).   In  this  case 

interchanging  the  two  parameters  a  and  3  in  the  character- 
istic equation  (3-4)  yields  the  previous  case. 


All  the  rest  of  the  cases,  which  may  appear  can 
be  considered  to  belong  in  Case  I,  and  can  be  analyzed  in 
a  similar  way. 

b.   f3(s)  =  0 

By  assuming  that  f,(s)  =  0,  equation  (3-2) 
becomes : 

f1(s)  <x  +  f2(s)  3  =*  Q  (3-13) 
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Since  f , (s)  and  f?(s)  are  not  equal  to  zero,  dividing  by 
f2(s)g  yields: 

-^ -+1  =  0  (3-14) 

f2Cs)   3 

or 

F,(s)  -  +  1  =  0  (3-15) 

i   g 

The  singular  conditions  for  this  case  as  it  can 
be  derived  from  equations  (3-5)  are  given  by  the  equation: 

/F-^s)  =  kn  (3-16) 

where  k  is  a  real  integer  or  zero.   It  can  be  shown  by  the 
same  procedure  as  was  done  in  Case  III  in  paragraph  2a(3), 
that  in  this  case  there  is  in  general  an  infinite  number  of 
singular  points.   Considering  equation  (2-5)  and  (2-16)  the 
equation  of  the  singular  line  of  any  singular  point  will  be 
in  the  form: 

-f-  =  C  (3-17) 

where  C  is  a  real  constant.   This  implies  that  although  in 
a  singular  case  there  will  be  an  infinite  number  of  pairs  of 
a  and  3  values  which  satisfy  equation  (3-15),  the  roots  of 
this  equation  will  be  fixed  for  any  of  these  pairs.   For 
this  reason  this  case  also  can  be  considered  as  a  trivial 
one  under  singular  considerations  and  not  useful  for 
applications . 
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3.    Summary  and  Conclusions 

In  the  preliminary  analysis  of  the  character- 
istic equation  of  a  control  system,  which  was  performed  in 
the  previous  paragraph,  it  was  found  that  there  is  a  form 
of  the  characteristic  equation,  which  in  general  has  an 
infinite  number  of  singular  points.   This  is  the  following: 

F,Cs)  a  +  F1n(s)  3  +  1  =  0  (3-18) 

or  in  the  simplest  case  where  n  =  2,  the  above  character- 
istic equation  becomes: 

F^s)  a  +  F12(s)  3  +  1  =  0  (3-19) 

where  F-,(s)  f    0,  and  it  is  either  a  polynomial  of  s  or  a 
ratio  of  two  polynomials  of  s. 

Therefore  equation  (3-18)  corresponds  to  the  most 
favorable  case  for  the  application  of  singular  line  theory 
in  control  system  design  and  compensation. 

The  next  step  is  to  find  the  structure  of  the 
system  (or  the  required  compensator  for  a  given  system) , 
which  has  a  characteristic  equation  corresponding  to  the 
form  of  equation  (3-18)  or  in  the  simplest  case  to  equation 
(3-19).   This  is  done  in  the  following  section. 
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IV.    ONE  INPUT  ONE  OUTPUT  LINEAR  FEEDBACK  SYSTEMS 

A.    PRELIMINARY  ANALYSIS 

In  figure  (4-1)  the  block  diagram  of  a  unity  feed- 
back control  system  is  shown.   Suppose  that  it  is  desired 
to  compensate  the  system  by  making  use  of  singular  line 
theory.   Then  according  to  the  results  of  the  previous 
section,  the  problem  is  to  find  the  compensation  structure 
or  structures  which  provide  a  system  with  a  characteristic 
equation  in  the  general  form: 

1  +a  F(s)  +g  F2(s)  =  0  (4-1) 

i.e.,  a  singular  system. 

Considering  a  compensation  structure  then  the  adjustable 
parameters  are  in  general  all  the  compensator  parameters, 
i.e.,  the  poles,  zeros  and  gain  values,  and  the  gain  of  the 
plant.   Since  the  poles  and  zeros  of  a  compensator  can  be 
either  real  or  complex,  a  preliminary  analysis  is  needed 
to  choose  the  two  parameters  a  and  6  from  among  the  gain 
values. 

In  the  following  sections  which  include  a  detailed 

analysis  of  the  singular  systems  it  is  shown  how  one  of 

these  parameters,  i.e.,  the  parameter  a,  was  introduced  in 

the  transfer  function  of  a  singular  cascade  compensator  in 

such  a  way  that  its  value  affects  its  poles.   (See  Section 
VI  C  and  equation  (6-18)). 
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In  the  preliminary  analysis  of  the  problem  of  finding 
the  compensation  structure  or  structures  which  provide  a 
system  with  characteristic  equation  in  the  form  of  equation 
(4-1),  i.e.,  a  singular  system,  the  five  most  commonly  used 
compensation  structures  were  analyzed  as  follows: 

1.    First  Compensation  Structure 

The  first  compensation  structure  which  was  con- 
sidered is  that  of  figure  (4-2). 

The  characteristic  equation  of  this  system  is: 
1  +  G(s)  +  G(s)  H(s)  =  0  (4-2) 

Comparison  of  equations  (4-1)  and  (4-2)  reveals 
that  in  order  to  obtain  a  singular  system  the  following 
relations  must  be  satisfied. 
G(s)  =  a  Fg(s) 

H(s)  =  ±   Fg(s)  C4"3) 

i.e.,  the  feedback  compensator  must  have  poles  and  zeros 
identical  with  those  of  the  plant,  and  the  two  parameters 
a  and  $  must  be  related  with  the  gain  values  as  is  shown  in 
equations  (4-3).   Introducing  equations  (4-3)  into  equation 
(4-2)  yields: 

1  +  aFg(s)  +BFg2(s)  =  0  (4-4) 

which  is  the  characteristic  equation  of  a  singular  system. 
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2 .    Second  Compensation  Structure 

The  second  compensation  structure  which  was 
considered  is  that  of  figure  (4-3) . 

The  characteristic  equation  of  this  system  is: 

1  +  GCs)  GcCs)  =  0  (4-5) 

Comparison  of  equations  (4-1)  and  (4-5)  reveals 
that  in  this  case  there  is  no  way  by  which  a  singular  system 
can  be  obtained  by  defining  the  two  parameters  a  and  3 
between  the  plant  and  cascade  compensator  gain  values. 

In  this  case  it  was  also  noticed  that  when  G(s) 
and  G  (s)  are  defined  as: 


(s  +  z  )  (s  +  z  )  ...  (s  +  z  ) 
G(s)  =  kFa(s)  =  k 


s  (s+pj)  (s+p2)  ...  (s+pn) 


and 


G  fs)  =  a+B  kFg(s)  = 


a   sN(s+p1) . . . (s+pn)+gk(s+z1) . . . (s+zm) 


g  "'      sN(s+Pl)  (s+p2) . . . (s+pn) 


the  characteristic  equation  of  the  system  becomes: 

1  +  akF  (s)+gk2F  2(s)  =  0 
S  ° 

Since  the  argument  of  the  second  and  third  term  in 
the  above  equation  are  not  affected  by  the  appearance  of  the 
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2 
k  and  k  factors  in  these  terms  respectively,  the  above 

equation  is  also  singular. 

Based  on  the  fact  that  the  cascade  compensator 

which  was  thus  introduced  in  order  for  a  singular  system  to 

be  obtained  has  the  same  number  of  poles  at  the  origin  as 

the  plant,  this  singular  case  (like  the  singular  case, 

which  was  analyzed  in  Section  V)  has  the  basic  disadvantage 

that  the  type  number  of  the  plant  is  not  preserved,  except 

in  the  case  where  the  plant  is  type  zero. 

3.   Third  Compensation  Structure 

The  third  compensation  structure  which  was  con- 
sidered is  shown  in  figure  (4-4)  . 
Suppose  that  in  general: 

G(s)  =  kFg(s) 
Gc(s)  =  kcFc(s) 
H(s)  =  kH  FH(s) 

The  characteristic  equation  of  this  system  will 


be 


1+G(s)  Gc(s)  +  G(s)  Gc(s)  H(s)  =  0  (4-6) 


Comparison  of  equations  (4-6)  and  (4-1)  reveals 
that  in  order  to  obtain  a  singular  system  the  following 
relations  must  be  satisfied: 

G(s)  Gc(s)  =  kkc  Fg(s)  Fc(s)  =  aFg(s)  Fc(s)    (4-7) 
and 
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H(s)  *  kH  FH(s)  =  |  FgCs)  FcCs)  (4-3) 

Introducing  equations  (4-7)  and  C4-8)  into  equa- 
tion (4-6)  yields: 

l  +  aFg(s)  Fc(s)  +  3Fg2(s)  Fc2(s)  =  0  (4-9) 

which  is  the  characteristic  equation  of  a  singular  system. 

This  case  can  be  considered  (under  singular  consid- 
erations) identical  with  the  first  one  of  figure  (4-2), 
when  G(s)  G  (s)  is  assumed  to  be  one  transfer  function. 

4.  Fourth  Compensation  Structure 

The  fourth  compensation  structure  which  was  con- 
sidered is  that  of  figure  (4-5). 

The  characteristic  equation  of  this  system  will  be 
in  general: 

1+G(s)  H(s)  +  G(s)  Gc(s)  =  0  (4-10) 

Comparison  of  equations  (4-1)  and  (4-10)  reveals 
that  for  a  given  plant  i.e.,  given  G(s)  =  kFCT(s),  there  is 
an  infinite  number  of  ways  by  which  G  (s)  and  H(s)  can  be 
defined  in  order  for  a  singular  system  to  be  obtained. 
Based  on  the  fact  that  certain  types  of  feedback  compensa- 
tion are  in  common  use,  e.g.,  velocity  and  acceleration 
feedback,  only  the  following  cases  were  selected  to  be 
analyzed  as  those  which  have  the  highest  probability  of 
application  in  practice. 
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a.   H(s)  =  a 

When  H(s)  is  a  gain  adjustment  which  is 
defined  to  be  the  parameter  a,  then  in  order  for  a  sin- 
gular system  to  be  obtained  the  following  should  be  the 
transfer  functions  of  figure  (4-5). 

G(s)  =  kFgCs) 

Gc(s)  =  gFg(s)  (4-11) 

H(s)  =  a 

Introducing  equations  (4-11)  into  (4-10) 


yields : 


l  +  akFs(s)  +  gkFg2(s)  =  0  (4-12) 


Because  the  factor  k  does  not  affect  any  argu 
ment  relation,  equation  (4-12)  will  correspond  to  a  sin- 
gular system. 

b,   H(s)  -    as 

When  H(s)  is  defined  to  be  a  tachometer  feed- 
back then  in  order  for  a  singular  system  to  be  obtained 
the  following  should  be  the  transfer  functions  of  the 
blocks  in  figure  (4-5) . 


GCs)  -  kFg(s) 

Gc(s)  =  Bs2Fg(s)  (4-13) 


H(s)  -  as 
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Introducing  equations  (4-13)  into  equation 
(4-10)  yields; 

1  +  aks  Fa(s)  +  gks2  F  2(s)  =  0  (4-14) 

which  is  the  characteristic  equation  of  a  singular  system. 

It  should  be  noted  that  a  basic  steady  state 
performance  requirement  of  the  system  shown  in  figure  (4-5) 
is  that  the  cascade  compensator  G _(s)  must  transmit  signal 
at  zero  frequency,  i.e.,  its  transfer  function  must  not 
have  any  zero  at  the  origin.   But  in  this  case  it  can  hap- 
pen, if  and  only  if,  F„(s)  has  at  least  two  poles  at  the 
origin,  which  implies  that  the  plant  must  be  at  least  type 
two . 

2 

c.   HCs)  =  as 

When  H(s)  is  defined  to  be  an  acceleration  feed- 
back, then  in  order  for  a  singular  system  to  be  obtained 
the  following  should  be  the  transfer  functions  of  the  blocks 
shown  in  figure  (4-5). 

G(s)  =  kFg(s) 

Gc(s)  =  3s4  Fg(s)  (4-15) 

2 
H(s)  =  as 


yields: 


Introducing  equation  (4-15)  into  equation  (4-10) 


1  +  aks2  Fg(s)  +  gks4  Fg2(s)  =  0  (4-16) 


which  is  the  characteristic  equation  of  a  singular  system. 
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It  should  be  also  noted  in  this  case  that  in 
order  for  the  cascade  compensator  G  (s)  to  transmit  signal 
at  zero  frequency  its  transfer  function  must  not  have  any 
zero  at  the  origin,  i.e.,  F  (s)  must  have  at  least  four 
poles  at  the  origin,  which  implies  that  the  plant  must  be 
at  least  type  four.   Because  such  a  plant  usually  cannot  be 
found  in  practice  this  case  will  not  be  further  considered. 


5.   Fifth  Compensation  Structure 

The  last  compensation  structure  which  was  con- 
sidered is  that  of  figure  (4-6) . 

By  considering  the  product  G(s)  G  2(s)  as  one 
transfer  function  then  under  singular  considerations  this 
case  turns  to  be  identical  with  the  previous  one. 

B .    SUMMARY 

In  this  section  a  preliminary  analysis  was  performed 
on  the  most  commonly  used  compensation  structures  in  order 
to  find  which  of  these  can  provide  a  singular  system,  and 
under  what  conditions,  i.e.,  a  system  with  characteristic 
equation  in  the  form  of  equation  (4-1). 

The  results  of  this  analysis  have  indicated  that  from 
a  theoretical  point  of  view  there  are  many  ways  by  which  a 
singular  system  can  be  obtained  by  properly  selecting  the 
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compensation  structure  and  by  properly  defining  the  com- 
pensators and  the  two  parameters  a  and  3 . 

The  conclusion  also  from  this  analysis  was  that  three 
singular  cases  must  be  considered  for  a  further  analysis, 
because  they  may  be  useful  and  applicable  in  control 
system  design  and  compensation.   These  are  the  following: 
The  first  singular  case  corresponds  to  the  compensation 
structure  of  figure  (4-2),  where  the  feedback  compensator 
and  the  two  parameters  a  and  3  are  defined  by  equation 
(4-3).   The  second  singular  case  corresponds  to  the  com- 
pensation structure  of  figure  (4-5) ,  where  the  compensators 
and  the  two  parameters  a  and  3  are  defined  by  equation 
(4-11) .   The  third  singular  case  corresponds  again  to  the 
compensation  structure  of  figure  (4-5),  where  the  com- 
pensators and  the  two  parameters  a   and  3  are  defined  by 
equation  (4-13) ,  where  the  plant  is  considered  to  be  at 
least  type  two. 

In  the  following  three  sections  which  refer  to  these 
three  singular  cases  respectively,  a  detailed  analysis  is 
performed. 
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V.    FIRST  SINGULAR  CASE  ANALYSIS 

In  the  previous  section  the  <<  first  singular  case  >> 
was  defined  as  that  which  corresponds  to  the  compensation 
structure  of  figure  (4-2) ,  where  the  associated  transfer 
functions  and  the  two  parameters  a  and  3  are  defined  as 

follows : 

(s+z,)    (s+z,)    ...    (s+zm) 
G(s)    =    aFg(s)    =    a-      LL± li . Ei-_  (5-1) 

sN(s+Pl)    (s+p2)    ...    (s+pn) 

and 

H(.J    -  i  Fg(s)    -   i      (S^''     (5*''}     •••    CS^  CS-2J 

01     sNfs+Pj)(s+pz)    ...    Cs+p„) 


a 


n' 


A.    SYSTEM  OUTPUT  AND  ERROR  ANALYSIS 


The  system's  transfer  function  is: 


T(.)    -   —  -    "V53  (5-3) 

R(s)         1    +aFg(s)    +3Fg2(s) 


or 


T(s)    =    s 


m  n 

a  Q  Cs+Zr3  P  ^2kj 


n  2       m_S_  -JL-  JL 


2NriL,,  ,2       N 

s 


J-n(s+pk)4asNp]Cs+Pk)p|(s  +  z;j)+Sp^(s  +  z.  f 
k=i  k=i  j  =i  j  =i 


(5-4) 
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The  output  of  the  system  is: 

C(s)  =  R(s)  T(s)  (5-5) 

and  the  error  of  the  system  is: 

E(s)  =  R(s)  {1  -  T(s)}  (5-6) 

When  the  final  value  theorem  was  applied  to  equations 
(5-5)  and  (5-6)  the  steady  state  output  and  error  were 
evaluated,  assuming  that  all  the  poles  of  the  transfer 
function  T(s)  were  on  the  left  half  of  the  s-plane. 

The  input  signals  considered  are  those,  which  are 
the  most  commonly  used,  i.e.,  step,  ramp  and  parabolic 
input . 

1.    Plant  Type  Zero 

In  this  case  where  the  plant  is  type  zero,  i.e., 

N 
N  =  0,  the  factor  s  which  appears  in  equation  (5-4)  be- 
comes unity. 

a.   Step  Input 

When  a  step  input  R(t)  =  Au(t)  or  R(s)  =  A/s 

was  considered  then: 

m  n 

Act         '     z  •     '     '     n, 

C  (t)    =   C      (t)    =    lim   s— T(s)    =   ^ J    k=1    3S 

t-»-cc  ss  s->o        s  n  m  n  m 

•  n  *  2+«n  z.  nPi  +enz  2 
k=i k    j=i  j  k=i k  j=i  j 
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and 

n  m 

Anp  2+u  n  z 2 

•  Ct)    =  essct)   -   lim  *A-  ll-TCsj}-  -fai    4,       Iff      1      = 

npk!.0n2nVBn,! 


b .      Ramp    Input 

When  a  ramp  input  RCt)  =  BtuCt) 

2 
or  R(s)  =  B  /  s  was  considered  then: 


COO  =  C   ft)  ■  lim  s "-i-  T(s)  = 
t^oc      "       s+o     sz 


and 


e(t)  =  ess(t)  =  lim  s  -i^Q  -  T(s)}  = 
t-*1  s-»-o     s 


c.   Parabolic  Input 


2 


When  a  parabolic  input  R(t)  =  ct  u(t) 


3 
or  R(s)  =  2c/s  was  considered  then: 


C(t)  =  C   ft)  =  lim  s  l£-  T(s) 
t-^  s+o   s 


and 

e(t)  =  e   ft)  =  lim  s  2c   {i-T(s)}  =  « 
t-*-  s+o    s-5 


2.    Plant  Type  One 

In  this  case  where  the  plant  was  considered  to 

be  type  one,  i.e.,  N  =  1,  the  steady  error  and  output  of 
the  system  were  found  as  follows: 


52 


a.   Step  Input 


C(t)  =  C   ft)  =  lim  s  £  T(s)  =  0 

t+oc        "  S  +  0     S 


and 


e(t)  =  e   ft)  =  lim  s  h   {l-T(s)}  =  A 

t">cc        S5         S  +  0     b 


b.   Ramp  Input 

m      n 

R        Bo.  L   .Z.   '  '  p. 

C(t)  =  C   (t)  =  lim  s  —  T(s)  =  1^ L^ — £- 

.KJ  ss  v  '  2  m 

t->-«  S  +  O     S^  1 r     2 

6  I  I  Z. 
j-i   J 


and 


e(t)  =  ess(t)  =  lim  s  L,  {l-T(s)}  =  <* 
t-*-«  s+o 

c.   Parabolic  Input 


C(t)  =  Css(t)  -  lim  s  ^f  T(s)  = 
t-*-*  s+o    s 


and 


e(t)  =  e   (t)  -  lim  s  ^£  {l-T(s)}  =  « 
t-*-^  s->o   s^ 


3.    Plant  Type  Two 

When  a  type  two  plant  was  considered,  i.e.,  N  =  2, 
the  steady  state  error  and  output  of  the  system  were  found 
as  follows: 
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a.   Step  Input 


C(t)  =  C   ft)  =  lim  s  £  TCs)  =  0 
t-»-«  s->o   s 


and 


e(t)  =  e   ft)  =  lim  s  |  (1-Tfs)}  =  A 

t->  cc  --  S->0 


b.   Ramp  Input 


C(t)  =  Css(t)  =  lim  s  -  T(s)  =  0 

t-K*  S  +  0     S 


and 


eft)  =  e   ft) 


=  lim  s  -9  {1-Tfs)}  =  oc 


c.   Parabolic  Input 


Cft)  =  C   (t)  -  lim  s  '■?!  T(s) 


m 


n 


2ca  I.  I  Z. 
122 1 


'k 


m 

1  =  1 


and 


eft) 

t+cc 


e   ft)  =  lim  s  ^  (1-Tfs)}  = 


ss 


s-*o   s 


4.    Plant  Type  Three  or  Over  Three 

In  all  such  cases  C   ft)  =  0,  and  the  error  be- 

ss  J  ' 

comes  finite  or  infinity  depending  on  the  input  signal 
which  is  considered. 
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5 .    Remarks 

From  the  above  results  which  are  summarized  in 
Table  I  of  Appendix  E,  it  can  be  seen  that  in  general  the 
type  of  the  plant  changes  due  to  the  form  of  the  feedback 
path  which  is  introduced  in  order  for  a  singular  system  to 
be  obtained.   Therefore,  compensation  of  a  given  system  by 
this  method  yields  a  system  which  in  general,  does  not 
respond  to  a  given  input  as  might  be  expected. 

Although  some  of  the  above  cases  which  were 
analyzed  may  have  a  specific  application  in  control  sys- 
tem design  and  compensation,  only  the  case  where  both 
output  and  error  have  finite  values  at  steady  state  seems 
to  be  of  a  more  general  interest. 

The  example  which  follows  is  referred  to  that 
case. 

B.    EXAMPLE  (5-1) 

Assume  that  in  figure  (4-2)  the  transfer  functions  of 
the  plant  and  the  feedback  compensator  are  as  follows: 

G(s)  =  aF  (s)  =  a — (5-7) 

g        (s+i)  (s  +  2)  (s  +  "0 

and 

H(s)  =  |  Fg  (s)  =  £  _^___ 


_  g     5+3 (5-8) 
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then  the  characteristic  equation  of  the  compensated  sys- 
tem is: 

1+  aFg(s)  +  6Fg2(s)  =  0  (5-9) 

or  after  substituting  equations  (5-7)  and  (5-8)  into  (5-9) 
yields : 

s6  +1455  +(77+a)s4  +(212  +  10a)s3  +(308  +35a+3)  s2  + 

+  (224  +  50a  +  66)  s  +  (64  +  24a+9g)  =  0  (5-10) 

By  using  the  method  of  paragraph  6c,  section  II,  the 
singular  points  of  the  system  were  found  for  several  values 
of  £.   They  are  listed  in  Appendix  A.   One  of  these  sin- 
gular points,  i.e.,  the  singular  point  (t   =  0.60,  to   = 
3.1104)  was  further  considered.   Then  from  equation  (2-6) 
the  Chebyshev  functions  of  the  second  kind  were  evaluated 
and  from  equation  (2-5)  the  coefficients  B, ,  C,  and  D, 
were  found  to  be: 


B,  =  -5.3085 
G,  =  0.6746 
Dx  =  41.7741 


Then  from  equation  (2-16)  the  equation  of  the  correspon- 
ding singular  line  was  found  to  be: 
3  =  -  61.9254  +  7.8693a  (5-11) 
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With  a  computer  program  then,  the  roots  of  the 
characteristic  equation  (5-10)  were  found  for  several 
different  parameter  pair  values,  which  satisfy  equation 
(5-11).   They  are  listed  in  Appendix  B.   By  examining 
these  roots  it  can  be  seen  that  the  singular  root  and  a 
real  root  are  fixed  and  the  other  three  roots  are  varying 
for  different  pair  values. 

In  figure  (5-1)  the  parameter  plane  diagram  of  this 

singular  system  is  shown  for  some  arbitrarily  selected 

values  of  C  and  to  . 

n 

In  figure  (5-2)  the  parameter  plane  diagram  of  the 
same  singular  system  is  shown  for  one  constant  c  and  one 
constant  u  curve,  which  correspond  to  the  singular  values 
of  the  previously  considered  singular  point,  i.e.,  c   = 
0.60  and  w   =  3.1104.   On  the  same  diagram  the  singular 
line  of  this  singular  point  is  also  shown.   It  was  noticed 
that  these  two  curves  and  the  singular  line  are  tangent  to 
each  other  on  the  same  point  of  the  parameter  plane. 

For  a  further  analysis  of  this  system  and  of  any 
other  system  which  is  referred  to  the  <<  first  singular 
case  >>,  the  same  procedure  which  is  described  in  the  next 
section  and  which  is  referred  to  the  "second  singular  case", 
can  be  followed. 
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VI.    SECOND  SINGULAR  CASE  ANALYSIS 

In  Section  IV,  the  "second  singular  case"  was  defined 
as  that  which  corresponds  to  the  compensation  structure  of 
figure  (6-1)  where  the  associated  transfer  functions  and 
the  two  parameters  a  and  3  are  defined  as  follows: 


(s  +  ZJCs  +  ZJ  ...  (s  +  Z  ) 

G(s)  =  kF  (s)  =  k  -Tr— - 5L_  (6-1) 

s  (s+p)(s+p2)  ...  (s+pn) 


(s  +  Z^Cs  +  Z  )     ...  (s  +  Z  ) 

Gfs)  =  3F_Cs)  -  3  -jj— - 2—         (6-2) 

g        sw(s+Pl)(s+p2)  ...  (s+pn) 


H(s)  =  a  (6-3) 

A.    SYSTEM  OUTPUT  AND  ERROR  ANALYSIS 

The  transfer  function  of  the  system  is: 

R(-SJ    1  +  G(s)  H(s)  +  G(s)Gc(s) 

Introducing  equations  (6-1)  through  (6-3)  into  equa- 
tion (6-4)  yields: 
m 

n  rs+7  .« 


3k  i.  I  (s  +  Z.) 

T(s)  =  — 1^_: — r 


?N    n  ,       M  n  m  m  2 

s   n(s+pk}  +aks  ncs+p^ncs+z^+Bkpics+ZjD2 

k=i  k=i  j=i  j-i 

(6-5) 
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The  output  of  the  system  is: 

C(s)  =  R(s)  T(s)  (6-6) 

and  the  error  of  the  system  is: 

E(s)  =  R(s)  -  C(s)  =  R(s)  U-T(s)}  (6-7) 

When  the  final  value  theorem  was  applied  to  equations 
(6-6)  and  (6-7)  the  steady  state  output  and  error  values 
were  evaluated,  assuming  that  all  the  poles  of  the  trans- 
fer function  (6-5)  lie  on  the  left  half  plane.   This 
analysis  was  based  again  on  the  three  usually  considered 
input  signals,  i.e.,  step,  ramp  and  parabolic  input. 

1.   Plant  Type  Zero 

In  this  case  where  the  plant  was  type  zero,  i.e., 
N  =  0,  the  steady  state  error  and  output  of  the  system 
were  as  follows: 

a.   Step  Input 

When  a  step  input  R(t)  =  Au(t)  or  R(s)  =  A/s 
was  considered  then: 


m 
Agk    I    \z 


n  2 


A  TfO    =  J=1    3 


C(t)    =   C      (t)    =    lims-T(s)    =   ., 

^v    '  ss      J  s      *■    J        n  n  m  m 

I     |Pk  +akl     lp,  I     lz.    +gkl      lz.2 
k=1  k-i  >i  J         j  =  i  J 
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and 


aCt)  =  ess( 

t-*oc 


n        n     m 

•  rn  2   ~~ i  n  ) 

a         AvJ  ,'p,  +ak,  ,  p,  I.  »z .  J 

t)  =  lim  sA{i-T(s)}=-i=kzirk ^\l"    JJ 

Q  n        n     m       m 

s-»-o   s         l — r      t — r   t — r     t — r 


k=  i 


Pk 


+  ak|  |  p.  I  |z.+ 
k=i   Kj=i  J 


m       m 

3  k 


J 


z  . 
=  i  ^ 


b.   Ramp  Input 


When  a  ramp  input  R(t)  =  Btu(t) 

2 
or  R(s)  =  B/s   was  considered  then: 


C(t)  =  C   ft)  =  lim  s  — 2  T(s)  = 

t">«        '  S  +  O     S 


and 


e(t)  =  e   (t) 


B 


=  lim  s  =-2    (l-T(s)}  =  « 


s+o 


c.   Parabolic  Input 


When  a  parabolic  input  R(t)  =  ct  u(t) 

3 
or  R(s)  =  2c  /  s  was  considered  then: 


C(t)  =  C   ft)  =  lim  s  ^|  T(s) 
t-*a  "'  s->o   s 


and 


e(t) 

t+<= 


e   (t) 

ss  v  J 


lim  s  ^|  {l-T(s)}  =  = 
s-*-o   s 


60 


2.    Plant  Type  One 

In  this  case  where  the  plant  was  considered  to 
be  type  one,  i.e.,  N  =  1,  the  steady  state  error  and  out- 
put of  the  system  were  found  to  be  as  follows : 

a.   Step  Input 

When  a  step  input  was  considered  then: 

C(t)  =  C   ft)  =  lim  s  |  T(s)  =  A 
t-**     "  ~J  s->-o 


and 


e(t)  =  ess(t)  =  lim  s  |  U-T(s)}  =  0 
t-**     "      s->-o 


b .   Ramp  Input 


When  a  ramp  input  was  considered  then 


C(t)    =    Css(t)    =    lim   s    ^    T(s)    = 
t-*x  ^  s-*-o        s 


and  n 

v>  Ba    '     '    Pi 


(t)    =   ess(t)    =    lim   s   -2    {1-T(s)}  = 


t-**  S->o  s 


m 


z  . 
3 


i  =  i 


c.   Parabolic  Input 

When  a  parabolic  input  was  considered  then: 
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Get)  =  C   ft)  =  lim  s  ^   T(s)  = 
t-»-<*  s-*-o   s 


and 


e(t)  =  e   ft)  =  lim  s  H  'U-T(s)}  =  * 
t-*-<*  s° 


3.    Plant  Type  Two 

In  this  case  where  the  plant  was  considered  to 
be  type  two,  i.e.,  N  =  2  the  steady  state  error  and  output 
of  the  system  were  evaluated  as  follows: 

a.   Step  Input 

When  a  step  input  was  considered  then: 

C(t)  =  Css(t)  -  lim  s  |  T(s)  =  A 

t+cc        """  S  +  O 


and 


e(t)  =  ess(t)  =  lim  s  |  U-T(s)}  =  0 
t-*=  s+o 


b .   Ramp  Input 


When  a  ramp  input  was  considered  then: 


C(t)  =  C   ft)  =  lim  s  \   T(s)  = 

t-**        '  S  +  O     S 


and 


e(t)  =  e   ft)  =  lim  s  \    {l-T(s)}  =  0 
t^     ss      s+o   s 
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c.   Parabolic  Input 


When  a  parabolic  input  was  considered  then: 


C(t)  =  e   ft)  =  lim  s  H   T(s)  = 

t+cc        bb         S  +  O     S 


and  n 

.  -  2ca  J.  J  pv 

eft)  =  e   ft)  ■  lim  s  -±§  {l-T(s)}  =  Jc"1   k 

t+cc     -"■■                          s  3  jh    z 


4 .    Remarks 

As  a  first  step  toward  a  further  study  of  the 
"second  singular  case",  an  error  and  output  analysis  was 
performed.   From  these  results  which  were  obtained  and 
which  are  summarized  in  Table  II  of  Appendix  E,  it  was 
noticed  that  in  contradiction  with  the  "first  singular 
case"  the  type  of  the  plant  is  preserved  in  this  case, 
where  the  compensation  structure  and  the  associated  com- 
pensators of  figure  (6-1)  are  introduced  in  order  for  a 
singular  system  to  be  obtained. 

B.    EXAMPLE  (6-1) 

In  this  example  a  type  one  plant  was  considered,  with 
transfer  function. 

G(s)  =  kF  (s)  =  100     'L  (6"8) 


g  s(s  +  l)(s  +  5) 
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Suppose  that  this  plant  is  compensated  by  the  com- 
pensation structure  shown  in  figure  (6-1) ,  then  according 
to  equations  (6-2)  and  (6-3)  the  transfer  functions  of  the 
compensators  are: 


Gfs)  =  3  F  (s)  =  3 (6-9) 

c         g        s(s+l)(s+5) 


and 

H(s)  =  a  (6-10) 

The  characteristic  equation  of  the  compensated  system 
is : 

1  +  G(s)  H(s)  +  G(s)  Gc(s)  =  0  (6-11) 

Introducing  equations  (6-8),  (6-9)  and  (6-10)  into  equation 
(6-11)  yields: 

s6  +  12s5  +  46s4  +(60+100a)s3  +  (25+600a)s2 

+500  as+lOOg  =  0  (6-12) 

Using  the  computer  method  of  paragraph  6c,  Section  II, 
the  singular  points  of  the  system  were  found  for  several 
values  of  r, .      They  are  listed  in  Appendix  C. 

One  of  these  singular  points  was  further  considered, 
i.e.,  the  singular  point  (?  =  0.40,  w   =  0.97096).   Then, 
corresponding  to  the  damping  ratio  value,  Chebyshev 
functions  of  the  second  kind  were  evaluated  from  equation 
(2-6)  and  from  equations  (2-5)  and  (2-16)  the  equation  of 
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the  singular  line  which,  corresponds  to  the  singular  point 
which  was  considered,  was  evaluated  to  be: 

3  =  -  4.92497a  +  0.25924  (6-13) 

Considering  the  last  two  terms  of  equation  (6-12)  a 
necessary  condition  for  stability  according  to  Routh's 
first  criterion  is  that  a>0  and  $>0.   Since  a  and  3 
parameters  are  related  by  equation  (6-13)  ,  it  implies  that 
a  necessary  condition  for  stability  is: 

0<  a  <    0.05265  (6-14) 

By  arbitrarily  selecting  a  =  0.05  from  the  range  of 
values  defined  by  (6-14),  the  corresponding  value  of  $ 
was  found  from  equation  (6-13)  to  be,  g  =  0.013.   Then  the 
roots  of  the  characteristic  equation  (6-12)  were  evaluated, 
for  this  pair  of  parameter  values,  by  using  a  computer 
program  with  the  ZPOLR  subroutine.   These  roots  are: 

s1  =  -  0.3884  +  jO .8899  s2  =  -  0.3884  -  J0.8899 

s2  =  -  5.2141  s4  =  -  5.0132 

s5  =  -  0.9270  s6  =  -  0.0701 

i.e.,  one  of  the  roots  is  the  singular  root  (S  =  0.40,w  = 


ns 


0.97096)  or  s  =  -  0.3884  ±  j  0 .8899 
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In  figure  (6-2)  the  parameter  plane  diagram  of  this 
singular  system  is  shown  for  some  arbitrarily  selected 
values  of  z,   and  w 


n 


C.    HARDWARE  CONSIDERATIONS 

In  figure  (6-3)  the  block  diagram  of  a  unity  feed- 
back system  (plant)  is  shown.   The  open  loop  transfer 
function  of  this  system  is: 

Cs+.Zj  )(s  +  z2)   ,.,.,..    (s  +  z    )  f6-151 

G(s)    =   kF    (s)    =   k  -^ 2 (0    "J 

s  s1N(s+Pl)(s+p2)    ...    (s+pn) 

Compensation  of  this  system  according  to  the  "second 
singular  case"  scheme,  yields  the  compensated  singular 
system  of  figure  (6-4) .   This  compensation  structure  will 
be  in  general  complex  since  it  introduces  a  cascade  and 
a  feedback  compensator.   Therefore,  a  simplified  version 
may  also  be  desired.   This  modification  is  shown  in  figure 
(6-5),  where  the  feedback  compensator  (gain  adjustment) 
has  been  taken  from  the  plant  and  has  been  introduced  at 
the  cascade  compensator  G _(s) .   Also  the  Root  locus  gain 
3  from  G  (s)  has  been  set  as  a  gain  adjustment  in  front  of 
the  minor  loop. 

The  characteristic  equation  of  the  systems  in  figure 
(6-5)  and  (6-6)  are  respectively: 

1  +  ak  F  (s)  +  3k  F  2(s)  =  0  (6-16) 

o  6 
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and 

1  +  a  Fg(s)  +  3k  Fg2(s)  =  0  (6-17) 

Comparison  of  these  equations  reveals  that  both  are 
singular  and  they  have  the  same  singular  points,  but  they 
will  be  in  general  different,  except  when  k  =  1,  in  which 
case  they  become  identical.   By  this  modification  the 
whole  compensation  structure  has  been  introduced  in  cas- 
cade with  the  plant,  as  is  shown  in  figure  (6-5). 

An  equivalent  system  of  figure  (6-5)  is  shown  in 
figure  (6-6)  where  the  whole  compensation  structure  of 
figure  (6-5)  is  replaced  by  an  equivalent  cascade  com- 
pensator. 

The  following  considerations  can  now  be  stated  con- 
cerning the  compensation  structures  of  figures  (6-4)  and 
(6-5)  or  (6-6).   First  both  are  singular.   Second  with 
the  compensation  structure  of  figure  (6-5)  or  (6-6),  the 
compensation  problem  has  been  simplified  and  better 
organized  since  it  will  be  required  to  build  a  filter  from 
distinct  components  or  a  filter  in  an  integrated  form 
which  will  be  introduced  in  cascade  with  the  plant.   This 
is  a  simpler  and  in  general  a  more  realistic  approach  to 
the  problem.   In  the  case  also  where  K  =  1,  i.e.,  the  gain 
of  the  plant  is  equal  to  unity  then  these  structures  are 
also  equivalent.   When  K  f   1  then,  although  they  will  have 
the  same  singular  points  (because  the  singular  conditions 
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to  be  satisfied  are  identical,  in  both  cases),  the 
corresponding  equations  of  the  singular  lines  will  be 
different.   This  means  that  for  the  same  singular  point, 
different  points  of  the  parameter  plane  (a,  B) ,  will  be 
associated  and  also  the  non-singular  characteristic  roots, 
will  be  in  general  different. 

The  cascade  compensator  G   (s)  which  was  introduced 
in  figure  (6-6)  has  the  transfer  function: 

3  (s  +  *1.).C*+*20  •••  (s  +  z  )  (6"18) 

G   (s)  =  -j, 2 

s   (s+PiHs+pz)  ...  (s+pn)  +a(s  +  Zl)  ...  (s  +  zm) 

where  the  subscripts  s  c  stands  for  <<  Singular  compen- 
sator >>.   Assuming  that  a  pair  of  a  and  $  parameter  pair 
values  which  satisfies  problem  specifications  has  been 
selected,  then  this  compensator  and  its  transfer  function 
can  be  precisely  determined. 

This  modified  singular  compensation  structure  of 
figure  (6-5)  or  (6-6)  was  further  analyzed  as  the  most 
applicable  due  to  its  simplicity  under  hardware  considera- 
tions, in  comparison  with  that  of  figure  (6-4).   It  must 

r 

be  also  noticed  that  in  a  specific  application,  where  a 
gain  feedback  adjustment  at  the  plant  is  possible  or 
desired,  e.g.,  when  it  is  desired  to  eliminate  the  effect 
of  some  nonlinearities  that  the  plant  may  have,  then  the 
initial  compensation  structure  of  figure  (6-4)  may  be 
considered  for  application. 
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D.    STABILITY  ANALYSIS 


1.    Introduction 


One  of  the  basic  requirements  of  any  control 
system  is  stability.   The  stability  concept  and  related 
definitions  can  be  found  in  reference  [7 ] .   The  first 
question  which  arises  when  stability  considerations  are 
made  for  the  singular  structure  of  figure  (6-6)  is:   What 
effect  has  the  cascade  singular  compensator  which  has 
been  introduced,  on  the  stability  of  the  initial  uncom- 
pensated system?   That  is,  in  addition  to  the  fact  that 
this  compensator  makes  the  system  singular,  does  it 
improve  or  impair  the  stability  of  the  initial  system? 
In  either  case  the  question  which  seeks  an  answer  is:   How 
is  it  possible  from  an  initially  stable  or  unstable 
system  to  obtain  a  stable  singular  system?   The  interpre- 
tations of  curves  and  singular  lines  on  a  parameter  plane 
diagram  and  determination  of  stability  is  not  in  general 
a  simple  process. 

Stability  analysis  by  using  the  Root  locus 
concept  seems  to  be  relatively  easy  and  can  give  a  more 
complete  picture  of  the  behavior  of  the  system.   The 
stability  analysis  as  it  can  be  performed  by  using  the 
Root  locus  concept  for  the  case  of  the  compensated  singular 
system  is  described  below. 
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2 .    Stability  Analysis  by  the  Root  Locus  Concept 

By  substituting  F  (s)  from  equation  (6-15)  into 
equation  (6-17)  and  after  some  manipulations  the  character 
istic  equation  of  the  compensated  singular  system  becomes: 


a  (s+z1)(s+z2)  ...(s+zm)s  (s+p!((s+p2)  . . . (s+pn) +3k(s+Z! )  . 
(s+zm)2+  s2N  (s+Pl)2  (s+p2)2  ...(s+pn)2  =  0         (6_19) 


Assuming  that  the  singular  points  of  the  singular 
compensated  system  have  been  found  and  that  one  of  them 
has  been  selected  as  the  operating  point  of  the  system, 
then  the  equation  of  the  corresponding  singular  line  can 
be  evaluated  from  equations  (2-5)  and  (2-16)  .   This 
equation  will  be  in  general  in  the  form: 

3  =  Aa+B 

where  A  and  B  are  real  constants.   By  substituting  the 
parameter  3  from  this  equation  (6-19)  and  after  some 
manipulations,  the  characteristic  equation  of  the  singular 
compensated  system  becomes: 

2  2  2 

AkCs+Zj)    (s+z2)    ...(s+zj   +sN(s+p1)(s+p2)...(s+pn)(s+z1)(s+z2)...(s+zin) 

a  ^  ^  2      2M  ~2  1  "  ~2  :"1 

Bk(s+Z!)    (s+z2)    ...(s+zm)   +s     Cs+Pi)    (s+P21    ...(s+pn) 

(6-20) 
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Since  A,  B  and  k  are  known  constants  the  numerator 
and  denominator  of  the  above  equation  can  be  factored, 
which  in  general  yields: 


Cs+Zj'  )(s  +  z^  ).......  (s  +  z-  ) 

a 2 —  =  -  i        (6-21) 


(s+p{  )(s+p2*  )  ...  (s+p^  ) 


where  z '.   and  p.!   may  be  real  or  complex.   By  treating  the 
parameter  a  as  the  Root  locus  variable  parameter  the  Root 
locus  of  the  compensated  system  can  be  drawn,  either  by 
the  classical  graphical  methods  or  by  the  use  of  a  com- 
puter.  In  the  last  case  where  a  computer  is  used,  usually 
it  is  possible  to  proceed  directly  from  equation  (6-19)  to 
draw  the  Root  locus  of  the  compensated  singular  system. 
Since  it  is  known  that  as  the  Root  locus  variable 
parameter  a  of  equation  (6-21)  varies  from  zero  to  infinity 
the  roots  of  the  characteristic  equation  move  on  the  Root 
locus  segments,  starting  from  the  poles  and  terminating 
to  zeros  or  at  infinity  when  zeros  are  not  available,  and 
that  the  singular  root  and  maybe  one  or  more  other  roots 
of  the  characteristic  equation  (6-21)  are  fixed,  indepen- 
dent of  the  variations  of  the  parameter  a,  it  implies  that 
some  zeros  and  poles  of  equation  (6-21)  have  to  be 
identical.   Therefore,  the  numerator  and  denominator  of 
equation  (6-21)  must  consist  of  a  number  of  common  factors 
equal  to  the  number  of  the  fixed  roots  that  the  charac- 
teristic equation  of  the  singular  compensated  system  has. 
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This  number  is  at  least  two,  corresponding  to  the 
system's  complex  singular  root.   This  fact  makes  the  pro- 
cess of  drawing  the  Root  locus  of  the  compensated 
singular  system  by  the  classical  graphical  methods,  much 
simpler  than  it  initially  looks.   Thus  the  Root  locus  of 
the  singular  compensated  system  consists  of  two  or  more 
fixed  points  of  the  s-plane  plus  the  classical  Root 
locus  segments. 

Considering  now  system  stability  all  the  related 
questions  can  be  answered  from  the  Root  locus  study.   For 
example,  if  any  Root  locus  segment  or  isolated  point  lies 
entirely  on  the  right  half  of  the  s-plane,  then  there  is 
no  pair  of  a  and  £  parameter  values  for  which  the  compen- 
sated system  can  be  stable.   If  on  the  other  hand  all  the 
Root  locus  segments  lie  on  the  left  half  of  the  s-plane 
as  well  as  all  the  isolated  points  then  all  pairs  of  a 
and  3  parameter  values  will  yield  a  stable  system.   Finally, 
if  some  Root  locus  segments  cross  the  imaginary  axis  in  one 
or  more  points,  then  the  range  of  the  parameter  a  values 
for  which  the  compensated  singular  system  is  stable  can  be 
determined.   The  corresponding  values  of  the  parameter  3 
are  then  found  from  the  equation  of  the  singular  line. 

The  following  example  is  an  application  of  the  above 
procedure  and  concepts  which  were  introduced  for  the  case 
of  the  singular  compensated  system. 
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3.    Example  (6-2) 

In  this  example  the  uncompensated  system  (Plant) 
of  figure  (6-3),  was  considered  to  have  transfer  function: 

G(s)  =  kF  (s)  =  k  1 ^6"22) 

g        s(s+l)  (s+5) 

The  characteristic  equation  of  this  system  is: 

s3  +  6s2  +  5s  +  k  =  0  (6-23) 

The  Root  locus  of  this  system  is  shown  in  figure  (6-7). 
The  gain  at  the  stability  limit  is  k=30.   When  this  system 
is  compensated  by  a  singular  cascade  compensator  the  block 
diagram  of  the  compensated  system  will  be  that  which  is 
shown  in  figure  (6-5)  or  the  equivalent  of  figure  (6-6)  . 
The  singular  cascade  compensator  which  was  thus  introduced 
has  a  transfer  function: 

s  +6s  +5s+a 

which  was  obtained  from  equation  (6-18)  . 

The  characteristic  equation  of  the  compensated 
singular  system  is: 


1+G(s)  Gsc(s)  =  0  (6-25) 
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Introducing  equations  (6-22)  and  (6-24)  into  equation 
(6-25)  yields: 

s6+12s5+46s4+(60+a)s3+(25+6ct)s2+5as+k3  =  0  (6-26) 

Two  values  of  the  gain k were  arbitrarily  selected, 
i.e.,  k =20  and k =35,  by  which  the  uncompensated  system 
becomes  stable  and  unstable  respectively.   These  cases 
were  separately  analyzed  as  follows: 

a.   Uncompensated  System  Stable,  k  =  20 

For  this  value  of  k,  the  uncompensated 
system  is  stable  and  the  characteristic  equation  (6-26) 
of  the  compensated  singular  system  becomes: 

s6+12s5+46s4+(60+a)s3+(25+6a)s2+5as+20g  =  0  (6-27) 

In  figure  (6-8),  the  parameter  plane  diagram 
of  the  compensated  singular  system  is  shown  for  arbitrarily 
selected  values  of ?  and  u^  .   By  choosing  cs=0.40  the 
corresponding  value  of  to   was  found  in  the  computer  by  using 
the   <<  Singular  Point  Program>>  as  cjns  =  0.97096.   Therefore, 
a  singular  root  of  the  compensated  singular  system  is: 


s  =  "  Cs  ^ns  ±   S"ns/  1-Cs  =  "  0.3384  ±  j  0.8899  T6-28) 

From  equations  (2-5)  and  2-16),  the  equation  of  the 
corresponding  singular  line  was  found  to  be: 

3  =  0.24621a  -  1.21242  (6-29) 

i.e.,  if  the  pair  of  a  and g  parameter  values  satisfies 
equation  (6-29),  then  the  characteristic  equation  (6-27) 
of  the  compensated  singular  system  will  have  a  fixed  complex 
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root  at  the  point  s  ■  -  0.3884  +j  0.8899. 

Considering  now  the  coefficients  of  the  zero 
and  first  power  of  s  terms  of  the  characteristic  equation 
(6-27),  a  necessary  condition  for  stability  is  that  both 
these  are  greater  than  zero,  which  implies  that  <*>  o  and 
g>0.   Based  on  this  observation  the  roots  of  the  character- 
istic equation  (6-27)  were  found  for  a  specific  range  of  the 
parameter  a  and  3  values  and  with  each  parameter  pair 
satisfying  equation  (6-29),  by  using  the  subroutine  ZPOLR 
in  a  computer  program.   These  roots  are  listed  in  Appendix 
D.   From  this  list  several  general  remarks  can  be  made 
about  the  range  of  the  a  and  6  parameters  in  the  stable 
region  and  the  variation  of  the  nonsingular  roots  in  this 
region,  i.e.,  the  stable  region  in  terms  of  the  parameter 
a  value  is  bounded  between  the  values  of  5  and  40  approx- 
imately, and  that  within  this  region  the  variations  of 
the  nonsingular  roots  is  such  that  the  dominance  and 
sensitivity  concepts  are  considerably  affected,  for  dif- 
ferent parameter  values. 

Substituting  the  value  of  3  from  equation  (6-29) 
into  equation  (6-27)  the  characteristic  equation  of  the 
compensated  singular  system  becomes: 


s  +12s5+46s4+(60+a)s3+(25+6a)s2+5as+(4.9242a  -24.2485)  =0      (6-3  0) 


Treating  then  the  parameter  a  as  the  Root  locus  variable 


75 


parameter,  equation  (6-30)  can  be  written  as: 


a 


s5  +  6s2  +  5s+  4.9242 =  -  _  (6-31) 

s6+12s5+46s4+60s3+25s2  -24.2485 


or  by  factoring  numerator  and  denominator  yields: 


(s+5. 2230) (s2+0. 7763s  +  0.9428) 


a 


(s+5. 2230) (s2+0.7768s+0.9428) (s+4.7195) (s+1.8458) (s-0.5653) 


(6-32) 

The  numerator  and  denominator  of  the  left  part 
of  equation  (6-32)  have  common  factors  corresponding  to 
the  real  root  s  =  -5.2230  and  the  complex  singular  root 
s  =  -0.3884  +  JO. 8899.   Therefore,  the  characteristic 
equation  of  the  compensated  singular  system  has  a  fixed 
real  root  and  a  fixed  complex  root  at  these  points  of 
the  s-plane. 

The  Root  locus  of  this  compensated  singular 
system,  which  was  obtained  from  equation  (6-30).   By  the 
use  of  the  R00TL0  subroutine  in  a  computer  program  is 
shown  in  figure  (6-9). 

At  this  point  it  must  be  noticed  that,  the 
remarks  stated  in  paragraph  two  about  the  fixed  roots 
that  a  singular  system  has,  the  analytical  results  shown 
in  Appendix  D,  and  the  Root  locus  of  figure  (6-9)  are 
in  agreement.   It  must  also  be  noticed  that  the  fixed 
roots  of  the  characteristic  equation  (6-30)  are  shown  by 
dark  cross  marks  in  a  circle,  in  figure  (6-9) . 
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b.   Uncompensated  System  Unstable, k  =  35 

When  the  gain  of  the  uncompensated  system 
(Plant)  has  the  value  of  35,  then  the  uncompensated  system 
is  unstable  since  a  pair  of  complex  roots  of  the  character- 
istic equation  (6-23)  lie  on  the  right  half  of  the  s,-plane. 

In  this  case  the  characteristic  equation  (6-26) 
of  the  compensated  singular  system  shown  in  figure  (6-6) 
becomes: 

s6+12s5+46s4+(60+a)s3+(25+6a)s2+5as+353  -  0  (6-33) 

In  figure  (6-10)  the  parameter  plane  diagram 

of  this  compensated  singular  system  is  shown. 

For  the  same  value  of  r  as  in  the  previous 

s 

case,  i.e.,?   =  0.40,  the  corresponding  value  of  w   was 
found  by  using  the   "Singular  Point  Program"  as  u)  =0.97096. 
Therefore,  a  singular  root  of  the  compensated  singular 
system  is  s  =  -0.3884  -  j 0.8899.   The  fact  that  the  singular 
root  remains  the  same  as  in  the  previous  case,  although  the 
characteristic  equation  has  been  changed  comes  as  a  conse- 
quence of  the  singular  line  theory  which  was  developed  in 
section  II  and  III,  i.e.,  if  the  characteristic  equation 
of  the  compensated  singular  system  is  considered  to  be  in 
the  form  of  equation  (6-17),  then  in  both  cases  Fg  (s)  is 
the  same  and  only  the  value  of  k  has  been  changed.   Because 
the  singular  points  according  to  the  singular  line  theory 
are  the  points  of  the  s-plane,  which  satisfy  equation: 

/Fg(s)  =  n  n 
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where  n  is  a  real  integer  or  zero,  it  implies  that  since 
Fcr(s)  is  the  same  in  both  cases  the  singular  points  will 
remain  the  same. 

From  equations  (2-5)  and  (2-16)  the  equation  of 
the  corresponding  singular  line  was  found  to  be: 

3  =  0.14069a  -  0.69281  (6-34) 

For  the  same  values  of  the  parameter  a,  as  those 
which  were  used  in  the  previous  case,  the  corresponding 
values  of  the  parameter  g were  found  from  equation  (6-34). 
For  these  formed  pairs  the  roots  of  the  characteristic 
equation  (6-33)  were  evaluated  by  using  the  subroutine 
ZPOLR  in  a  computer  program.   They  are  also  listed  in 
Appendix  D.   By  comparing  the  roots  of  the  characteristic 
equation  (6-27)  and  (6-23),  which  are  listed  in  Appendix  D, 
it  can  be  seen  that  they  are  identical  for  the  same  value 
of  the  parameter  a.   The  reason  for  this  fact  is  that  when 
the  singular  point  s  =  -0.3S84  +  j 0.8899  or  (C   =  0.40, 

co    =  0.97096),  is  considered,  the  equation  of  the 

ns  '  *  >       n 

corresponding  singular  line  as  a  function  of  k,  according 
to  equation  (2-5)  and  (2-16)  is: 

6  -  |  (4.9242a  -  24.2485)  (6-35) 

When  the  value  of  the  parameter  6  from  equation  (6-35)  is 
substituted  in  the  general  form  of  the  characteristic 
equation  of  the  compensated  singular  system,  i.e.,  into 
equation  (6-17)  yields: 

1+a  Fg(s)+  (4.9242a  -24.2485)  Fg2(s)  =  0 
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or  after  substituting  F  [s)  from  equation  C6-22)  and 
making  the  necessary  manipulations  yields: 


s6+12s5+46s4+(60+a)s3+(25+6a)s2+5as+ 

(4.9242a  -  24.2485)  =  0  (6-36) 

i.e.,  the  characteristic  equation  of  the  compensated 
singular  system  is  independent  of  the  value  of  k  when  the 
same  singular  point  is  considered.   Therefore,  for  the 
same  singular  point  and  for  the  same  value  of  the  parameter 
a  the  compensated  singular  system  will  have  the  same 
characteristic  roots  independently  of  the  value  of  k.   As 
a  consequence  of  these  remarks  not  only  for  k  =  20  but 
for  any  value  of  k  >0  the  Root  locus  of  the  compensated 
singular  system  will  be  identical  with  the  Root  locus 
shown  in  figure  (6-9)  when  the  singular  point  s  =  -0.3884  + 
j  0.8899  is  considered. 

4.    Summary  and  Remarks  on  the  Stability  of  the 
Compensated  Singular  System 

When  the  plant  to  be  compensated  by  a  singular 
cascade  compensator  has  no  zeros  in  its  transfer  function 
then  the  parameter  6  and  the  Plant's  gain  k  will  appear 
only  in  the  zero  power  of  s  term  of  the  characteristic 
equation  of  the  compensated  singular  system  as  can  be 
seen  from  equation  (6-19).   In  addition,  only  the  product 
of  3  and  k  will  appear  in  this  term,  plus  a  constant  when 
the  plant  is  type  zero.   When  a  specific  singular  operating 
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point  of  the  system  is  considered  (Ss>uns)>  then  the  para- 
meters a  and  3  will  be  related  by  the  equation  of  the 
singular  line  which  corresponds  to  the  above  mentioned 
singular  point.   According  to  equation  [2-5)  and  (2-16), 
the  equation  of  the  singular  line  will  be  in  the  form: 

3  =  |  (A  +  Ba)  (6-37) 

where  A  and  B  are  constants.   From  this  equation  it  will 
be  observed  on  the  parameter  plane  that  in  the  case  under 
consideration  when  the  value  of  k  varies  it  will  change 
the  slope  of  the  singular  line  which  passes  from  the 
selected  singular  point.   From  this  equation  also  it  can  be 
seen  that  the  product  3k  can  be  expressed  as  a  function  of 
the  parameter  a  alone.   By  substituting  this  product  in  the 
characteristic  equation  (6-19)  of  the  compensated  singular 
system  the  whole  equation  will  be  expressed  as  a  function 
of  the  parameter  a.   Thus  the  roots  of  the  characteristic 
equation  will  depend  on  the  value  of  a  only. 

When  the  plant  to  be  compensated  by  a  singular 
cascade  compensator  has  one  or  more  zeros  in  its  transfer 
function  then  the  parameter  Band  the  plant's  gain  k  will 
appear  in  three  or  more  power  of  s  terms  of  the  character- 
istic equation  of  the  compensated  singular  system  as  can 
be  seen  from  equation  (6-19).   Because  they  will  always 
appear  together  in  a  product  form  the  coefficient  C.  of 
equation  (2-5)  can  be  written  in  the  form: 

Cl  =  kE  (-1)J  C'k%  j  Uj-iCO        (6-38) 
j=o 
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where  C ',  is  the  constant  coefficient  of  the  3k  product, 
k 

in  the  j power  of  s  term  of  equation  C6-19).   When  a  specific 

operating  singular  point  of  the  compensated  singular  system 

is  considered,  i.e.,  (  r     .m      ),  then  the  parameters  a  and 

hs      ns 

3  will  be  related  by  the  equation  of  the  corresponding 
singular  line  of  the  above  mentioned  singular  point.   Then 
according  to  equation  (2-5)  and  (2-16),  the  equation  of  the 
singular  line  will  be: 

3  =  |  (A1  +  B'oi)  (6-39) 

where  A'  and  B'  are  constants.   By  continuing  with  the  same 
reasoning  as  in  the  previous  case,  the  same  remarks  were 
found  to  be  applied  here. 

Therefore  it  can  be  stated,  that  in  general  the 
stability  of  a  compensated  singular  system  when  a  specific 
singular  point  is  considered,  is  independent  of  the  value 
of  the  plant  gain  k  and  depends  on  the  value  of  the  para- 
meter a  alone  or  in  other  words  the  stability  of  the  compen- 
sated singular  system  can  be  interpreted  in  terms  of  the 
parameter  a  alone.   The  stable  region  also  (if  any)  on  the 
parameter  plane  will  be  bounded  between  constant  a  lines. 
Thus  when  the  compensated  singular  system  is  stable  for 
some  range  of  the  parameter  a  values,  it  will  also  be 
stable  for  the  same  range  of  the  parameter  a  values  and 
for  any  value  of  the  plant  gain  k.   Therefore,  when  stability 
analysis  is  performed  by  the  Root  locus  method  where  the 
parameter  a  is  used  as  the  Root  locus  variable  parameter 

then  the  value  of  k  can  be  set  to  any  desired  value  as  will 
be  dictated  from  other  design  factors. 
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E.    HOW  TO  OBTAIN  A  SINGULAR  ROOT  AT  A  DESIRED  POINT 
OF  THE  S- PLANE 

1 .    Introduction 

Compensation  of  a  system  by  using  a  singular 
cascade  compensator  will  be  in  general  useful  if  the 
compensated  system  provides  a  singular  root  where  the 
specifications  of  the  problem  dictate.   Therefore,  when 
designing  by  using  the  singular  line  theory  the  next 
problem  which  seeks  an  answer  is:  How  a  singular  root  can 
be  obtained,  where  the  system's  specifications  dictate? 

The  general  objective  in  all  cases  will  be  to 
make  use  of  the  advantage  that,  by  varying  the  parameter 
values  the  singular  root  remains  fixed  and  as  other  roots 
are  varying  the  system's  performance  may  be  moved  toward 
the  system's  specifications. 


2.    Selection  of  the  Singular  Point 

Assume  that  a  plant  with  transfer  function: 

Cs  +  Zl)Cs  +  z2)     ...     (s  +  z    )  f6-4fn 

G(s)    =    kF    (s)    =   k— « 2—  ^    4UJ 

5  s    C.s+p1)(s+p2)    ...    (s+pn) 

is  going  to  be  compensated  by  a  "singular  cascade  compensa- 
tor" whose  transfer  function  according  to  equation  (6-18)  is 
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(s  +  z  )  (s  +  z  )......  (s+z  ) 

sc        sN(s+p  )(s+p2)  .  ..(s+pn)+  a  (s  +  z^Cs  +  z^  ...(s  +  zm) 

(6-41) 
then  the  characteristic  equation  of  the  compensated 
singular  system  will  be: 

1  +  aF  (s)  +  gk  F  2(s)  =  0  (6-42) 

According  to  equation  (3-10)  the  singular  points 
of  the  system  will  be  the  points  of  the  s-plane  which 
satisfy  the  following  equation: 

/Fg(s)  =  kxH  (6-43) 

where  k,  is  a  real  integer  or  zero.   Therefore,  the  location 
of  the  singular  points  on  the  s-plane  is  determined  only 
by  the  zeros  and  the  poles  of  the  plant,  i.e.,  the  singular 
points  are  a  function  of  F  (s)  only,  and  in  consequence 
they  can  be  affected  only  by  changing  F  (s) . 

In  the  case  where  among  the  singular  points  there 
is  such  a  point  which  satisfies  problem  specifications  as  the 
operating  point  of  the  compensated  singular  system,  then 
there  is  no  need  to  change  F  (s) .   Considering  now  the 
case  where  such  a  singular  point  is  not  available  then  F  (s) 
has  to  be  changed  in  order  to  obtain  a  singular  point  where 
it  is  dictated  by  the  problem's  specifications.   The  fol- 
lowing examples  were  selected  for  illustration  of  the  design 
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procedure  in  the  above  case, 
a.   Example  (6-3) 

Assume  the  plant  to  be  compensated  has  the 
transfer  function: 


G(s)  =  kF  (s)  =  i (6-44) 

&  s(s  +  l)(s  +  5) 

and  that  problem  specifications,  i.e.,  damping  ratio, 
natural  undamped  frequency,  settling  time  and  peak 
overshoot  require  a  dominant  singular  root  at  the 
point  (?s  ■  0.4,  tons  =  0.97096)  or  ss  =  -  0.3884  ±  J0.8899 

In  figure  (6-11)  the  poles  of  Fg  (s)  have  been 
plotted  on  the  s-plane  as  reference  points,  in  order  to 
measure  the  angles  involved.   The  desired  singular  point 
is  also  located  on  the  s-plane.   By  the  use  of  a  spirule 
it  was  found  that: 

/Fg(s)  =  -  {*1  +  <t>2  +  <j>3}  =  -  180° 
therefore  the  point  S   is  actually  a  singular  point  since 
for  this  point  equation  (6-43)  is  satisfied  fork-,  =  -1. 

This  point  also  was  found  to  be  a  singular 
point  by  the  use  of  the  "Singular  point  program"  in 
example  (6-1)  . 

It  must  be  noticed  that  another  way  to  check 
if  the  point  (c   =  0.4,w   =0.97096)  is  singular  for  the  . 
compensated  singular  system,  is  to  draw  the  Root  locus  of 
the  uncompensated  system.   Then  if  the  desired  operating 
point  lies  on  the  above  Root  locus  (according  to  the 
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analysis  which  was  developed  in  paragraph  2a(3),  Section 
III),  it  is  singular.   Because  this  method  is  more 
laborious  its  use  in  general  is  not  justified,  except  in  the 
case  where  the  desired  operating  point  is  not  singular  and 
some  trade  offs  have  to  be  done  in  order  to  decide  if  Fg(s) 
must  be  changed  [by  the  appropriate  initial  compensator) 
or  another  operating  point  which  lies  on  the  Root  locus  of 
the  uncompensated  system  can  be  selected  as  the  operating 
point  of  the  compensated  singular  system. 

b.   Example  (6-4) 

Assume  now  that  the  plant  to  be  compensated 
has  transfer  function: 

rr     S  1  -c     r     ->     i         s  +  5  (6-45) 

G(s)  =  kFp(s)  =  k  -* 

&  sz(s+10) (s+20) 

and  that  problem  specifications  dictate  for  an  operating 
singular  point  (C   =0.6,w   =1),  or  S  =  -0.6±j  0.8.   As  in 
the  previous  example  the  poles  and  zeros  of  Fg(s)  were  plot- 
ted on  the  s  plane  as  shown  in  figure  (6-12).   Then  by  the 
use  of  a  spirule  the  argument  of  Fg(s)  was  evaluated  and  it 
was  found  that: 

/Fg(s)  =  81  -(2^+f  +<!>3)  =  22.5°  -  (2xl27°  +  10.5°  +  5°)  =-247° 

(6-46) 
Therefore  the  desired  point  is  not  singular  since  it  does 

not  satisfy  equation  (6-43).   This  implies  that  Fg(s)  must 

be  changed  somehow  in  order  for  the  desired  point  to  become 
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singular.   From  a  theoretical  point  of  view  it  can  be 
achieved  as  follows: 

(1)   Introducing  an  Additional  Zero.   Suppose 
that  an  additional  zero  is  introduced  at  the  point  s=-1.82 
Then  the  angle  contribution  of  this  zero  will  be  +67°  as 
it  can  be  seen  in  figure  (6-13a).   Then,  by  defining  the 
transfer  function  of  the  new  obtained  system  after  the 
above  zero  was  introduced  G  (s) ,  yields: 


Gz(s)  =  kFz(s)  =  k{(s+1.82)  Fg(s)} 


and 


L 


Fz(s)  =  -  180' 


which  implies  that  the  desired  operating  point  (c   =  0.6, 


=  1)  is  singular  for  the  system  whose  transfer  function 


ns 

G  (s)  -  k  (f5)  Cs  +  1.82) 
Z        sZ(s+10)  (s+20) 

(2)   Introducing  an  Additional  Pole  .   Suppose 
that  instead  of  a  zero  an  additional  pole  is  introduced 
at  the  point  s=-0.58,  then  the  angle  contribution  of  this 
pole  will  be  -113°  as  it  is  shown  in  figure  (6-13b) .   Then 
by  defining  the  transfer  function  of  the  new  system  after 
the  pole  was  introduced,  G  (s)  yields: 


and 


G  (s)  =  k  FpCs)  =  k  {iT^Fg(s)} 


A  Is)  ". 


360' 
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which  implies  that  the  desired  operating  point  (?  =0.6, 

u  =  1)  is  singular  for  the  system  with  transfer  function: 

s  +  5 


G  (s)  =  k  -n 

P        s  (s+0.58) (s+10) (s+20) 

(3)    Introducing  an  Additional  Pole  and  Zero  . 

Suppose  now  that  an  additional  pole  and  a  zero  are 
introduced  at  the  transfer  function  of  the  initially 
given  plant.   This  pair  of  points  can  be  chosen  in  such 
a  way  that  their  net  angle  contribution  will  make 
the  desired  singular  point  (Cs  =  0.60,wns=  1)  actually 
singular.   For  example  if  the  pole  is  chosen  at  the  point 
s  =  -  16  and  the  zero  at  the  point  s  =  -  1 . 66  as  it  is  shown 
in  figure  (6-13c),  then  by  defining  the  transfer  function  of 
the  new  system  after  this  pole  and  zero  have  been  introduced, 
Gzp(s)  yields: 


and 


Gz:p(s)  =  k  F7  fs)  =  k  {  S  +  1-66  F  (s)} 
Zp         zp  s  +  16    g 

/Fg   (s)  =/Vs)   +  73°  "  6°  =  "  180° 

which  implies  that  the  desired  operating  point  is  singular 

for  the  system  with  transfer  function: 

(s+5)  (s+1. 66) 

G   (s)  =  k 

zp 

s2(s+10)  (s+16)  (s+20 
3.   Summary  and  Remarks 

Considering  in  general  an  uncompensated  plant  with 
open  loop  transfer  function  G(s)  =  k  Fg(s),  where  Fg(s) 
is  a  ratio  of  any  two  polynomials  of  s,  then  the  argument 
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of  F  (s)  evaluated  at  any  point  of  the  s-plane  will  always 
lie  between  the  limits: 

k1n  </f8Cs)  <  (k^i)  n 

where  k,  can  be  any  real  integer  or  zero.   Therefore,  the 
maximum  angle  correction  which  F  (s)  may  require  in  order 
to  satisfy  equation  (6-43)  for  a  desired  point  of  the 
s-plane  will  be  less  than   180°  ,  which  implies  that  a 
single  zero  or  a  single  pole  or  a  zero  and  a  pole  or  finally 
any  number  of  zeros  and  poles  can  be  properly  selected  and 
introduced  to  F„(s)  in  order  to  make  the  argument  of  F„(s) 
satisfying  the  angle  requirement  for  singular  conditions, 
for  a  specific  singular  point,  which  is  desired  to  be  the 
operating  point  of  the  compensated  singular  system. 

These  poles'  and  zeros  in  general  may  be  real  or 
complex.   In  consequence  it  implies  that  when  such  a  change 
of  Fg(s)  is  required  an  initial  cascade  compensator  must  be 
introduced,  before  a  "singular  cascade  compensator"  is 
introduced. 

Since  a  compensator  with  only  zeros  actually  can 
not  be  built,  the  case  of  introducing  one  or  more  zeros 
alone  has  to  be  ignored  as  unrealistic.   Therefore,  among 
the  rest  of  the  possibilities  available  an  optimal  selection 
of  the  initial  cascade  compensator  has  to  be  done.   This 
selection  can  be  based  on  noise  considerations,  problem 
specifications  and  the  "dominance  of  the  singular  root" 
concept  which  is  described  in  the  next  part  of  this  section, 
i.e.,  except  the  fact  that  the  Root  locus  of  the  uncompensated 
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system  will  be  forced  to  pass  through  the  point  where  the 
desired  singular  root  is  located,  the  shape  also  of  the 
resulting  Root  locus  must  be  such  that  it  promises  some 
range  of  dominancy  of  the  singular  root.   (See  Example 
(6.6)). 

It  must  also  be  noticed  that  as  it  will  be  shown 
in  the  following  analysis  when  an  initial  cascade  compensa- 
tor is  required,  then  this  requirement  will  actually  affect 
only  the  transfer  function  of  the  singular  cascade  compensa- 
tor and  neither  the  complexity  of  the  rest  of  the  design 
problem  nor  the  compensation  structure  (See  Example  (6-6) 
and  figure  (6-28c)).  Under  these  considerations  an  initial 
cascade  compensator  with  only  zero's  may  also  be  considered. 

F.    DOMINANCE  OF  THE  SINGULAR  ROOT 
1 .   Dominant  Root  Concept 

The  concept  of  the  dominant  root  has  wide  applica- 
tions in  analysis  and  design  of  linear  systems.   Considering 
a  system  with  a  characteristic  equation  of  order  n,  then 
there  will  be  n  characteristic  roots  for  this  system.   If 
one  of  these  roots  dominates  the  transient  response,  then 
this  root  will  be  called  dominant.   When  the  dominant  root 


of  a  system  is  a  complex  root,  s  =-500  ±joi   /l-e;  ,  then  the 
overall  response  of  the  system  for  a  given  input  may  be 
closely  approximated  from  the  corresponding  standard  second 


order  system  graphs  (i.e.,  the  second  order  system  whose 

the  ch 
zeros) 


2  2 

the  characteristic  equation  is  s  +2?w  s+w  =0  and  has  no 

n  n   n 
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2 .  Dominance  of  the  Singular  Root 

When  a  singular  root  has  been  determined  from 
the  problem  specifications,  then  this  singular  root  may 
be  available  from  the  given  plant  or  not.   In  the  last 
case  an  initial  cascade  compensator  will  be  required,  in 
order  for  the  desired  singular  root  to  become  available. 
The  design  process  of  this  initial  cascade  compensator 
was  previously  illustrated. 

In  either  case  it  will  also  be  desired  for  the 
above  considered  singular  root  to  be  dominant  for  one  or 
more  (depends  from  the  specific  application)  parameter  pair 
values.   The  design  process  in  all  possible  cases  toward 
the  dominancy  of  the  singular  root  will  be  illustrated  by 
the  following  examples: 

3 .  Example  (6-5) 

In  this  example  the  Plant  of  example  (6-2)  was 
again  considered  with  the  Root  locus  gain  k  =  20,  i.e.: 

G(s)  -  k  Fg(s)  -  20  sCstl}(st5) (6-47) 

Since  the  value  of  the  plant's  gain  does  not 
affect  the  location  of  the  singular  roots  on  the  s-plane, 
as  was  shown  in  example  (6-2),  the  plant  under  considerati 
when  it  is  compensated  by  a  "singular  cascade  compensator" 
will  have  the  same  singular  roots  with  these  which  are 
listed  in  Appendix  C  for  the  example  (6-1)  and  (6-2).   S 

of  these  roots  were  considered  and  their  dominancy  was 

investigated. 
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In  figure  (6-14a)  and  (6-14b)  the  compensated 
singular  system  is  shown,  where  in  figure  C6-14a)  the 
singular  compensator  is  considered  to  be  in  integrated 
form,  and  in  figure  (6-14b)  the  same  singular  compensator 
is  considered  to  consist  of  distinct  components. 

The  singular  cascade  compensator  which  thus  has 
been  introduced  according  to  equation  (6-18)  has  transfer 
function: 
r  fsi  -      3  C6"48) 

scLSJ    s(s+l)  (s  +  5)+a 

The  characteristic  equation  of  the  compensated 
singular  system  was  found  to  be: 

s6+12s5+46s4+(60+a)s3+(25+6a)s2+5as  +206  =  0       (6-49) 

a.   Singular  Root  (?   =  0.40,wns=  0.97096) 

When  the  singular  root  (?s=  0.40,cons=  0.97096) 

was  considered,  then  the  problem  which  was  investigated  was: 

Are  there  any  pairs  of  parameter  a  and  g  values  for  which  this 

singular  root  becomes  dominant? 

From  equations  (2-5)  and  (2-16)  the  equation 

of  the  corresponding  singular  line  was  found  to  be: 

3  *  0.24631a  -  1.21242  6-50) 

Introducing  equation  (6-50)  into  equation 
(6-49)  the  characteristic  equation  of  the  compensated 
singular  system  in  terms  of  the  parameter  a   becomes: 
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s6+12s  +  46s  +(60+a)s  +(25+6a)s  +  5as 

(6-51) 
+(4. 9242a  -  24.2485)  =  0 


Treating  then  the  parameter  a  as  the  Root  locus 
variable  parameter,  the  Root  locus  of  the  compensated 
singular  system  for  the  singular  point  under  consideration 
was  drawn  in  the  computer.   It  is  shown  in  figure  (6-9). 
The  compensated  singular  system  has  three  fixed  roots,  i.e., 
the  singular  root  s  =-0.3884  ±j  0.8899  and  a  real  root 
s=-5.2230,  which  are  shown  with  dark  circle  marks  at  the 
corresponding  locations  of  figure  (6-9).   The  compensated 
singular  system  also  has  three  other  roots  whose  location 
varies  on  the  corresponding  Root  locus  segments  according 
to  the  value  of  the  parameter  a . 

Examining  the  roots  of  the  characteristic  equation 
C6-51)  in  relation  with  the  Root  locus  of  figure  (6-9) 
it  can  be  seen  that  the  value  of  the  plant  gain  does  not 
affect  the  location  of  the  roots  of  the  characteristic 
equation  of  the  compensated  singular  system  and  that  the 
location  of  the  three  varying  roots  is  affected  only  by 
the  value  of  the  parameter  a .   The  dominancy  also  of  the 
singular  root  can  be  affected  only  by  the  other  two 
characteristic  roots  which  lie  on  the  same  Root  locus  seg- 
ments with  the  singular  root.   These  roots  are  either  real 
or  complex  according  to  the  value  of  the  parameter  a  . 
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It  was  found  that  in  the  stable  region  the  value 
of  the  parameter  a  is  approximately  betwen  the  values  of 
five  and  fourty.   When  a=  9.85  then  the  two  varying  roots 
in  the  vicinity  of  the  singular  root  become  identical  with 
the  singular  root,  i.e.,  for  the  value  of  the  parameter 
a  =  9.85  the  singular  root  is  repeated.   It  was  also  found 
that  a=9.85  not  only  the  singular  root  but  and  the  real 
fixed  root  ats=-5.2230  was  also  repeated. 

In  an  attempt  to  find  the  values  of  the  parameter 
a  for  which  the  singular  root  is  dominant  (or  in  other 
words  the  system's  performance  can  be  well  predicted  from 
the  corresponding  to  the  singular  root  second  order  model) 
several  runs  with  different  parameter  pair  values  each, 
were  performed  in  the  computer  using  a  DSL  (Digital 
Simulation  Language)  program.   In  this  program  the  compensa- 
ted singular  system  shown  in  figure  (6-14b)  was  implemented. 

When  the  second  order  model  corresponding  to  the 

singular  root  (£   =0.40,w   =0.97096  was  considered  with 
to  s        ns 

transfer  function: 

giSl  =  ^ =  0-9428  (6.52) 

Rts)    s2  +  2cswnss  +  ojns2    s2  +  0.77685s  +  0.9428 

its  transient  response  characteristic  values  were  calculated. 
These  values  together  with  those  which  were  obtained  from 
the  computer  for  the  compensated  singular  system  are 
listed  in  Table  I  of  Appendix  F. 
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By  examining  this  table  the  following  remarks  can  be  stated: 

(1)  There  is  no  pair  of  parameter  values  for 
which  all  the  transient  response  characteristics  of  the 
compensated  singular  system  can  be  well  predicted  from 

the  values  corresponding  to  the  singular  root  second  order 
model . 

(2)  For  all  the  parameter  pairs  both  the 
singular  root  and  the  other  root  in  the  vicinity  of  the 
singular  root  contribute  to  the  transient  response  of  the 
compensated  singular  system. 

C3)   Considering  that  in  practice  the  only 
transient  response  characteristic  values  for  which  the 
designer  is  in  general  interested,  are  the  settling  time 
and  the  percent  overshoot,  then  from  this  point  of  view, 
it  can  be  said  that  the  singular  root  is  dominant  for 
a  =7.8  and  6 =  0.708.   Assuming  also  that  there  is  some 
tolerance  in  the  system  performance  then  these  two  parameter 
pair  values  can  vary  within  some  range  according  to  the 
tolerance  values. 

The  conclusion  from  this  example  was 
that  when  dominant  mode  design  is  used  for  compensation  of 
the  given  plant  then  the  singular  line  theory  provides  in 
general  an  acceptable  solution.   In  figure  (6-15)  the  time 
response  of  the  compensated  singular  system  is  shown  for 
a  unit  step  input  and  for  the  two  parameters  having  the 
values  of  a=  7.8  and  g=  0.71.   In  figure  (6-16)  the  time 
response  of  the  compensated  singular  system  is  shown  for 
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the  parameter  pair  values  o.~    9.85  and  8  "  1.21  for  which  the 
singular  root  and  the  fixed  real  root  are  repeated.   This 
fact  as  can  be  seen  from  this  figure  causes  a  very  high 
overshoot  value.   In  figure  [6-17)  the  time  response  of 
the  compensated  singular  system  is  shown  for  such  a  pair 
of  values,  i.e., a  =20  and  8  =  3.71,  for  which  the  contribution 
of  the  singular  root  and  the  other  root  in  the  vicinity  of 
the  singular  root  can  be  recognized.   In  figure  (6-18)  the 
time  response  of  the  compensated  singular  system  is  shown 
for  the  case  where  the  characteristic  root  in  the  vicinity 
of  the  singular  root  lies  on  the  real  axis.   In  this  case 
because  this  root  is  closer  to  the  origin  and  to  the  imaginary 
axis  than  the  singular  root,  it  is  dominant  and  the  transient 
response  of  the  compensated  singular  system  has  an  exponential 
form. 

b.   Singular  Root  (£s  =  0  .  30,tons  =1 .  15268) 

The  second  singular  root  from  the  table  of 
Appendix  C  which  was  examined  for  dominancy  was  the  singular 
root  (?s  =  0.30,wns=  1.15268)  or  s  =-0.34580  ±  j  1.09959. 

From  equations  (2-5)  and  (2-16)  the  equation 
of  the  corresponding  singular  line  was  found  to  be: 

3  =  0.35266a  -  2.48732  (6-53) 

Introducing  equation  (6-53)  into  equation 
(6-49)  the  characteristic  equation  of  the  compensated 
singular  system  in  terms  of  the  parameter  abecomes: 

s  +12s5+46s4(60+a)s3+(25+6a)s2+5as+7. 0532a  -  49.7464  =  0 

(6-54) 
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Treating  then  the  parametera  as  the  Root  locus 
variable  parameter,  the  Root  locus  of  the  compensated 
singular  system  for  the  singular  root  under  consideration 
was  drawn  in  the  computer.   It  is  shown  in  figure  (6-19). 
This  Root  locus  as  can  be  seen  from  this  figure  consists 
of  three  points  of  the  s-plane  where  the  three  fixed  roots 
of  the  compensated  singular  system  lie  and  the  classical 
Root  locus  segments.   The  above  mentioned  three  fixed  roots 
are  shown  with  dark  circle  marks  at  the  corresponding 
locations ,  i.e., 

s,  =  -  5.3084 

s2    =    -    0.34580+    jl. 09959 

s3  =  -  0.34580  -jl. 09959 

In  an  attempt  to  find  the  values  of  the  parameter 

a   for  which  the  singular  root  is  dominant  several  runs  with 

different  parameter  pair  values  each,  were  performed  in  the 

computer  using  the  DSL  simulation  program.   The  transient 

response  characteristic  values  which  were  obtained  from 

the  computer,  together  with  these  corresponding  to  the 

singular  root  (t   =0.30, to   =  1.15268)  second  order  model 

s        ns 

are  listed  in  Table  II  of  Appendix  F. 

By  examining  this  table,  the  same  general  remarks 
and  conclusions  were  derived  in  this  case  as  in  the  previous 
one.   The  basic  remark  again  was  that  there  is  a  pair  of 
parameter  values,  i.e.,  a=  10.6,  $  =  1.25,  for  which  the 
singular  root  is  dominant  from  a  practical  point  of  view. 
In  figure  (6-20)  the  time  response  of  the  compensated 
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singular  system  is  shown  for  this  pair  of  parameter  values, 
In  figure  also  (6-21)  through  (6-23),  the  time  response  of 
the  compensated  singular  system  is  shown  for  several 
characteristic  pairs  of  parameter  values. 

c  Singular  Root  (5   =  0.75 ,^ns  =  0.60667) 

The  third  singular  root  which  was  examined 

for  dominancy  was  (£   =  0.75, to   a  0.60667)  or 

s  =  -0.45500  ±  j  0.40127. 

From  equations  (2-5)  and  (2-16)  the  equation 

of  the  corresponding  singular  line  was  found  to  be: 

3  -  .0. 09367a  -  0.17547  (6-55) 

Introducing  equation  (6-55)  into  equation 

(6-49)  the  characteristic  equation  of  the  compensated 

singular  system  in  terms  of  the  parameter  a  becomes : 


s6  +  12s5  +  46s4+(60  +  a)s3+(25  +  6a)s2  +  5as  +  1.8734a-3.5094  =  0  (6-56) 


Treating  then  the  parameter  aas  the  Root  locus 
variable  parameter,  the  Root  locus  of  the  compensated 
singular  system  for  the  singular  root  under  consideration 
was  drawn  in  the  computer.   It  is  shown  in  figure  (6-24)  . 
This  Root  locus  (as  can  be  seen  from  figure  (6-24)  also 
consists  of  three  points  of  the  s-plane  where  the  three 
fixed  roots  of  the  compensated  singular  system  lie  and  the 
classical  Root  locus  segments.   The  three  points  where  the 
three  fixed  roots  are  located  are  shown  with  dark  circle 
marks  at  the  corresponding  locations,  i.e.: 
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s,=  -  5.0899 

s2=  -  0.45500+  JO. 40127 

s3=  -  0.45500-  JO. 40127 

By  the  same  procedure  used  in  the  previous  cases 
the  same  general  remarks  and  conclusions  were  derived  in 
this  case  also.   The  basic  remark  again  was  that  there  is 
a  pair  of  parameter  values  a  and  3  for  which  the  singular 
root  is  dominant  under  a  practical  point  of  view.   This  pair 
as  can  be  seen  from  the  Table  III  of  Appendix  F  is:  ot  =  3.6, 
3  =  0.162.   In  this  table  the  transient  response  character- 
istic values  which  were  obtained  from  the  computer  together 
with  these  corresponding  to  the  singular  root  (c   =  0.75, 

w   =  0.60667)  second  order  model  are  listed.   In  figure 
ns 

(6-25)  the  time  response  of  the  compensated  singular  sys- 
tem is  shown  for  the  previous  mentioned  parameter  pair 
values . 

d.   Remarks 

It  must  be  noted  here  that  the  proceeding 
design  was  available  for  the  particular  system.   It  may 
not  be  generally  available  or  may  not  provide  an  accept- 
able solution  to  some  problems. 

4 .   Root  Loci  Correlation 

In  the  previous  example  if  the  desired  singular 
root  of  the  compensated  singular  system  could  not  be  made 
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dominant  by  a  pair  of  parameter  values,  then  an  initial 
cascade  compensator  would  be  required,  as  has  been  already 
stated  in  paragraph  two. 

Since  the  initial  cascade  compensator  is  applied 
to  the  uncompensated  system,  i.e.,  the  system  before  the 
singular  cascade  compensator  has  been  introduced,  in  order 
to  affect  the  compensated  singular  system,  i.e.,  the 
system  after  the  singular  cascade  compensator  has  been 
introduced,  a  prior  knowledge  of  this  effect  is  necessary 
in  order  to  avoid  an  undetermined  number  of  trials. 

This  prior  required  knowledge  or  the  criterion 
for  the  selection  of  the  initial  cascade  compensator  can 
be  based  on  the  correlation  between  the  Root  locus  of  the 
uncompensated  system  and  the  Root  locus  of  the  compensated 
singular  system,  as  can  be  seen  from  the  following  analysis 

When  an  uncompensated  linear  control  system  with 
unity  feedback  is  considered,  its  characteristic  equation 
will  in  general  have  the  form: 

(s  +  z  )  (s  +  2  -)  ...  (s  +  z  ).  .  .  . 

k  — jp-J J 2 =  k  F  (s)  =  -  1       (6-57) 

sN(s+Pl)(s+p2)  ...  (s+pn)       g 

where  k  Fg(s)  is  the  open  loop  transfer  function  and  k 
is  the  Root  locus  gain.   The  Root  locus  of  this  uncompen- 
sated system  which  is  drawn  with  the  Root  locus  gain  k  as 
the  Root  locus  variable  parameter  consists  of  all  the 
points  of  the  s-plane  which  satisfy  the  angle  requirement: 
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A  (s)  =  C2n±i)n 


Therefore  any  point  SRT  of  the  Root  locus  of  the  uncom- 
pensated system  as  it  can  be  seen  from  equation  (6-57) 
and  the  previous  angle  relation  satisfies  the  following 
two  relations: 


k  Fg(SRL)  =  -  1 


and 


z 


Fg(SRL)  =  C2n±l)n 


or 


FgCSRL)  <   0  (6-58) 


and       /Fg(SRlP  =  ^2n±l)  n  (6-59) 


where  n  is  a  real  integer. 

When  the  system  which  was  considered  is  com- 
pensated by  a  singular  cascade  compensator  the  charac- 
teristic equation  of  the  compensated  singular  system 
according  to  equation  (6-17)  will  be: 


1+a  Fg(s)  +3  k  Fg2(s)  =  0  (6-60) 


Then  according  to  equation  (3-12)  a  point  of 
s-plane  is  a  singular  point  for  the  compensated  singular 
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system  if  and.  only  if  it  does  not  lie  on  the  real  axis 
and  satisfies  equation: 


L 


F  fs)  =  nk  (6-61) 


where  n  is  a  real  integer  or  zero.   Therefore  all  the 
points  of  the  Root  locus  of  the  uncompensated  system,  which 
do  not  lie  on  the  real  axis  (this  restriction  is  implied 
from  the  definition  of  the  singular  point  in  Section  II, 
A)  are  singular  points  because  they  satisfy  equation  (6-59) 
and  in  consequence  equation  (6-61) . 

Suppose  that  one  of  these  points,  i.e.,  the 
point  (c  ,w   )  is  selected  as  the  operating  point  for  the 
compensated  singular  system.   For  this  point  equation 
(6-58)  becomes: 

V^s'V,)  <C  °  (6-62) 

The  equation  of  the  corresponding  singular  line 
(i.e.,  the  equation  which  defines  the  relation  between  the 
two  parameters  a  and  3  for  which  the  selected  singular 
point  will  always  be  a  root  of  the  characteristic  equation) 
can  be  found  either  by  using  equations  (2-5)  and  (2-16) 
or  by  direct  substitution  of  the  selected  singular  point 
in  equation  (6-60),  which  yields: 

1+a   FgC?.s,wns)    +6   k   Fg2'CCs,unsD    =    0  (6-63) 
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Solving  this  equation  for  3k  and  substituting 
it  in  equation  (6-60)  yields: 


1+a  Fg  CC    ,w      ).  ? 

l+«      Fg(s) 2 5-J1S F^z(s)    m    Q  (6_64) 

Fa2C?    ,<o      ) 
g    ^   s'   ns; 


Equation  (6-64)  is  the  characteristic  equation 
of  the  compensated  singular  system  as  a  function  of  the 
parameter  a,  when  the  singular  point  {£    ,.oi   )  is  considered 
to  be  the  operating  point  of  the  compensated  singular  sys- 
tem.  This  equation  can  be  written  in  the  form: 


0  {w%s)Fg(s)}  {F8cgS'%s>  ;  FgCs)}  = .  i 


(6-65) 


or 

(6-66) 

FK(C  ,a>   )(s  +  z  )(s  +  z  )...(s+z  ) 
a   5   s  ns i 2 m n 

Fg  (C  .w   )sN(s+p  )(s+p  )  .  .  .  (s+p  )  +  (s  +  z  )(s  +  z  )...(s+z  ) 
&   s  ns       12        n      i     2        n 


Considering  the  Root  locus  of  the  uncompensated 
system,  the  Root  locus  of  the  compensated  singular  system 
and  the  corresponding  equations,  the  following  remarks 
and  conclusions  can  be  stated. 
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a.  Number  of  Fixed  Roots 

From  equation  (6-65)  ,  the  common  factor  in 
the  numerator  and  the  denominator  on  the  left  side  of 
equation,  implies  that  the  compensated  singular  system 
has  a  number  of  fixed  roots  equal  to  the  order  of  the 
characteristic  equation  of  the  uncompensated  system,  e.g., 
when  the  order  of  the  characteristic  equation  of  the 
uncompensated  system  is  three,  there  will  be  three  fixed 
roots,  two  of  which  have  to  be  the  complex  pair  of  the 
singular  roots. 

b.  Root  Loci  Correlation 

The    angle   relation  which   the   points    of   the 
Root   locus   of   the   compensated   singular   system   satisfy   is: 


FgU    ,w      )F    (s) 

s  =    (2n±l)n  (6-67) 


WvP+Vs) 


Then  according  to  equations  (6-58)  and  (6-62)  all  the 
points  of  the  Root  locus  of  the  uncompensated  system  satis 
fy  equation  (6-67)  .   Therefore  all  the  points  of  the  Root 
locus  of  the  uncompensated  system  will  be  also  points  of 
the  Root  locus  of  the  compensated  singular  system. 

Second,  from  equation  (6-57)  and  (6-66)  it 
is  observed  that  they  have,  at  the  expression  on  the  left 
side,  identical  zeros.   They  also  have  the  same  number  of 
poles . 
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The  conclusion  from  the  above  remarks  is 
that  the  Root  locus  of  the  compensated  singular  system  will 
be  identical  with  the  Root  locus  of  the  uncompensated  sys- 
tem with  the  only  difference  that  the  Root  locus  segments 
of  the  compensated  singular  system  may  be  further  extended 
from  the  pole  side  end  points,  where  the  corresponding 
Root  locus  segments  of  the  uncompensated  system  terminate, 
until  they  reach  the  poles  of  the  compensated  singular 
system. 

Examining  the  figures  (6-7),  (6-9),  (6-21) 
and  (6-25)  they  confirm  all  the  above  derived  remarks  and 
conclusions.   Therefore,  if  an  initial  cascade  compensator 
is  introduced  to  the  uncompensated  system  this  will  affect 
the  Root  locus  of  the  uncompensated  system  and  the  Root 
locus  of  the  compensated  singular  system  in  the  same  way, 
i.e.,  the  prior  knowledge  of  the  effect  of  the  initial 
compensator  on  the  compensated  singular  system  can  be 
obtained  by  studying  the  effect  of  the  initial  compensator 
on  the  uncompensated  system. 

5.   Range  of  Dominancy  of  the  Singular  Root 

When  a  singular  compensated  system  is  considered 
and  an  operating  singular  point,  then  the  parameters  a  and 
3  will  be  related  by  the  equation  of  the  corresponding 
singular  line.   In  general  there  will  be  two  ranges  of 
a  and  $  parameter  values  of  primary  interest.   The  first 
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one  will  be  that  for  which  the  compensated  singular  system 
is  stable  and  the  second  one  will  be  that  within  the  pre- 
vious one,  in  which  the  singular  root  is  dominant.   In 
general  the  wider  these  ranges  the  more  desirable  will  be 
the  solution.   Therefore  when  an  initial  cascade  compensator 
will  be  introduced  to  the  uncompensated  system  in  order  to 
make  the  singular  root  dominant,  it  must  be  selected  in 
such  a  way  that  it  will  also  give  a  wide  range  of  dominancy 
of  the  singular  root. 

It  must  be  realized  that  in  a  self-adaptive 
singular  system  where  the  two  parameters  a  and  3  are 
automatically  adjusted  so  that  the  singular  point  will 
remain  on  the  singular  line  during  the  operation  of  the 
system,  the  range  of  dominancy  of  the  singular  root  will 
be  very  important. 

The  following  example  is  mainly  referred  to  the 
range  of  dominancy  of  the  singular  root. 

6 .    Example  (6-6) 

In  this  example  the  same  plant  which  is  shown 
in  figure  (6-26)  was  again  considered.   Its  transfer 
function  is: 


GCs)  "  s(s  +  l)(s  +  5) 
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For  this  example,  it  was  assumed  that  problem  specifica- 
tions require  an  operating  point  for  the  system,  i.e.,  an 
operating  dominant  singular  root,  with  damping  ratio 
£  =  0.60,  and  natural  frequency  to  =1.2  rad/sec.   A  wide 
range  of  dominancy  was  also  assumed  to  be  required  for  the 
above  singular  root. 

The  Root  locus  of  the  plant  with  the  Root  locus 
gain  used  as  the  Root  locus  variable  parameter  is  shown  in 
figure  (6-29),  where  the  location  of  the  desired  singular 
root  is  also  shown.   By  studying  this  figure  and  taking 
into  account  the  remarks  which  were  stated  in  part  VI  E 
and  VI  F,  the  conclusion  was  that  an  initial  cascade  com- 
pensator is  needed  for  the  following ■ two  reasons:   The 
first  is  that  since  the  desired  singular  root  does  not 
lie  on  the  Root  locus  of  the  plant  this  Root  locus  must  be 
reshaped  in  order  to  pass  through  the  point  where  the 
desired  singular  root  is  located,  and  the  second  reason  is 
that  the  shape  of  the  Root  locus  must  be  changed  in  order 
for  this  singular  root  to  become  dominant  with  a  wide  range 
of  dominancy,  as  it  is  required  from  the  problem  specifi- 
cations . 

After  the  above  remarks  and  taking  into  account 
the  noise  problem  it  was  decided  to  compensate  the  plant 
initially  by  a  cascade  lead  filter  with  transfer  function: 

Ms)  -  S  +  Z 


c^J    ~    s  +  20  (6-68) 
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The  compensated  system  is  shown  in  figure  (6-27) 
The  characteristic  equation  of  this  system  is: 

s4+26s3+125s2+(100+k)  s+2k  =  0  (6-69) 

where  k  =  20.   Equation  (6-69)  can  be  written  in  the  form: 
,.  s  +  2 


s(s  +  l)  (s  +  5)(s  +  20) 


=  -  1  (6-70) 


The  Root  locus  of  this  system  with  k  used  as  the 

Root  locus  variable  parameter  is  shown  in  figure  (6-30). 

At  z  =  0.6  the  corresponding  Root  locus  point  is 

U  =  0.60,  cj   =  1.232586).   The  value  of  w   although  it 
v  n  n       6 

could  be  measured  directly  from  the  Root  locus  diagram  of 

figure  (6-30)  was  found  in  the  indicated  accuracy  by  the 

use  of  the  "singular  point  program".   Therefore  this 

initially  compensated  plant  provides  the  desired  singular 

root  (assuming  that  the  small  difference  between  the  desired 

value  of  a)  ,  i.e.,  w  =  1.2  rad/sec  and  the  actual  value 
n'      '   n 

of  a)   =  1.232586  is  within  the  system's  tolerance).   The 
n.  J 

shape  also  of  the  Root  locus  of  figure  (6-30)  is  such  that 
promises  some  range  of  dominancy  of  the  singular  root. 

Based  on  these  remarks  the  initially  compensated 
plant  was  next  compensated  by  a  singular  cascade  compen- 
sator, the  transfer  function  of  which  according  to  equation 
(6-18)  is: 

G   (s)  =  - g  C5  +  2) <6-71) 

sc      s(s+l)  C.S-+5)  (s  +  20)  +a  (s  +  2) 
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The  compensated  singular  system  is  shown  in 
figure  (6-28a) •   The  characteristic  equation  of  this 
system  was  found  to  be: 

s8+52s7+926s6+(67Q0+a)s5+(2Q825+2  8a)s4+ 

C25000+177a)s3+(10000+350a  +20B)s2+  (6-72) 

(200a+806)s+806  =  0 

When  the  singular  root  CC   =  0.60,  u   =  1.232586) 
to  s       '   ns  J 

was  considered,  the  equation  of  the  corresponding  singular 
line  was  found  from  equations  (2-5)  and  (2-16)  to  be: 

6  =  3.31245248a  -  219.4468287  (6-73) 

Introducing  equation  (6-73)  into  equation  (6-72) 
and  after  some  manipulations  yields: 

s5  +  28s't  +  177s3+416.  249s  2  +  464.996s  +  264.  996 


a 


s8  +  52s7  +  926s6  +  6  700s5  +  20  8  2  5s1|  +  2  5000s3  +  5611.06  35s2 


-=-1  (6-74) 


-17555. 746s-  17555.746 


or  after  factoring  numerator  and  denominator  yields: 

a      (s2  +  1.479s  +  1.519)  (s+20.204)  (s  +  4  .  516)  (s+2) 

(s2  +  1.479s  +  1.5l9)  (s  +  20.204)  (s  +  4. 316)  (s+19.786) 


(s+5.642)    (s+1.4129)    (s-0.840) 


-=^1  (6-75) 
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which  is  in  agreement  with  the  conclusions  made  in  para- 
graph four,  i.e.,  the  number  of  the  fixed  roots  is  four 
which  is  also  the  order  of  the  characteristic  equation  of 
the  system  before  the  singular  cascade  compensator  was 
introduced  and  that  the  Root  locus  of  the  compensated 
singular  system  with  the  parameter  a  used  as  the  Root  locus 
variable  parameter  has  identical  zeros  and  the  same  number 
of  poles  as  the  uncompensated  system,  i.e.,  the  system 
before  the  singular  cascade  compensator  was  introduced. 

The  Root  locus  of  the  compensated  singular  sys- 
tem was  drawn  by  using  a  computer  program  in  which  equation 
(6-74)  was  implemented.   This  is  shown  in  figure  (6-31). 
In  this  figure  the  four  fixed  roots  are  shown  with  dark 
circles  at  the  corresponding  locations. 

a.   Range  of  Dominancy  of  the  Singular  Root 

By  the  use  of  a  DSL  simulation  program  in 

which  the  system  of  figure  (6-28b)  was  implemented,  the 

range  of  dominancy  of  the  singular  root  (c   =  0.60, 

to   =  1.232586)  was  investigated.   More  than  one  hundred 
ns  J  • 

properly  selected  parameter  pair  values  were  used  within 
the  system's  stable  region.   The  conclusion  from  this 
analysis  was  that  the  singular  root  has  a  very  wide  range 
of  dominancy.   The  detailed  results  which  were  interpreted 
on  both,  the  s -plane  and  the  a- 3  parameter  plane,  have  as 
follows: 
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(1)   S-plane.   The  Root  locus  of  the  com- 
pensated singular  system  is  shown  in  figure  (6-31) .   This 
Root  locus  has  been  drawn  with  the  parameter  a   used  as  the 
Root  locus  variable  parameter.   According  to  equation 
(6-72)  and  (6-75) ,  the  compensated  singular  system  has 
eight  characteristic  roots,  four  of  which  are  fixed  at  the 
points : 

s1  =  -4.316      s3  =  -0.7396+J0.9861 
s2  =  -20.204     s4  =  -0.7396-j0.986l 

i.e.,  two  fixed  real  roots  and  the  fixed  complex  pair  of 
the  singular  roots.   All  these  fixed  roots  are  shown  in  the 
corresponding  locations  of  figure  (6-31)  with  dark  circle 
marks.   The  location  of  each  of  the  other  four  character- 
istic roots  is  varying  on  the  corresponding  Root  locus 
segments  as  a  function  of  the  parameter  a  value.   It  must 
be  noticed  that  a  specific  value  of  the  parameter  a  defines 
a  unique  pair  of  parameter  a  and  $  values  since  for  each 
Value  of  a  corresponds  one  and  only  one  value  of  3  as  can 
be  seen  from  equation  (6-73),  which  is  the  equation  of  the 
operating  singular  line. 

Interpreting  the  computer  results  which  were 
obtained  from  a  DSL  simulation  and  a  polynomial  root  finding 
program,  the  following  remarks  and  conclusions  can  be 
stated  in  correlation  of  the  parameter  a  value  and  the 
Root  locus  of  figure  (6-31) . 
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(a)   0  <  a  <■  66.25.   Within  this  range  of 
the  parameter  a  value ,  the  compensated  singular  system  is 
unstable  because  there  is  a  real  root  on  the  right  half  of 
the  s-plane.   Based  on  the  need  for  a  specific  name  to  be 
given  to  this  root,  since  it  was  found  that  after  the 
singular  root  this  is  the  only  root  which  affects  the 
transient  response  of  the  compensated  singular  system  and 
the  fact  that  this  root  lies  on  the  same  Root  locus  segment 
with  the  "singular  root",  it  was  called  the  "non  singular 
root". 

Cb)   a  =  66.25.   For  this  value  of  the 
parameter  a  the  non  singular  root  lies  on  the  origin  of 
the  s-plane  and  the  compens-ated  singular  system  is  still 
unstable  or  conditionally  stable. 

(c)   66.25  <  a  <  81.30.   Within  this  range 
of  the  parameter  a  value,  the  non  singular  root  lies  on 
the  real  axis  in  the  left  half  of  the  s-plane.   The  compen- 
sated singular  system  is  stable  and  the  non  singular  root 
is  dominant.   In  figure  (6-32)  the  transient  response  of  the 
compensated  singular  system  is  shown  for  a  unit  step  input 
and  with  an  arbitrarily  selected  value  of  the  parameter  a 
within  this  range,  i.e.,  a  =  75  for  which  the  corresponding 
value  of  3  was  found  from  the  equation  of  the  operating 
singular  line,  (equation  £6-73))  to  be  g  =  28.987. 
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(d)  a  =  81.30:    For  this  value  of  the 
parameter  a  the  non  singular  root  i$  repeated  at  the 
point  where  the  Root  locus  breaks  away  from  the  axis 
of  reals,  i.e.,  at  this  point  the  non  singular  root 
meets  the  other  real  root  which  was  moving  to  the  right 
on  the  real  axis.   This  root  is  again  the  dominant  root. 
In  figure  (6-33)  the  transient  response  of  the  compensa- 
ted singular  system  is  shown  for  this  value  of  the 
parameter  a.   As  the  value  of  the  parameter  a   further 
increases,  these  two  roots  form  a  complex  pair  of  roots. 
In  this  case  by  the  name  "non  singular  root"  this  complex 
pair  will  be  meant. 

(e)  81.30  <  a    95.   Within  this  range 
of  the  parameter  a  value  the  non  singular  root  is 
complex  and  although  it  has  moved  away  from  the  real 

axis  it  is  still  closer  to  the  origin  and  to  the  imaginary 
axis  than  the  singular  root.   Within  this  range  the  con- 
tribution of  the  singular  root  to  the  transient  response 
of  the  compensated  singular  system  was  negligibly  small, 
and  the  non  singular  root  was  dominant.   In  figure  (6-34) 
the  transient  response  of  the  system  is  shown  for  a  =  90. 
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(f)   95  <:  a  <  105.   Within  this  range  of 
the  parameter  a  value  the  non  singular  root  has  ap- 
proached closer  to  the  singular  root.   The  contribution 
of  both  of  these  roots  to  the  transient  response  of  the 
compensated  singular  system  is  significant  and  none  of 
these  can  be  said  to  be  dominant. 

When  the  second  order  model  cor- 
responding to  the  singular  root  (t,      =    0.6,w    =  1.2326) 

was  considered,  with  transfer  function: 


ns 


1.5193 


2  2 

S  +  2c  <*>  „  S  +  0), 


s  ns 


s  +0.7396s+1.5193 


(6-76) 


ns 


the  corresponding  transient  response  characteristics  were 
calculated  and  they  were  found  to  be  as  follows: 


Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot 

Transient  oscillatory  frequence 


t 

p 

=    3,19    sec 

t 
s 

=    5.41    sec 

Mp 

=    9.48% 

u+ 

=0.99    rad/sec 
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Within  this  range  of  parameter  a   value  the  transient 
response  characteristics  of  the  compensated  singular  system 
were  found  to  be  as  follows: 


Time  to  peak  of  first  overshoot: 


Settling  time 


Percent  overshoot 


Transient  oscillatory  frequency 


5    sec . 


5    <:    ts     <    8    sec. 


8    s   Mp    <    111 


0 . 6    s    oi      <:    1    rad/sec 


Because  in  practical  applications  the 
transient  response  characteristics  which  are  of  the  basic 
interest  are  the  settling  time  and  the  percent  overshoot 
it  seem  from  the  above  results  that  an  acceptable  solution 
to  the  problem  can  be  derived  within  this  range  of  the 
parameter  a.   In  figure  (6-35)  the  transient  response  of 
the  system  is  shown  for  a  =  100. 

(g)   105  <  a  <  131.5.   Within  this  range  of 
the  parameter  a  value  the  singular  and  the  non  singular 
root  are  very  close  to  each  other.   The  contribution  of 
both  of  these  roots  to  the  transient  response  of  the  com- 
pensated singular  system  is  again  significant,  and  none 
of  these  can  be  called  dominant.   The  transient  response 
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was  in  general  characterized  by  long  settling  time,  which 
was  approximately  tg  =8.5  sec  and  high  percent  overshoot 
values,  which  were  found  to  lie  in  the  range  from  11%  up 
to  18.5%.   In  figure  (6-36)  the  transient  response  of  the 
compensated  singular  system  is  shown  for  a  =  120. 

(h)   a  =  131.5.   For  this  value  of  the 
parameter  a  the  non  singular  root  becomes  identical  with 
the  singular  root,  i.e.,  the  singular  root  is  repeated. 
In  figure  (6-37)  the  transient  response  of  the  system  is 
shown  for  this  value  of  the  parameter  a. 

(i)   131.5  <  a  <  285.   Within  this  range 
of  the  parameter  a  value  the  non  singular  root,  although 
it  is  further  away  from  the  origin  and  the  imaginary  axis 
than  the  singular  root,  its  contribution  to  the  transient 
response  of  the  compensated  singular  system  is  still 
significant  and  none  of  these  roots  can  be  called  dominant 
The  transient  response  characteristics  of  the  compensated 
singular  system  were  found  to  be  as  follows: 


Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot 

Transient  oscillatory  frequency 


2.8  < 


* 


<  4  sec 


7  <  t_  <  8.8  sec 

v     b    V 

18.96  <  Mp  <  23.6% 
0.98  <  "t  <  1.1  rad/sec 
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In  figure  (6-38)  the  transient  response  of  the  compensated 
singular  system  is  shown  for  a  =  200. 

(j)   285  <  a  <  385.   Within  this  range  of 
the  parameter  a  value  the  non  singular  root  has  moved 
further  from  the  imaginary  axis  and  from  the  origin,  and 
the  singular  root  has  become  the  dominant  root  but  with 
low  degree  of  dominancy.   Within  this  range,  the  transient 
response  characteristics  of  the  compensated  singular  system 
were  found  to  be  as  follows: 


Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot 

Transient  oscillatory  frequency 


v  2-8 

6.6  <  t„  <  7  sec 

s  v 

14.22  <  Mp  <  18.96% 
0.99   <  o)t  <1.1  rad/sec 


i.e.,  the  percent  overshoot  is  from  50%  up  to  100%  over  the 
corresponding  to  the  singular  root  second  order  model 
value,  and  the  settling  time  from  22%  up  to  30%.   In  figure 
(6-39)  the  transient  response  of  the  compensated  singular 

system  is  shown  for  a  =  330. 

'■-, 

(k)   385  <  a  <  1520.   Within  this  range 
of  the  parameter  a  value  the  non  singular  root  is  so  far 
away  from  the  origin  in  comparison  with  the  singular  root, 
that  its  contribution  to  the  transient  response  of  the 
compensated  singular  system  is  either  very  little  or 
neglegible  and  the  singular  root  is  the  dominant  root. 
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Within  this  range  of  the  parameter  a  value  the  transient 
response  characteristics  of  the  compensated  singular  sys- 
tem were  found  to  be  as  follows: 

Time  to  peak  of  first  overshoot:  2.8  <  t   <3.1  sec 

Settling  time  :  5.6  <  ts  <  6.6  sec 

Percent  overshoot  :  11.85.  $  Mp  $  14.22% 

Transient  oscillatory  frequency:  0.98  < u  <0.99  rad/sec 

i.e.,  the  percent  overshoot  from  10%  up  to  50%  over  the 
second  order  model  value  and  the  settling  time  no  more 
than  22%  over  that  corresponding  to  the  second  order  model 
value.   In  figure  (6-40)  through  (6~44)  the  transient 
response  of  the  compensated  singular  system  is  shown  for 
several  values  of  the  parameter  a  within  the  range  which 
was  considered.   It  was  noticed  in  this  case  that  although 
the  parameter  values  vary  within  a  wide  range,  the  transient 
response  of  the  compensated  singular  system  remains  almost 
unaffected. 

CI)   1520  <  a  <  1922.   Within  this  range 
of  the  parameter  a  value  the  non  singular  root  appears  again 
to  affect  appreciably  the  transient  response  of  the  com- 
pensated singular  system  due  to  the  fact  that  although  it 
is  far  away  from  the  origin  it  is  very  close  to  the  imagin- 
ary axis.   The  transient  response  is  characterized  with 
long  settling  time  and  high  frequency  of  oscillations.   In 
figure  (6-45)  through  (6-47)  the  transient  response  of  the 
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compensated  singular  system  is  shown  for  three  selected 
values  of  the  parameter  a  within  the  range  which  was 
considered. 

(m)   a  =  1922.   For  this  value  of  the 
parameter  a  the  non  singular  root  lies  again  on  the 
imaginary  axis  and  the  system  is  unstable  or  conditionally 
stable.   The  corresponding  value  of  w  for  the  non  singular 
root  at  that  point  is  w  =  8.67  rad/sec. 

(n)   a  >  1922.   Within  this  range  of  the 
parameter  a  value  the  non  singular  root  lies  on  the  right 
half  of  the  s-plane  and  the  compensated  singular  system 
is  unstable.   In  figure  (6-48)  the  transient  response  of 
the  system  is  shown  for  such  a  value  of  a, i.e.,    a=2000. 

The  conclusion  from  the  above  analysis  is 
that  there  is  in  general  a  wide  range  of  the  parameter  a 
values  for  which  the  singular  root  is  dominant,  and  which 
may  be  restricted  according  to  the  tolerance  of  the 
transient  response  characteristic  values  which  are  usually 
determined  with  problem  specifications. 

(2)   Parameter  Plane..   In  figure  C6-49)  the 
parameter  plane  diagram  of  the  compensated  singular  system 
is  shown  for  the  singular  values  of  Cg  and  oons  together 
with  the  singular  line  which  is  associated  with  them.   As 
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can  be  seen  in  this  figure  all  the  above  mentioned  lines 
are  tangent  at  the  same  point  of  the  parameter  plane. 
In  figure  (6-50)  the  same  parameter  plane 
diagram  is  shown  for  different  scale  values.   On  the  same 
diagram  the  compensated  singular  system  stable  region  and 
dominant  singular  root  region  have  been  interpreted  in 
terms  of  singular  line  segments. 

b.   Hardware  implementation  of  the  compensated 

singular  system 

According  to  Part  C  of  this  Section,  the  com- 
pensated singular  system  can  be  implemented  as  shown  in 
figure  (6-28a)  or  (6-28b)  or  (6-28c),  all  of  which  are 
equivalent. 

Choice  of  the  configuration  to  be  used  in 
practice  will  be  in  general  determined  from  the  specific 
application  and  the  type  of  compensator  available,  i.e., 
integrated  form  or  distinct  components. 

When  it  is  required  to  maintain  the  two  para- 
meters a  and  3  adjustable  or  when  the  compensator  of  the 
system  is  going  to  be  built  from  distinct  components  then 
the  configuration  of  figure  (6-28b)  seems  to  be  the  most 
representative  one  for  such  a  case. 

Considering  also  the  advance  of  the  technology 
today  in  integrated  circuits,  it  must  be  realized  that 
any  required  compensator  can  be  built  in  integrated  form. 
In  such  a  case  either  figure  ^6-28a)  or  figure  (6-28c) 
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can  be  implemented,  after  the  two  parameter  values  have 
been  determined. 

c.   Remarks 

When  the  singular  line  theory  is  considered  as 
a  design  method  for  compensation  of  linear  systems  then 
the  following  remarks  can  be  stated  concerning  the  singular 
case  which  was  considered  in  this  section. 

CI)   It  can  be  applied  to  any  type  of  plant. 

(2)  The  design  problem  has  been  simplified  and 
better  organized  since  in  general  the  basic  design  problem 
is  not  to  find  the  required  compensator  but   the  appropriate 
values  of  the  two  parameters  a  and  3  which  have  to  be  used 
with  the  singular  compensator  which  was  defined  by  equation 
(6-18)  . 

(3)  The  existence  of  the  singular  line  pro- 
vides increased  flexibility  in  the  design  procedure. 

It  must  also  be  noticed  that  when  the  two 
paramaters  a  and  3  are  considered  as  time  variant  quantities 
of  an  adaptive  system  and  its  compensator  respectively,  then 
by  adjusting  their  initial  value  and  their  rate  of  change 
in  such  a  way  that  the  operating  point  moves  along  the 
corresponding  to  the  singular  root  dominant  root  region, 
system  performance  will  be  essentially  unchanged.   In 
this  case  other  characteristic  roots  of  the  system  will 
change  but  since  they  are  not  dominant,  system  performance 
can  not  be  affected.   When  the  above  mentioned  rate 
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adjustments  are  made  automatically  the  compensated  singular 
system  becomes  self-adaptive.   This  type  of  application 
of  the  singular  line  theory  in  adaptive  systems  needs  to 
be  further  investigated. 

In  the  previous  example  the  results 
which  were  obtained  were  related  with  a  mathematical 
accuracy  which  in  general  can  not  be  achieved  in  practice. 
In  this  case  a  reasonable  question  will  be:   How  is  the 
performance  of  the  compensated  singular  system  affected 
when  the  values  of  the   two  parameters  a  and  3  are  set 
approximated  with  respect  to  their  theoretical  values? 
or  what  is  the  sensitivity  of  the  singular  dominant  root 
with  respect  to  each  or  both  of  the  two  parameter  values? 
This  problem  is  next  analyzed. 

G.   SENSITIVITY  ANALYSIS 

1.   Introduction 

The  characteristics  of  a  physical  component  are 
subject  to  change  for  various  reasons.   Perfect  control 
systems  with  stable  components  which  can  be  built  or  can 
maintain  their  design  values  under  any  conditions  are 
difficult  to  obtain  in  practice.   Take  for  example  such 
a  passive  component  as  a  resistor.   It  is  almost  impossible 
to  obtain  a  resistor  with  the  exact  design  resistance 
value.   For  this  reason  a  tolerance  usually  ranging  from 
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II  up  to  30%  will  accompany  its  specifications.   Even  if 
this  "ideal"  component  is  obtainable,  factors  like  tem- 
perature variations,  radiation  etc.,  may  eventually  cause 
the  resistance  value  to  vary  and  produce  undesirable 
results  on  the  control  system  performance.   These  general 
rules  are  also  applied  when  an  active  or  solid  state  de- 
vice is  considered. 

A  measure  of  the  effects  of  parameter  variation 
upon  the  performance  of  a  control  system  has  in  general 
been  demoted  as  sensitivity.   Basically  there  are  two 
typesof  parameter  variations  studied  by  various  inves- 
tigators.  The  first  one  which  is  the  most  popular  and 
the  most  common  in  practice  is  the  incremental  variation 
case,  where  the  change  in  parameter  values  is  assumed  to 
be  very  small,  i.e.,  inf inite-simal .   The  corresponding 
sensitivity  study  in  this  case  is  called  <<Microscopic 
Sensitivity  Analysis  >>.   The  second  case  is  that  of 
relative  large  variations  in  parameter  values  i.e., 
finite  changes.   In  this  case  the  corresponding  sensitivity 
study  is  called  <<  Microscopic  Sensitivity  Analysis  >>. 

The  definition  also  of  the  sensitivity  coefficient 
varies  among  investigators.   Generally  the  sensitivity 
coefficient  S  can  be  defined  as  the  ratio  of  the  amount  by 
which  a  control  system's  performance  characteristic  deviates 
from  its  nominal  (or  original)  value  when  one  or  more  sys- 
tem parameters  vary  to  its  nominal  value.   From  this 
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definition  the  sensitivity  coefficient  is  mathematically 
given  by: 

S  =  —  C6-77) 

c 

where  C  denotes  the  nominal  value  of  the  system  charac- 
teristic which  is  studied  and  Ac  the  deviation  from  this 
value  when  one  or  more  system  parameters  which  are  con- 
sidered vary. 

There  have  been  many  sensitivity  functions 
defined  for  the  effects  of  parameter  variations  on  the 
transfer  functions  of  linear  control  systems,  some  of 
which  are  the  following: 

a.  The  Kokotovic  Method 

The  root  sensitivity  functions  developed  by 
Kokotovic  and  Siljak  provide  most  general  solution  to  the 
sensitivity  problem  for  small  parameter  variations  in 
linear  control  systems  and  the  most  applicable  to  parameter 
plane  methods . 

b.  The  F.H.  Hollister  Method 

The  macroscopic  root  sensitivity  functions 
defined  by  F.H.  Hollister  are  applicable  for  large  para- 
meter variations. 

c.  Bode  Method 

Bode  defines  a  logarithmic  or  normalized 
sensitivity  function  of  a  variable,  x.  (t,  d.),  with 
respect  to  a  parameter   a.,  as: 
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3  In   a. 

SX1  =  L 

dj     9  Jin  x. 


In  general  an  analytic  investigation  of  the 
sensitivity  problem  of  a  linear  control  system  is  neither 
simple  nor  easy.   By  computing  sensitivity  coefficients 
limits  of  individual  parameter  variations  may  be  deduced 
and  tolerance  levels  can  be  formulated. 

2 .   Sensitivity  Analysis  of  a  Singular  System 

The  sensitivity  analysis  of  a  singular  system 
can  be  performed  by  any  method  applicable  to  a  linear 
feedback  control  system.   By  such  a  method  the  sensitivity 
of  any  system  performance  characteristic  due  to  variations 
of  one  or  more  parameter  values  can  be  analyzed.   Especially 
when  compensation  of  a  control  system  is  designed  by  the 
use  of  the  "singular  line  theory",  then  the  transient 
response  characteristics  and  the  two  variable  parameters  a 
and  3  are  of  the  main  interest  for  the  singular  system 
which  results,  at  least  for  the  first  steps  of  the  design 
process.   Therefore  it  will  be  very  important  to  know  how 
transient  response  characteristics  will  be  affected  and 
how  tolerance  levels  can  be  established  concerning  varia- 
tions of  the  two  parameter  a  and  $  values.   In  terms  of 
sensitivity  this  could  be  stated  as:   How  sensitive  are 
the  transient  response  characteristics  of  a  singular  system 
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to  variations  of  the  two  parameter  values  and  how  can  the 
corresponding  tolerance  levels  be  determined? 

The  goal  of  the  present  analysis  is  to  give  an 
answer  to  these  questions  by  a  relatively  simple  and 
accurate  way,  in  which  computer  use  will  be  optional. 

Generally  when  designing  by  the  use  of 
singular  line  theory,  specifications  on  transient  response 
characteristics  determine  the  desired  operating  singular 
point  (see  paragraph  VI,  F,  6).   From  this  singular  point 
or  singular  root  the  equation  of  the  corresponding  singular 
line  will  be  evaluated.   Some  other  specifications  or 
physical  limitations  of  the  system  may  also  restrict  the 
permissible  operating  region  on  the  singular  line  to  a 
finite  line  segment.   Then  when  the  operating  point  moves 
or  stays  fixed  on  this  singular  line  segment,  performance 
requirements  on  transient  response  characteristics  are 
satisfied.   The  design  process  up  to  this  point  has  been 
described  in  detail  in  the  previous  part  of  this  section. 

It  will  be  assumed  now,  that  the  operating  point 
moves  away  from  the  singular  line  due  to  variations  of 
either  one  or  both  of  a  and  3  parameter  values  from  their 
design  values.   This  will  be  in  general  the  case  for  an 
actual  system.   In  such  a  case  the  acceptable  tolerance 
on  the  transient  response  characteristics  will  determine 
how  far  away  from  the  singular  line  the  operating  point 
is  permitted  to  move. 
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Since  there  is  at  least  one  singular  line  for 
each  value  of  r,   for  a  singular  system  it  will  be  expected 
that  for  a  singular  system  there  is  an  infinite  number  of 
singular  lines  on  the  parameter  plane,  each  of  which 
corresponds  to  a  different  value  of  g.   Therefore  accord- 
ing to  the  tolerances,  which  are  given  for  the  transient 
response  characteristics  of  the  compensated  singular 
system,  a  permissible  region  can  be  established  on  the 
parameter  plane  for  the  operating  point,  which  will  be 
bounded  between  two  such  singular  lines.   This  concept 
not  only  was  applied  in  the  following  example  but  was 
also  further  investigated  from  a  theoretical  point  of  view 
and  it  was  found  that  it  is  true  and  valid,  which  will  be 
shown  later  in  the  example. 

Determination  of  these  two  singular  lines  on  the 
parameter  plane  will  not  only  give  the  answer  to  the  pre- 
vious question,  but  it  will  give  also  a  lot  of  other 
information  concerning  sensitivity  and  performance  of  a 
singular  system  for  parameter  variations.   It  will  also 
give  a  precise  way  for  tolerance  values  establishment  for 
a  singular  control  system. 

Although  different  specifications  or  different 
systems  may  give  rise  to  different  boundaries  for  the 
operating  point  on  the  parameter  plane  or  different  cri-.  \ 
teria  in  order  to  determine  these  boundaries,  the  procedure 
which  must  be  followed  will  be  in  general  the  same.   This 
will  be  illustrated  by  the  following  example. 
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3 .   Example  (6-7) 

Consider  again  the  plant  of  example  C6-6)  with 
transfer  function: 


GCs)  ■     20 


s(s  +  l)  (s  +  5) 

which  is  compensated  by  an  initial  cascade  compensator 
with  transfer  function: 

Gc(s)  -    'J£S- 

c       s  +  20 

The  Root  locus  of  the  resulting  system  which  is 
drawn  with  the  gain  of  the  plant  used,  as  the  Root  locus 
variable  parameter,  is  shown  in  figure  (6-51). 

It  is  already  known  (see  paragraph  VI,  F,  6) 
that  any  point  of  the  above  Root  locus,  which  does  not 
lie  on  the  real  axis  is  a  singular  point  for  the  compen- 
sated singular  system.  It  will  be  assumed  that  the 
specifications  on  the  transient  response  characteristics 
for  the  compensated  singular  system  require  the  operating 
point  to  be  the  Root  locus  point  which  is  located  at 
(5   =  0.60,  a)   =  1.23)  and  that  according  to  the  given 

►tolerance  on  the  transient  response  characteristics  which 
are  interpreted  in  terms  of  the  second  order  model,  this 
singular  point  is  permitted  to  move  within  the  finite  Root 
locus  segment  which  is  bounded  between  the  Root  locus 
points : 
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C  Cl   =  0.65,  w'   =  1.03) 
s  ns 


and 


(  5"   =  0.55,  a)"    =  1.58) 
v   s        '   ns        J 


These  points  can  be  determined  by  the  intersection 

of  the  Root  locus  with  the  corresponding  constant  £  lines 

(see  figure  (6-51))  or  can  be  evaluated  by  using  the 

"Singular  point  program"  in  the  computer.   In  the  second 

case,  although  more  precise  values  can  be  obtained  for 

to   '  and  w   "  it  will  be  seen  later  in  this  example  that 
ns       ns  r 

such  an  accuracy  is  not  actually  needed. 

Since  the  above  defined  Root  locus  segment  has 
been  established  by  interpreting  problem  specifications 
in  terms  of  second  order  model,  will  be  actually  the 
permissible  region  for  the  operating  singular  root,  if 
and  only  if,  this  singular  root  has  a  high  degree  of 
dominancy  anywhere  within  this  region. 

By,  either  studying  the  Root  locus  of  figure 
(6-51)  or  taking  into  account  the  results  obtained  from 
example  (6-6)  ,  a  basic  remark  concerning  the  dominancy  of 
the  singular  root  in  the  case  which  is  considered  is  that 
as   £   decreases,  the  degree  of  dominancy  of  the  singular 
root  also  decreases  because  it  moves  away  from  the  origin 
and  from  the  imaginary  axis,  and  the  contribution  of  the 
non-singular  root  (see  example  (6-6))  becomes  more 
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significant.   Based  on  this  remark  and  the  need  for  the 
given  specifications  on  the  transient  response  charac- 
teristics to  be  met,  it  was  decided  the  permissible 
region  for  the  operating  singular  root  locus  to  be  restric- 
ted between  the  points: 


(5   =  0.60,  uj   =  1.23) 
*•  s  ns 


and 


U'   =  0.65,  id'   =  1.03) 
v  s        '   ns       J 


which  are  shown  in  figure  (6-51) . 

At  this  point  it  must  be  noticed  that  such  a 
decision  is  based  absolutely  on  personal  judgments  the 
correctness  of  which  can  be  checked  by  simulation  of  the 
compensated  singular  system  at  the  limiting  acceptable 
values  for  the  two  parameters  a  and  6,  which  were  obtained 
from  this  analysis,  as  it  will  be  shown  later  in  this 
example. 

By  introducing  a  singular  cascade  compensator, 
the  characteristic  equation  of  the  compensated  singular 
system  becomes: 

s8+52s7+92  6s6+(6700+a)s5+(20  8  2  5+28a)s4+ 

(2  5000+17  7a)s3+(10000+350a+20B)s2+       (6-7  8) 

(200a+803)s  +  803  =  0 
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When  the  singular  point  (c   =  0.60,lo    =  1.23) 
&     f  ^s  »   ns       J 

was  considered,  the  equation  of  the  corresponding  singular 
line  was  found  from  equation  (2-5)  and  (2-16)  to  be: 

3   =  3.304a   -  218.355  (6-79) 

Comparison  of.  equations  (6-79)  and  (6-73)  reveals 

that  the  use  of  the  computer  for  evaluation  of  the  polar 

coordinates  of  a  singular  point  and  from  these  the  equation 

of  the  corresponding  singular  line  to  be  calculated,  does 

not  give  any  significant  improvement  on  the  accuracy  of 

the  equation  of  the  singular  line,  for  analysis  and  design 

purposes . 

When  the  singular  root  (£f  =  0.65,  w '   =  1.03) 
&  *•  s  ns       J 

was  considered  the  equation  of  the  corresponding  singular 
line  was  found  from  equation  (2-5)  and  (2-16)  to  be: 

3  =  2.444a   -  138.097  (6-80) 

The  above  two  singular  lines  were  drawn  on  the 
parameter  plane  diagram  which  is  shown  in  figure  (6-52). 
After  taking  into  account  the  results  of  example  (6-6)  ,  the 
characteristic  equation  of  the  compensated  singular  system, 
and  equations  (6-79)  and  (6-80)  of  the  boundary  singular 
lines,  the  stable  region  A  B  C  A  B  C  A,  and  the  region 
D  E  F  G  in  which  the  varying  singular  root  has  a  high 
degree  of  dominancy  were  approximately  established  on  the 
parameter  plane  diagram  of  figure  (6-52)  . 
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When  the  operating  point  lies  at  the  point  A 
where  the  two  singular  lines  intersect,  then  the  compen- 
sated singular  system  has  both  the  corresponding  singular 
roots  as  characteristic  roots.   When  the  operating  point 
lies  within  the  region  D  E  F  G  where  the  varying  singular 
root  is  dominant  with  a  high  degree  of  dominancy,  then 
the  specifications  on  the  transient  response  characteris- 
tics for  the  compensated  singular  system  are  met,  i.e: 


Peak  overshoot  Mp 

Settling  time  t 

&  s 

Trans,  oscil.  freq.  co 


From  6.8%  up  to  12.6% 

From  4.6  sec  up  to  5.9  sec 

From  0.78  rad/sec  up  to  1.32rad/sec 


This  was  actually  verified  in  the  computer  by  a 
DSL  simulation  program  of  the  compensated  singular  system, 
in  which  several  pairs  of  a  and  6  parameter  values  which 
define  points  in  the  region  D  E  F  G  of  the  parameter  plane 
diagram  of  figure  (6-52)  were  considered.   In  figure  (6-53) 
through  (6-56)  the  time  response  of  the  system  for  some 
of  these  points  is  shown. 

a.   Sensitivity  Study 

Each  operating  point  within  the  region  D  E  F  G 
of  figure  (6-52)  lies  on  a  singular  line  (it  will  be  shown 
later  in  this  example) ,  which  corresponds  to  a  singular 

root  (c  ,o)  )  .   These  polar  coordinates  of  this  singular 

v  s'  ns;         r  ° 

root  can  be  approximated  for  any  such  point  from  the 
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corresponding  values  of  the  boundary  singular  lines.   For 

example  the  point  M,  and  any  other  point  which  lies  half 

way  between  the  boundary  singular  lines  will  correspond 

to  a  singular  root  fc   ■  0.625,  u   =  1.13).  When  such  a 
&  v  s        '   ns       J 

point,  i.e.,  the  point  M,  (a  =  900,  8  =  2400)   is 
considered  the  corresponding  transient  response  charac- 
teristics can  be  evaluated  from  the  second  order  model 

which  corresponds  to  the  singular  root  (?   =  0.625,  to 

r  6  v  s        '   ns 

1.13). 

Considering  now  that  due  to  variations  of 

either  one  or  both  of  the  parameter  a  and  8  values  the 

operating  point  moves  to  another   location  within  the 

region  D  E  F  G  of  the  parameter  plane  of  figure  (6-52), 

then  the  corresponding  new  c,      and  w   values  can  be 

r     6      s      ns 

approximated  and  from  these  the  corresponding  transient 

response  characteristics.   Therefore  the  sensitivity 

coefficient  which  relates  the  transient  response 

characteristics  and  the  variations  of  the  two  parameter 

a  and   8  values  can  be  evaluated  and  from  these  sensitivity 

studies  can  be  made  as  it  is  shown  below: 

(1)   Sensitivity  of  the  Peak  Overshoot  due  to 

Variations  of  the  Parameter  a.   Let  the 

operating  point  move  from  NL  (a  =  900,  8  =  2400)  to  M- 

(a  =  990,  8=  2400)  due  to  a  +10%  change  of  the  parameter  a 

value.   The  corresponding  singular  points  were  approximated 

as  (c   =  0.625,  w   =  1.13)  and  (c  '   =  0.64,  w'   =  1.07) 
v^s        '   ns  s        '   ns 
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respectively.   The  corresponding  percent  peak  overshoots 
as  they  can  be  found  from  the  second  order  model  are 
Mp  =  8%  and  Mp   =7.3%.   Therefore  according  to  the  defi- 
nition of  the  sensitivity  coefficient  which  was  given  in 
paragraph  one: 


£   -  ^F2  "  JL¥S-  '    -°-0875 


(2)   Sensitivity  of  the  Peak  Overshoot  due 
to  Variations  of  the  Parameter  g.   Let 
now  the  operating  point  more  from  M,  (  a  =  900,  g  =  2400) 
to  M3(a  =  900,  g  =  2640)  due  to  a  +  10%  change  of  the  para- 
meter g  value.   The  corresponding  singular  points  were 
approximated  as  (z;   =  0.625,  u   =  1.13)  and  (c;   =  0.61, 
u>"   =  1.19)  respectively.   The  corresponding  percent 
overshoot  was  evaluated  from  the  second  order  model  to  be 
Mp  =  8%  and  Mp"  =  8.9%.   Therefore  the  sensitivity  co- 
efficient which  relates  the  peak  overshoot  and  the  +101 
variation  of  the  parameter  g  value  from  its  value  at  the 
initial  point  Mj_  will  be: 


S^P=  4    -  %  -  8'9-8     .  0.1125 
g      Mp        8 


Comparison  of  the  above  evaluated  two  sensi- 
tivity coefficients  reveals  that  when  the  operating  point 
Mi  is  considered,  then  the  peak  overshoot  of  the  compensated 
singular  system  is  more  sensitive  to  variations  of  the 
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parameter  6  than  it  is  to  variations  of  the  parameter  a 
(almost  1.28  times  more  sensitive),  and  that  the  effect 
of  the  variations  of  these  two  parameters  on  the  Peak  over- 
shoot of  the  compensated  singular  system  is  opposite. 

It  must  be  noticed  that  although  the  values 
obtained  for  the  sensitivity  coefficients  are  approximate, 
the  conclusions  which  can  be  derived  from  a  relative  com- 
parison must  be  considered  as  very  accurate  for  analysis 
and  design  purposes. 

In  the  same  way  any  sensitivity  problem  which 
relates  the  transient  response  characteristics  of  the  com- 
pensated singular  system  and  variations  of  either  one  or 
both  of  the  parameter  a  and   3  values  can  be  studied. 

b.   Selection  of  the  Operating  Point 

In  general  the  operating  point  M  of  a  compen- 
sated singular  system  must  be  selected  where  the  singular 
root  has  the  highest  degree  of  dominancy  [which  results  in 
the  best  possible  approximation  of  the  transient  response 
characteristics  of  the  compensated  singular  system  by  the 
corresponding  second  order  model)  and  on  a  point  halfway 
between  the  boundary  singular  lines  (which  will  give  the 
maximum  tolerance  for  variations  of  the  two  parameters  a 
and  3,  when  it  is  considered  that  both  can  vary  simul- 
taneously) .   Under  these  considerations  the  best  location 
for  the  operating  point  in  the  case  which  was  studied  in 
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this  example  seems  to  be  at  the  point  MCa  =  1310  ,  3  = 
3500) ,  which  is  shown  in  figure  (6-52) ,   It  must  be  noticed 
that  some  other  specifications  of  system's  performance,  or 
other  limitations  may  require  another  location  for  the 
operating  point  M  to  be  determined  within  the  region 
D  E  F  G. 

c.   Tolerance  Levels  Establishment  on  Parameter  a 
and  3  Variations 

After  the  operating  point  M  has  been  chosen, 
then  tolerance  levels  on  the  two  parameters  a   and  3 
variations  must  be  established.   These  tolerance  levels 
must  be  determined  in  such  a  way  that  the  operating  point 
will  be  restricted  to  lie  within  a  circle  which  has  center 
the  point  M  itself  and  is  tangent  to  the  two  boundary 
singular  lines  and  to  any  other  boundary  line  in  the 
neighborhood  of  this  point.   For  the  case  of  the  example 
which  was  studied  the  required  tolerance  levels  for  both 
the  a  and  3  parameters  are  ±7.51  from  the  corresponding 
nominal  values  of  a  =  1310  and  3  =  3500,  as  it  can  be 
evaluated  from  figure  (6-52),  i.e.,  for  no  more  than  7.5% 
variations  towards  both  directions  of  each  parameter  a 
and  3  value,  the  operating  point  will  lie  within  the  circle 
shown  in  figure  (6-52)  whose  center  is  at  the  point  M, 
and  the  performance  of  the  compensated  singular  system  will 
be  such  that  the  specifications  will  be  met. 
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4.   Comments 

In  establishing  tolerance  levels  in  practice  some 
trade  offs  can  be  made  by  considering  each  parameter 
separately  since  the  nature  of  each  is  different,  i.e.,  the 
parameter  3  is  always  associated  with  a  gain  value  while 
the  parameter  a  is  associated  either  with  another  gain 
value  or  with  pole  locations  (see  figure  (6-28a  ,b,c) ) .   If 
for  example,  in  the  previous  case,  which  was  studied  above, 
the  poles  of  the  singular  compensator  can  be  set  very 
accurately,  due  to  the  nature  of  the  singular  compensator 
which  is  going  to  be  used,  then  its  gain  3  will  not  be 
required  to  be  set  with  a  higher  than  ±151  accuracy,  in 
order  for  the  operating  point  to  be  restricted  within  the 
desired  circle  which  was  mentioned  above,  as  can  be  seen 
from  figure  (6-52). 

Another  general  remark  which  was  derived  from 
figure  (6-52)  is  that  as  far  away  from  the  intersection  of 
the  two  boundary  singular  lines  the  operating  point  M 
is  located,  the  less  sensitive  the  compensated  singular 
system  becomes  to  the  parameter  a  and  3  variations. 

The  above  comments  together  with  the  results 
obtained  up  to  this  point  from  the  sensitivity  analysis 
reveal  that  design  of  compensation  by  using  singular  line 
theory  leads  to  very  accurate  setting  of  tolerance  levels 
for  the  two  parameter  a  and  3  values  and  very  easy  and 
accurate  way  of  determining  the  optimal  operating  point  for 
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the  compensated  singular  system  either  on  the  s -plane  or 
on  the  a-8  plane. 

Although  all  the  required  information   in  order  to 
make  a  conclusion  about  the  sensitivity  of  a  compensated 
singular  system  transient  response  characteristics,  due  to 
variations  of  the  parameter  a  and  B  values,  have  not  been 
derived,  based  on  the  results  of  the  analysis  which  was 
performed  up  to  this  point,  it  is  believed  that  in  general 
compensation  of  a  system  by  using  singular  line  theory 
may  lead  to  control  systems  which  are  not  sensitive  or  at 
least  very  sensitive  to  parameter  values  variations. 

In  the  previous  example  it  was  shown  that  when 
the  operating  point  moves  within  the  region  D  E  F  G  of 
figure  (6-52)  then  the  transient  response  characteristics 
of  the  compensated  singular  system  can  be  well  approximated 
by  the  corresponding  second  order  model.   Actually  there  is 
no  question  about  it,  when  the  operating  point  lies  on  a 
singular  line,  like  the  two  boundary  singular  lines  D  E 
and  G  F  of  the  previous  mentioned  operating  region.   In  the 
case  where  the  operating  point  lies  between  these  boundaries, 
although  the  results  from  example  (6-7)  reveal  that  the 
method  is  still  valid,  some  reasonable  questions  which  may 
rise,  in  this  case  from  a  theoretical  point  of  view  are: 
Is  the  compensated  system  still  singular?   Do  singular 
conditions  exist?   Does  the  operating  point  still  lie  on 
a  singular  line?   or  have  all  the  characteristics  of  the 
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singular  conditions  been  lost?   These  questions  were 
answered  as  follows: 

When  a  linear  control  system  with  negative  unity 
feedback,  whose  open  loop  transfer  function  is  GCs)  = 
kF  (s)  is  compensated  by  a  singular  cascade  compensator,  ' 
the  characteristic  equation  of  the  compensated  singular 
system  will  be: 


1  +a  F  (s)  +B  kFg2(s)  =  0 


It  has  been  shown  (see  paragraph  III,  C,  3)  that  such  a 
system  has  in  general  an  infinite  number  of  singular  points, 
i.e.,  points  of  the  s-plane  which  satisfy  equation: 


L 


Fg(s)  =  nn 


where  n  is  an  real  integer  or  zero.   When  such  a  singular 
point  has  been  specified,  the  equation  of  the  corresponding 
singular  line  can  be  determined  either  by  direct  substitu- 
tion of  the  value  of  the  singular  point  s  in  the  character- 
istic equation  of  the  compensated  singular  system  or  by 
using  equations  (2-5)  and  (2-16).   From  the  equation  of 
this  singular  line  then,  an  infinite  number  of  a  and  8 
pair  values  corresponding  to  the  singular  point  under 
consideration,  can  be  determined.   In  the  opposite  case 
where  a  pair  of  a  and  g  values  has  been  specified  the 
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equation  of  the  corresponding  singular  line  and  correspon- 
ding to  this  equation  the  singular  point  can  also  be 
determined. 

This  will  be  illustrated  by  considering  the  follow- 
ing two  cases  from  example  (6-7)  where  the  operating  point 
moves  within  the  region  D  E  F  G.   In  the  first  case  the 
operating  point  arbitrarily  was  assumed  to  lie  at  the  point 
M,  (a=900,  6  =  2400)  of  figure  (6-52).   Then  the  characteristic 
equation  (6-78)  of  the  compensated  singular  system  becomes: 


s8+ 52s 7+9  26s6+ 7600s 5+46 02 5s4+18430  0s3+ 3730 0  0s2 
+372000S+192000  =  0 


By  the  use  of  a  computer  program  with  the  ZPOLR 
subroutine  the  roots  of  the  above  equation  were  found  to 
be : 

s±   =    -  0.7044558  +  jO. 8843676      s2  =  -  0  .  7044558- j 0 . 8843676 
s3  =  -  0.8651241  +  J5. 8923703      s4  =  -  0  .  8651241- j 5 . 8923703 

s5  =  -  22.1206845      ,         s6  =  "  20-1762447 
sy  =  -  4.4148436  sg  =  -  2.1490673 

Each  of  the  above  complex  roots  was  checked  and 
it  was  found  that  each  lies  on  the  Root  locus  of  figure 
(6-51).   Therefore  each  of  these  complex  roots  is  also  a 
singular  root,  (see  statement  after  equation  (6-61)).   A 
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more  general  theoretical  proof  of  this  statement  which 
is  always  true,  can  be  derived  from  paragraph  VI,  F,  4. 

Because  the  operating  region  of  the  parameter 
plane,  where  the  operating  point  M,  lies  has  been  designed 
as  dominant  root  region  only  the  first  singular  root,  i.e., 
the  singular  root  s,,„,  which  also  was  dominant,  was  further 
considered.   This  root  is: 


S1  2  =  -  0.7044558  ±j 


0.8843676 


or 


(C   =  0.623,   u    =  1.131) 

v  s        '    ns        J 


When  this  singular  root  was  considered,  from 
equation  (2-5)  and  (2-16) ,  the  equation  of  the  corresponding 
singular  line  was  found  to  be: 

3  =  2.8472a  -   162.1855  (6-81) 

Since  the  coordinates  of  M, ,  i.e.,  a  =  900  and 
3  =  2400  satisfy  the  above  equation  it  implies  that  this 
point  lies  on  a  singular  line  of  a  dominant  singular  root 

whose  coordinates  are  (s   =  0.623,  w   =  1.131). 

*•  s        '   ns        ' 

Two  points  are  interesting  here.   First,  the  above 
coordinates  of  the  dominant  singular  root  were  approximated 
from  the  corresponding  values  of  the  two  boundary  singular 
lines  of  figure  (6-52),  during  the  sensitivity  study  of 
example  (6-7),  (see  paragraph  G3  a  (1)  of  this  section), 
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and  second,  when  the  other  singular  root  was  considered, 
then  by  the  same  method,  it  was  found  that  the  operating 
point  M  lies  also  on  the  singular  line  corresponding  to 
this  root,  i.e.,  finally  the  operating  point  M,  lies  on 
the  intersection  of  two  singular  lines.  Of  course  under 
dominant  root  considerations  this  last  singular  line  is 
of  no  interest. 

It  was  next  assumed  that  the  operating  point  was 
moved  to  a  new  location  Mr  Cot  =  1430,8  =  4280)  of  the 
dominant  root  region  D  E  F  G  of  figure  (6-52).   Then  the 
characteristic  equation  of  the  compensated  singular 
system  was  evaluated  to  be: 


s8+5  2s7+926s6+8130s5+603  2  5s4+2  78110s3+596100s2 
+628400  s  +'342400  =  0 


By  the  same  method  as  in  the  previous  case  the  roots 
of  this  equation  were  found  to  be: 

S1  =  -  0.7257602  +  jO. 9470602  s2  =  -  0  .  7257602 -j 0 . 9470602 

s3  -  -  0.3981063  +  J7. 5192493  s4  =  -  0  .  3981063- j 7 . 5192493 

sc  =  -  23.1172585  s,  =  -  4.3550711 

s?  =  -  20.1934085  sg  =  -  2.0865288 

Each  of  the  above  complex  roots  was  also  checked 
and  it  was  found  that  it  lies  on  the  Root  locus  of 
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figure  (6-51).   Therefore  they  are  also  singular  roots 
In  this  case  the  dominant  singular  root  is: 

S,  -  -  -  0.7257602  ±j0. 9470602 

or 


(C   =  0.608,  u    =  1.193) 

*•  s        '   ns        J 


When  this  singular  root  was  considered,  from 
equation  (2-5)  and  (2-16),  the  equation  of  the  corresponding 
singular  line  was  found  to  be: 

3   =  3.1307a  -  196.1499  (6-82) 

Since  the  coordinates  of  the  operating  point  M_ 

which  was  considered,  i.e.,  a  =  1430  and  3  =  4280,  satisfy 

the  above  equation,  it  implies  that  the  operating  point 

Mr  lies  on  the  singular  line  which  is  defined  by  equation 

(6-82) .   It  must  also  be  noticed  that  equivalent  remarks 

like  these  which  were  stated  in  the  previous  case  about  the 

other  complex  root  s,  .    also  were  found  to  apply  and  in 

this  case. 

By  the  same  way  the  singular  lines  and  singular 

points  corresponding  to  any  point  within  the  region  D  E  F  G 

or  to  any  other  point  of  the  parameter  plane,  can  be 

evaluated. 

The  conclusion  from  these  results  is  that  indepen- 
dently of  where  the  operating  point  lies  on  the  a- 3  plane, 
singular  conditions  always  exist. 
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The  results  from  the  example  which  was  studied  are 
summarized  in  Table  I  of  Appendix  G. 

# 

H.   STEADY  STATE  ACCURACY  OF  A  SINGULAR  SYSTEM 

1 .  Introduction 

A  measure  of  control  system  accuracy  is  usually 
the  steady  state  error.   Related  concepts  and  definitions 
can  be  found  in  reference  [2] . 

When  steady  state  error  considerations  are  made 
for  a  singular  compensated  system  shown  in  figure  (6-5), 
the  question  which  arises  is:   What  is  the  effect  of  the 
singular  cascade  compensator  on  the  error  constant  K 
(which  is  also  called  the  DC  gain  or  zero-frequency  gain) , 
of  the  system  which  is  compensated  by  a  singular  cascade 
compensator? 

The  goal  of  the  present  study  is  for  an  answer  to 
be  given  to  this  question  and  for  the  required  design 
procedure  to  be  established. 

2 .  Steady  State  Error  Analysis  of  a  Singular  System 

It  is  known  that  if  a  finite  error  exists  in  steady 
state  operation  of  a  control  system  its  magnitude  is  deter- 
mined by  the  reciprocal  of  the  error  constant  K  . 
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When  a  unity  feedback  control  system  is  considered 
with  open  loop  transfer  function: 

(s  +  z    )  (s  +  z    )    ...    (s  +  z  0 
G(s)    =   k   F    (s)    =   k 


g  sN(s+PiKs+P2)     ...    (s+pn) 


then  the  error  constant  K  is 

A. 


or 


K   =  Aim  sN  G(s)  (6-83) 

x   s+o 


K   =  k  ZiS",zm  (6-84) 

PXP2  ...Pn 


When  the  above  considered  system  is  compensated 
by  a  singular  cascade  compensator  whose  transfer  function 
G   (s)  is  given  by  equation  (6-18),  then  the  error  constant 
of  the  compensated  singular  system  is: 


K        =    £im      sN   G      (s)    G(s)  (6-85) 

sx        s->o  SC 


Kp     Z     Z      ...      Z 


12  Jm        _      3 


K        =      l—l - —  =   -*-  K  (6-86) 

sx  ax 

a  p  p    ...    p 
i    2  n 
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i.e.,  the  error  constant  of  the  compensated  singular 
system  differs  from  the  error  constant  of  the  system 
before  the  singular  cascade  compensator  has  been  introduced 
by  a  factor  equal  to  the  ratio  3/a  of  the  two  parameter 
values.   Since  an  operating  singular  line  does  not  in 
general  pass  through  the  origin  of  the  a-3  plane,  the 
ratio  3/a  will  not  in  general  remain  constant  when  the 
operating  point  is  moved  along  a  singular  line. 

When  an  operating  region  on  the  a-3  plane  is 
considered,  which  actually  will  be  the  case,  i.e.,  the 
parameter  values  will  be   fixed,  then  the  error  constant 
of  the  compensated  singular  system  will  have  also  a  fixed 
value . 

Since  the  ratio  3/a  corresponds  to  the  slope  of  a 
straight  line  through  the  origin  of  the  a-3  plane,  it 
implies  that  the  locus  of  all  the  points  of  the  a-3  plane 
with  a  fixed  value  of  the  ratio  3/a  is  a  straight  line 
through  the  origin,  and  that  the  locus  of  all  the  points 
of  the  a-3  plane  with  a  specified  range  of  values  of  the 
ratio  3/a  is  a  region  of  the  a-3  plane  which  is  bounded 
between  two  such  straight  lines  through  the  origin. 

Therefore  interpretation  of  any  steady  state 
accuracy  specifications  in  terms  of  the  3/a  ratio  value 
leads  to  "steady  state  error  boundary  lines"  establishment 
on  the  a-3  plane  as  is  shown  in  the  following  example. 
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3.   Example  (6-8) 

In  this  example  it  was  assumed  that  steady  state 
accuracy  specifications  were  established  for  the  compensated 
singular  system  of  the  previous  example  (6-7),  and  that 
interpretation  of  these  in  terms  of  the  error  constant 
K   set  the  following  restriction  on  the  g/a  ratio  value, 
i.e.  : 


-2-  >  1  (6-87) 


The  corresponding  "steady  state  error  boundary 
line" for  this  case  is  the  straight  line  through  the  origin 
with  slope  equals  to  one.   This  was  drawn  on  the  a-g  plane 
shown  in  figure  (6-52) .   Then  all  the  points  of  the  a-g 
plane  above  this  line  satisfy  steady  state  accuracy 
specifications  of  the  problem.   Therefore  any  point  of  the 
dominant  singular  root  region  D  E  F  G  shown  in  figure  (6-52) 
will  also  give  an  acceptable  steady  state  error  constant 
for  the  compensated  singular  system. 

In  the  same  way,  when  another  system  or  different 
specifications  are  given  concerning  the  steady  state 
accuracy  establishment  of  the  corresponding  "steady  state 
error  boundary  lines"  or  line,  on  the  parameter  plane  the 
required  information  for  analysis  and  design  purpose  can 
be  derived. 


146 


VII.    THIRD  SINGULAR  CASE  ANALYSIS 


In  this  section  the  "third  singular  case"  which  was 
defined  in  Section  IV,  was  considered  for  analysis.   This 
singular  case  corresponds  to  the  compensation  structure  of 
figure  (4-5)  where: 

G(s)  =  kFg(s) 

Gc(s)  =  3s2Fg(s)  (7-1) 

H(s)  =  as 
and  the  plant  G(s)  =   kFg(s)  has  at  least  two  poles  at 
the  origin,  i.e.,  the  cascade  compensator  has  no  zero's 
at  the  origin. 

The  characteristic  equation  of  the  compensated 
singular  system  is: 

1  +  aks  F  (s)  +  3k  (sF  (s)2)  =  0 

A.   SYSTEM  OUTPUT  AND  ERROR  ANALYSIS 

The  transfer  function  of  the  compensated  singular 
system  is: 

Cfsl                         Bks2F   2(s) 
T(s)=  idii-  =   S (7-2) 

R(s)  1+aks   F    (s)+Bks2    F   2(s) 
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and  the  output  and  error  are  respectively; 


C(s)  =  R(s)  TCs) 


and 


E(s)  -  R(s)  {1  -  T(s)} 


(7-3) 


(7-4) 


1.   Plant  Type  Two 


When  the  plant  is  type  two,  i.e.: 


(s  +  z  )  (s  +  z  ).  .  .  (s  +  z  ) 

G(s)  =  kF  (s)  =  k— i 2- ^ 

s  (s+Pl)  (s+p2) .. . (s+pn) 

equations  (7-3)  and  (7-4)  become  respectively: 

m  2 

gk      [~](s  +  zj) 

j  =  i 
C(s)    = 


(7-5) 


-R(s) 


k=i 


3  =  1 


k=i 


s      _l    (s+Fk)    +aks  I      |(s  +  Zj)    [     |(s+pk)+3k  [[(s  +  zj) 

(7-6) 


3  =  i 


and 


E(s)    = 


n  m  n 

I     Ks+piJ    +aks    I      |(s  +  zi  )         Its+pk) 
k=i j  =  i  k=i 


n 


n 


-R(s) 


|(s+pk)    +aks   I      j(s  +  Zj)    I      j(s+pk)+3k  j      !    (s  +  zj) 
j=i  k=i  j=i 


k=i 


(7-7) 
When  the  final  value  theorem  was  applied  to  the 
above  equations  the  steady  state  output  and  error  were 
evaluated  assuming  that  all  the  poles  of  the  transfer  func- 
tion T(s)  lie  on  the  left  half  of  the  s-plane. 
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a.   Step  Input 

When  a  step  input  R(t)  ■  AuCt)  or  R[s)  =  A/s 
was  considered  it  was  found  that: 


and 


C(t) 

t-»-cc 


E(t) 


CssCt)  «  A 


E   Ct)  «  0 

ss  v  J 


b.   Ramp  Input 

When  a  ramp  input  RCt)  =  BtuCt)  or  R(s)  =  B/s' 
was  considered  it  was  found  that: 


and 


COO  =  C   ft)  -« 


n 


Ba 


ECt)  ■  E   (t)  = 

t-><x 


C=l  Pk 


m 


k=i   J 


c.   Parabolic  Input 

When  a  parabolic  input  R(t)  =  Ct  u(t)  or 

3 
R(s)  =  2C/s   was  considered  it  was  found  that: 


and 


C(t) 


E(t) 

t+cc 


Css(t)      ■     « 


EssCt)    =   « 
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2.   Plant  Type  Three 

In  the  same  way,  the  steady  state  output  and  error 
value  for  a  type  three  plant  were  found  to  be  as  follows: 

a.   Step  Input 

C(t)  =  C   ft)  »   A 


and 


ECt)  -  E   ft)  m    0 
t-»*      *b 


b.   Ramp  Input 


cct)  =  cssct)  = 


and 


ECt)  =  Ess(t)  = 

t->  = 


n 

BaHpk 

k=i  K 


m 

j«l  -> 


c.   Parabolic  Input 


C(t)  =  C   (t)  = 

t->a: 


and 


E(t)  =  E   (t) 

t+cc        " 
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3.   Remarks 

The  general  form  of  the  characteristic  equation 
of  the  compensated  singular  system  in  the  case  which  was 
studied  in  this  section  is; 

1+aks  F,Cs)  +  ek{s  F  Cs))  =  0  (7-8) 

where  the  transfer  function  of  the  plant  is: 

G(s)  -  kFg  (s)  (7-9) 

By  defining: 

F(s)  =  sFCT(s) 
the  characteristic  equation  (7-8)  becomes: 

l+akF(s)  +  3k  F2(s)  =  0  (7-10) 

which  is  the  characteristic  equation  of  a  singular  system. 

Comparison  of  equations  (7-10)  and  (6-16)  reveals 
that  the  design  process  which  has  to  be  followed  when 
compensation  is  performed  by  the  compensation  scheme, 
which  corresponds  to  the  "third  singular  case"  will  be  in 
general  the  same  with  that  corresponding  to  the  "second 
singular  case".   The  only  difference  will  be  that  instead 
of  using  the  Root  locus  of  the  uncompensated  system  as  the 
basic  design  tool,  the  Root  locus  of  the  uncompensated 
system  with  one  pole  less  than  the  actual  number  of  poles 
at  the  origin  has  to  be  used.   This  will  be  illustrated  in 
the  following  example. 
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B.    EXAMPLE  (7-1) 


In  this  example  the  plant  shown  in  figure  (7-1)  was 
considered  with  open  loop  transfer  function: 

1 


GCs)  =  kF  (s)  -  100  -2 


s  Cs+2) 


(7-11) 


This  plant  is  unstable  and  it  was  assumed  that  prob- 
lem specifications  dictate  for  a  stable  system  with  a 
dominant  root  with  polar  coordinates: 

C  =  0.40 
and 

2<w  <3 

The  compensated  singular  system  is  shown  in  figure 
(7-2).   As  can  be  seen  the  whole  design  problem  consists 
of  finding  the  appropriate  pair  of  a  and  g  parameter 
values  by  which  the  given  specifications  can  be  met.   This 
was  performed  in  three  steps  which  follow: 


1 .   Step  One 

G 

The  first  step  was  to  determine  the  required 
operating  singular  root  for  the  compensated  singular  sys- 
tem. 

The  Root  locus  of  the  uncompensated  system  with  a 
pole  less  than  the  actual  number  of  poles  at  the  origin 
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is  shown  in  figure  (7-3),  i.e.,  this  is  the  Root  locus  of 
equation: 


s(s  +  2) 


-B--1 


(7-12) 


By  the  same  reasoning  as  in  the  case  of  the  "second 
singular  case"  it  can  be  shown  that  any  point  of  the  above 
Root  locus,  which  does  not  lie  on  the  real  axis  is  a  sin- 
gular point  of  the  compensated  singular  system.   (See  Section 
VI,  F,  4,  b) .   Considering  next  the  specifications  the 
desired  singular  point  was  located  on  the  Root  locus  at  the 
point  C?  s  0.40,  ajn  =  2.5)  . 


2.  Step  Two 

The  second  step  was  to  find  the  equation  of  the 
singular  line,  corresponding  to  the  selected  singular 
point. 

From  equation  (7-9)  the  characteristic  equation 
of  the  compensated  singular  system  was  found  to  be: 

s4+4s3+(4+100a)s2+200aS+100g  =  0  (7-13) 


When  the  singular  root  (?_-  0.40, 


ns 


■  2.5)  of  the 


above  characteristic  equation  was  considered,  the  equation 
of  the  corresponding  singular  line  was  found  from  equation 
(2-5)  and  (2-16)  to  be; 


3  =  6. 25a-Q. 390625 


(7-14) 
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3.   Step  Three 

The  third  step  which  was  the  last  one  was  to  find 
the  required  operating  point  of  the  singular  line  which 
was  defined  by  equation  (7-14) . 

Because  it  was  also  desired  that  the  design  pro- 
cess for  this  problem  be  completed  without  using  a  computer 
this  step  was  performed  as  follows. 

Introducing  equation  (7-14)  into  (7-13)  the 
characteristic  equation  of  the  compensated  singular  system 
becomes : 


s4+4s3+(4+100a)s2+200aS+625a-39.0625  =  0 


(7-15) 


Although  according  to  the  results  which  were  stated 
in  Section  VI,  F,  4,  the  Root  locus  of  the  compensated 
singular  system,  i.e.,  the  Root  locus  of  equation  (7-15), 
with  the  parameter  a  used  as  the  Root  locus  variable  para- 
meter has  the  shape  of  the  Root  locus  of  equation  (7-12) 
shown  in  figure  (7-3),  and  therefore  the  drawing  of  the 
Root  locus  of  equation  (7-15)  is  not  actually  required,  it 
was  done  in  this  example  for  a  more  detailed  illustration. 

Equation  (7-15)  can  be  written  as: 


100s  +200S+625 
s4+4s3+4s2-39.0625 


=  -1 
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or  after  factoring  numerator  and  denominator; 


a- 


{s+(l  +  j2.291)}  {S+C1-J2.291)  } 


{s+(l  +  j2.291)  }  (S+C1-J2.291)}  {s+3.693}  {s-1.693} 


-s-1 


The  Root  locus  of  the  above  equation  is  shown  in 
figure  (7-4) .   This  root  locus  which  was  drawn  by  hand, 
consists  of  a  fixed  point  there  where  the  singular  root 
C?s=0 . 40,  wns=2.5)  or  s=-l±j2.291  is  located  and  the 
regular  Root  locus  segments. 

Comparison  of  figures  (7-3)  and  (7-4)  reveals  that 
the  two  Root  loci  have  actually  the  same  shape. 

Because  the  singular  root  and  the  complex  non- 
singular  root  (See  Section  VI,  F, 4, a, (1) , (c) ) ,  have  always 
equal  real  parts,  in  this  case,  in  order  for  the  singular 
root  to  be  dominant  its  residue  must  be  much  greater  than 
that  corresponding  to  the  non-singular  root,  i.e.,  the 
singular  root  must  lie  much  closer  to  the  origin  of  the 
s-plane  than  the  non-singular  root.   Therefore  the  non- 
singular  root  must  lie  far  away  above  the  singular  root 
which  is  shown  on  the  Root  locus  either  of  figure  (7-3)  or 
(7-4). 

By  arbitrarily  selecting  the  distance  of  the  non- 
singular  root  from  the  origin  to  be  relative  very  large, 
in  comparison  with  that  which  corresponds  to  the  singular 
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root,  dominance  of  the  singular  root  in  general  shall  be 
expected,  e.g.,  let  the  non-singular  root  be: 

s  =  -1  +  j  50 
Then  the  left  part  of  equation  (7-15)  must  be  equal  to  the 
product  of  the  singular  root  and  non-singular  root.   This 
yields  a  =  25.07.   From  equation  (7-14)  the  corresponding 
value  of  the  parameter  g,  was  found  to  be  $  =  156.31. 

Therefore  the  required  parameter  pair  values  or  the 
desired  operating  point  of  the  a-3  plane  is  (a=25.07, 
3=156.31) . 

4.   Remarks 

The  compensated  singular  system  shown  in  figure 
(7-2)  was  simulated  in  the  computer  by  using  a  DSL  program 
in  which  the  two  parameters  a  and  6  were  set  at  the  above 
indicated  values. 

From  this  program  it  was  found  that  the  time 
response  characteristics  of  the  compensated  singular  system 
were  exactly  the  same  as  these  corresponding  to  the  singular 
root  second  order  model.   It  is  shown  in  figure  (7-5). 

It  must  also  be  noticed  that  the  desired  perform- 
ance of  the  compensated  singular  system  can  be  achieved  by 
an  infinite  number  of  a  and  g  parameter  pair  values,  which 
correspond  to  a  different  selection  in  the  location  of  the 
non-singular  root  on  the  s-plane. 
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Criterion  for  the  selection  of  the  appropriate 
location  for  the  non-singular  root  can  be  based  either  on 
sensitivity  or  and  steady  state  accuracy  considerations  as 
follows : 

a.  Sensitivity  Considerations 

As  the  location  of  the  non-singular  root  moves 
away  from  the  axis  of  reals  on  the  Root  locus  of  figure  (7-4) 
the  value  of  the  Root  locus  variable  parameter  a  increases. 

Considering  next  the  parameter  a-3  plane  and 
the  singular  line  (equation  (7-14)),  which  corresponds  to 
the  operating  singular  root  (c   =  0.40,  oj    ■  2.5),  as  the 
value  of  the  parameter  a  increases  the  operating  point 
(a,B)  moves  away  from  the  origin  on  the  singular  line,  but 
according  to  the  results  obtained  from  the  sensitivity 
study  in  the  previous  section  it  yields  a  less  sensitive 
compensated  singular  system  with  respect  to  variation  on 
the  parameter  a   and  g  values.   These  remarks  can  be  used 
accordingly  in  each  case. 

b.  Stealy  State  Accuracy  Considerations 

When  the  error  constant  of  the  compensated 
singular  system  was  considered  it  was  found  to  be: 

K    -  -L- 
sx    ~ 
la 

When  problem  specifications  concerning  the 

error  constant  K   are  also  established,  then  an  additional 

sx  ' 
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criterion  for  the  selection  of  the  location  of  the  non- 
singular  root  will  be  the  desired  value  of  the  parameter  a 
for  which  the  ratio  g/a  is  such  that  the  optimal  error 
constant  coefficient  can  be  achieved. 

In  general  the  ratio  g/a  can  be  effected  either 
by  selecting  a  different  location  for  the  non-singular  root 
or  by  selecting  a  different  operating  singular  point,  i.e., 
operating  on  a  different  singular  line  or  by  introducing 
additional  gain  adjustments  at  the  forward  path  of  the 
minor  loop  of  figure  (7-2) . 

It  must  be  noticed  also  that  any  other  specifica- 
tions can  be  taken  under  consideration  in  each  design  step 
of  the  previous  illustrated  example. 
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VIII.    EXTENSIONS  OF  THE  SECOND  SINGULAR 
CASE  COMPENSATION  STRUCTURE 

A.    INTRODUCTION 

Among  the  three  singular  cases,  which  were  previously 
analyzed,  the  most  interesting  for  applications  in  design 
of  compensation  of  a  linear  feedback  system  by  the 
dominant  mode  method,  appears  to  be  the  "second  singular 
case"  (See  Section  VI) . 

The  basic  advantages  of  this  case  are  that  the  plant 
type  number  is  preserved,  it  can  be  used  for  any  type  of 
plant  and  for  any  plant  structure. 

On  the  other  hand  the  basic  limitation  of  the  "first 
singular  case"  (See  Section  V)  is  that  in  general  the  type 
number  of  the  plant  is  not  preserved,  and  for  the  "third 
singular  case"  [See  Section  VII)  is  that  the  plant  has  to 
be  at  least  type  two. 

For  the  above  reasons,  the  "second  singular  case"  was 
further  considered  and  some  extensions  of  its  compensation 
structure  were  studied. 
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B.    SECOND  SINGULAR  CASE  COMPENSATION  STRUCTURES 

1.    The  Basic  Structure 

The  basic  compensation  structure  of  the  "second 
singular  case"  is  shown  in  figure  (6-4)  where: 

(s  +  z  )  (S+Z  ). . . (s+Z  ) 

G(s)  =  kFg(s)  =  k_ 1 1 2 (8-1) 

N 
s  (s+Pl)  (s+p2). . . (s+pn) 

represents  the  transfer  function  of  the  plant  to  be  com- 
pensated. 

H(s)  =  a  (8-2) 

represents  a  gain  adjustment,  introduced  as  a  negative 
feedback  from  the  output  to  the  input  of  the  plant,  and 

(s  +  z  )  (s  +  z  )  .  .  .  (s  +  z  ) 
Gc(s)  =  gFg(s)  =  6 I 1 ^ (8-3) 

sN(s+Pi)  (s+p2) ... (s+pn) 

represents  the  transfer  function  of  a  cascade  compensator 
whose  poles  and  zeros  are  identical  with  these  of  the 
plant  and  its  Root  locus  gain  is  the  second  adjustment  of 
the  compensation  structure. 

This  compensation  structure  defined  above,  cor- 
responds to  a  singular  system  which  was  studied  in  Section 
VI  as  the  "second  singular  case". 
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Its  transfer  function  has  the  general  form: 


m 

2 


trcs)  -£&!  - 


ek  I  |  Cs+z.) 
'All -i 


B     r,  r    -\  i.   n  n       m         m 

rcs)       2Nrn/-_.  *?  ,  _ni 


ncs+PK^aksNr~ics+pj  rics+zo^krics+z,)' 

k=i    K  k=i    K  j=i    J    j=i    J 


(8-4) 


The  error  constant  K   .,  was  evaluated  and  it 

sx-B 


was  found  to  be 

m 


ksx-b  =  it— - t8"5) 

n 

n 


k=i 


The  characteristic  equation  also  of  this  basic 
singular  compensation  structure  was  found  to  be  in  the 
general  form: 


1+akF  (s)  +  gkF  2(s)  =  0  (8-6) 

From  this  basic  structure  six  modified  singular 
compensation  structures  were  derived.   They  were  analyzed 
as  follows: 
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2.   Modified  Structures 

a.   First  Modified  Structure 

The  first  modified  singular  compensation 
structure  is  shown  in  figure  (8-1) .   It  was  derived  from 
the  previous  basic  structure  and  has  been  analyzed  in 
Section  VI. 

In  figure  (8-1)  GCs)  denotes  the  transfer 
function  of  the  plant  to  be  compensated  and  G   (s)  the 
transfer  function  of  the  singular  cascade  compensator, 
which  according  to  equation  (6-18)  its  transfer  function 
is: 


m 

e I I  Cs+Zj) 

GsJs)    -.— Q in C8'7) 


N 
s 


r~i(s+p  )+0  n  cs+z.) 

k=l       K       j=l        J 


k= 

where  s  ,  p-.  and  z-,  represents  the  poles  at  the  origin, 
the  other  poles  and  the  zeros  of  the  plant  respectively. 
The  transfer  function  of  this  compensation 
structure  was  found  to  be: 

(8-8) 
m 

ccs)  ekrb^)2 

T  Cs)  =  —  = !zl__ _ _ 


US)  2NPH       2 

s 


n  n       m         m 


rics+pj2+asNrics+p  j  rics+z.o+ek  r~ics+z.)' 

k=i         K  k=i  j=i         J  j=i         J 
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The  error  constant  K     was  also  found  (See 

sx-1  v 


equation  (6-85))  to  be: 


K      =  _1_  Jal  '   -   3  K 


SX-1     a 


m 
3             j'-'i    J 

a 

np* 

k=i 

a    x 


(8-9) 


where  K   denotes  the  error  constant  of  the  plant,  before 

.A. 

it  was  compensated  by  the  singular  cascade  compensator. 

The  characteristic  equation  of  this  compensa- 
tion structure  will  be  in  the  general  form: 

l  +  aFg(s)+3kFg2(s)=  0  (8-10) 

The  basic  advantages  of  this  modified  singular 
compensation  structure  with  respect  to  the  previous  basic 
one,  are  the  following  two: 

(1)  Simpler  Structure.   The  compensation 
structure  is  simpler  because  only  a  single  cascade  filter 
is  required.   Historically  it  is  realized  that  series 
compensation  has  played  a  prominent  role  in  system  design. 

(2)  Greater  Error  Constant.  Comparison  of 
equations  (8-5)  and  (8-9),  reveals  that  the  error  constant 
of  this  structure  is  greater  than  the  error  constant  of 
the  previous  basic  one  by  a  factor  equal  to  the  Root  locus 
gain  k  of  the  plant. 
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b.   Second  Modified  Structure 

The  second  modified  singular  compensation 
structure  is  shown  in  figure  (8-2).   It  was  derived  from 
the  previous  one  by  changing  the  location  of  the  singular 
compensator  from  the  forward  path  at  the  feedback  path. 
Of  course  this  modified  singular  structure  is  no  longer 
a  unity  feedback  control  system. 

For  this  compensation  structure  the  transfer 
function  was  evaluated  and  it  was  found  to  be: 

n       m  m 


ksNr~i(s+pj  rics+z,)+ak  n  cs+z,y 

m  C(s2  k=i         K     j  =  i         J  j  =  i  J 

n  n  m 

ri(s+pj2+asNri(s+p  jri(s+z,)+ek  n  cs+z.)' 

k=i        K  k=i        K  j=i        J  j=i  1 


T2  i5)=  RTiJ  " n n 'In  ~ in 

2Nr~l.  2       N[ — l^         J      I  _     I — |  2 


(8-11) 

The  characteristic  equation  also  of  this 
compensation  singular  structure  was  found  to  be: 

1+G(s)    g'(s)    =    0 

or 

l+aFg(s)+3k  Fg2(s)    =    0  (8-12) 

In  order  for  steady  state  accuracy  to  be  in- 
vestigated the  following  considerations  were  made: 
The  error  of  the  system  is: 

E=R-C=(1-T2)R  (8-13) 
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Introducing  equation  (8-11)  into  equation  (8-13) 


yields: 


2N 


E(s)=- 


n  n  m  m 

n(s+pj2+(a-k)sNncs+pjncs+z-i)--(ek-ak)  n  cs+z.y 

k=i        K k=i        K  j  =  i        J j  =  i         J 


k=_ 
n 


n 


m 


m 


R(s) 


s  n(s+Pj2+asN  n  (Sip j  n  (s+z.)+ek  n  cs+z_o' 

k=i       j=i     J     j=i     J 


k=i 


rK' 


(8-14) 

When  a  type  one  plant  was  considered  and  a  ramp 

2 
input,  i.e.,  N=l  and  R(s)=B/s  ,  and  when  the  final  value 

theorem  was  applied  to  equation  (8-14) ,  assuming  a  stable 

system,  the  steady  state  error  was  evaluated  as: 

E(t)=Ess(t)  =  - 

t+oc 

For  a  physical  interpretation  of  the  above  re- 
sult the  following  considerations  were  also  made  on  the 
compensation  structure  which  is  shown  in  figure  (8-2) . 

For  a  type  one  plant,  and  a  step  input,  at 

steady  state  E  has  to  be  zero,  i.e.,  E   =  0.   Because 

E  =  R-C  ,  it  implies  that  R  =  C   .   Since  C   is  a 
'      r  ss    ss  ss 

constant  and  the  system  is  not  a  unity  feedback  system 

Ceo  t   C„.   Therefore  C   f   R. 
ss    ss  ss 

When  a  type  one  plant  was  again  considered 
with  a  ramp  input,  at  steady  state  E  has  to  be  constant, 
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i.e.,  E   =  Constant.   Because  E  =  R-C   this  implies  that 

R  =  C   .   Since  C   is  a  ramp  and  the  system  is  not  a 
ss  ss         ^  J 

•  •  *  •         • 

unity  feedback  system,  C„  f   C   .   Therefore  C   =  R, 

J  J  '   ss  '   ss  SS    ' 

which  implies  that  an  infinite  error  E   must  be  expected. 
This  result  was  also  derived  previously  by  analytical 
study  of  equation  (8-14)  . 

Based  on  these  remarks  and  results,  the  con- 
clusion is  that  such  a  system  does  not  reproduce  the  input 
signal  at  the  output  or  in  other  words  at  steady  state  the 
output  does  not  follow  the  input,  which  may  be  an  advan- 
tage or  a  disadvantage  depending  on  the  specific  applica- 
tion, and  which  can  be  eliminated  as  will  be  seen  by  the 
following  modified  singular  compensation  structure. 

When  a  type  one  plant  was  considered  again  and 
a  step  input  applied,  the  steady  state  error  was  found  to 
be: 

Ess=  (1-  |  )  A 

where  A  is  the  value  of  the  step  input,  i.e.,  at  steady 
state  the  output  differs  from  the  input  by  the  value  of 

Ess' 

In  order  for  all  the  above  mentioned  problems 
to  be  eliminated  when  this  is  desired,  the  following 
considerations  were  made  on  this  modified  singular 
structure. 
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When  R  was  considered  as  the  input  signal  which 
the  output  of  the  system  was  desired  to  reproduce,  i.e., 
to  follow,  then  the  error  will  be  given  by  equation: 

E  =  R-C  (8-15) 

By  considering  also  an  input  which  is  not  followed 
by  the  output  to  be  R  ,  then  according  to  the  notations 
used: 


C  =  T2R 


(8-16) 


Introducing  equation  (8-16)  into  equation  (8-15)  yields 


E  =  R-T2R 


(8-17) 

By  assuming  that  the  actual  input  to  the  system 
R  is  correlated  to  the  desired  input  R  by  a  function  Y, 
i.  e: 

r'=  YR  (8-18) 

then  equation  (8-17)  becomes: 

E  =  R-T2  YR  =  (1-T2Y)  R  (8-19) 

Introducing  equation  (8-11)  into  equation  (8-19)  yields: 


E(s)   = 


n  n  m  m 

NTl(s+pJ2+(a-kY)sNri(sX)  riCs+zO  +  Cek-akY)  1      Ics+z- 
_c=l kj*l j  =  l  J j  =  l  J 


2N n  n  m  m 

>  ri(s+pj2+asNri(s+pj  n  cs+z^+gk  n  cs+z.) 

k-l         K  k=i         K     j  =  i  J  3  =  1  J 


-R(s) 


(8-20) 
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In  order  for  the  steady  state  error  to  be 
finite  for  a  ramp  input  to  a  type  one  plant  or  zero  for  a 
step  input,  Y  has  to  be  defined  as: 


Y  -=   i 


C8-21) 


Introducing  equation  (8-21)  into  equation  (8-18)  yields: 


R  =  -  R 

a 


(8-22) 


After  the  above  result  was  obtained  a  third 
modified  structure  was  derived  and  is  analyzed  as  follows: 

c.   Third  Modified  Structure 

The  third  modified  singular  compensation 
structure  is  shown  in  figure  (8-3),  where  R  is  the  input  to 
the  system  which  is  to  be  followed  at  steady  state  by  the 
output  of  the  system. 

The  transfer  function  of  this  modified  struc- 
ture was  found  to  be: 


n 


ks 


N 


m 


m 


MCs+pJ  |   lcs+z.)-Hok  j   I  (s+Z.) 

T  ro_  c(s)  =  g   k*i    j=i   J   j=i   J 


n  n        m  m 

;2Nrks+pj2+asN  n  (s+pj  n  cs+z^+ekii  cs+zj* 

k=i        K  k=i  K    j=i  J  j=i  J 


(8-23) 
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The  characteristic  equation  was  also  found  to 
be  in  the  general  form: 

1+a  Fg(s)+Bk  Fg2(s)  =  0  (8-2 

By  the  same  procedure  as  in  the  case  of  the 


first  modified  structure,  the  error  constant  K    ,  for 

sx-3 

this  case  was  found  to  be: 

m 


nzi 


K   -«  -yiL.  k  -J£i =  -4^ K  (8-25) 

sx-3   2    „-,_     n       2   ,   x  K  J 


a    -3k 


k=i 


a  -gk 


where  K   denotes  the  error  constant  of  the  plant  before  it 
was  compensated.   For  the  examples  which  were  worked  out 

2 

in  Section  VI,  the  ratio  ag/{a  -3k}  turns  to  be  greater 
than  unity.   Therefore,  when  such  is  the  case  the  error 
constant  of  the  uncompensated  plant  will  be  improved  when 
the  plant  is  compensated. 

When  this  compensation  structure  was  compared 
with  the  first  modified  structure  it  was  found  that  they 
have  the  same  characteristic  equation  but  different  error 
constant  and  transfer  function.   Therefore  the  stable 
region  on  the  a-3  plane  will  be  the  same  in  both  cases. 
Also  for  the  same  operating  point  Ca>3)>  although  they 
will  have  the  same  characteristic  roots,  the  residues  of 
these  roots  will  be  different  in  each  case  due  to  the  dif- 
ferent numerator  of  their  transfer  functions.   Since  the 
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dominancy  of  a  characteristic  root  is  a  function  of  its 
residue,  in  addition  to  the  distance  of  the  location  of 
this  root  on  the  s-plane  from  the  imaginary  axis,  a  dif- 
ferent degree  of  dominancy  for  the  operating  singular 
root  must  be  expected  in  general  in  each  case. 

For  a  specific  plant  and  a  specific  operating 
point  on  the  a-g  plane  a  further  investigation  of  the 
dominancy  of  the  operating  singular  root  can  be  performed 
either  by  a  computer  simulation  of  each  compensation  • 
structure  and  comparing  their  time  response  characteristics 
with  those  corresponding  to  the  singular  root  second  order 
model  or  by  analytical  evaluation  of  the  residues  cor- 
responding to  the  characteristic  roots  of  each  case. 

d.   Fourth  Modified  Structure 

The  fourth  modified  singular  structure  is 
shown  in  figure  C8-4) .   It  was  derived  from  the  first  one 
by  properly  defining  the  feedback  compensator  in  order  for 
these  two  compensation  structures  to  have  the  same  open 
loop  transfer  function,  i.e.: 


G(s) 

GCs)  G   fs)  =  3 (8-26) 

1+GCs)  G   (s) 
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from. which  the  transfer  function  of  the  filter  GscCs)  was 
evaluated  to  be: 


n         n         m         m 

N  " 


GscO)  = 


ncs+pkMsN  n(s+Pi>cc  i  ics+z  )-b  rics^i)} 

k=i    K    k=i    K    3=1    J    j=l    J 


m      2 
3k  r~|(s+z.j) 

j  =  1 

(8-27) 


Because  GCs)  defined  by  equation  (8-1)  denotes  the 
transfer  function  of  a  physical  control  system,  it  implies 
that  N+n>m,  i.e.,  the  above  defined  filter  has  more  zeros 
than  poles.   In  practice  a  structure  of  this  type  is 
realized  by  introducing  additional  poles  as  needed,  but 
placing  them  well  outside  the  bandwidth  of  the  system. 
Whether  this  can  be  done  without  affecting  the  singularity 
characteristics  has  not  been  investigated.   Based  on  the 
argument  requirement  for  singular  conditions  (See  equation 
(3-12  and  (10-11)),  it  is  believed  that  if  these  additional 
poles  are  placed  on  the  s-plane  in  such  locations  that 
their  net  contribution  to  the  argument  of  G   (s)  is 
negligibly  small,  then  it  can  be  assumed  that  the  singu- 
larity characteristics  of  the  system  remain  unaffected. 
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The  transfer  function  of  this  structure  was  found 


to  be: 


T  CO  =   CCs)  -   G(s) 


1+G(s)+G(s)  Gsc(s) 


(8-28) 


or  after  some  manipulations 


m 


T„  (s)   - 


C(s) 


ek'; 


2-1 


Cs+Zj.) 


n  n  m  m 

n(s+Pk)2+asN  I      |(s+pk)  ri(s+z,)+  3k[      Ics+z)' 
k-i        K  k=i        K    j-i        J         j-i        J 


2N 


(8-29) 

The  characteristic  equation  of  this  compensation 
structure  was  found  to  be  in  the  general  form: 


l+aFg(s)+ek  FgZ (s)  =  0 


(8-30) 


By  the  same  procedure  also,  as  in  the  case  of  the 

first  modified  structure,  the  error  constant  K    .    for 

sx-4 


this  case  was  found  to  be: 


sx-4   a   x 


(8-31) 


Comparison  of  the  results  which  were  obtained  for 
the  first  modified  compensation  structure  with  these 
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obtained  above,  reveals  that  the  performance  of  the  system 
in  both  cases  is  the  same. 

e.   Fifth  Modified  Structure 

The  fifth  modified  singular  compensation 
structure  is  shown  in  figure  (8-5).   It  was  derived  from 
the  second  one  shown  in  figure  (8-2)  by  splitting  the 
feedback  path  into  two  parallel  paths  as  shown,  such  that 
their  sum  is  the  compensator  transfer  function: 

Gsc(s)  =    l+GgCCs.)  (8-32) 

from  which  the  transfer  function  of  the  minor  loop  feed- 
back filter  was  found  to  be; 

m         n 


.-a)rL+z.)-sNn 


G-_(s)  =  JO± J l£i £ (8-33) 

sc  n         m  v    J 


,..      m     (3-a)tj='(s+Zj)-s  y  (S+Pk) 
n         m 

SN  rics+pk)+«ri  (s+z. 

k=i   K   j=i    3 


The  transfer  function  of  this  compensation 
structure  is: 


T  (S)  =   CCS)  =       GCs) 
S       RCs)   1+G(s)+G(s)  G  sc(s) 


ii 
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or  after  some  manipulations  it  becomes 


T  (s)  = 


-  C(s)  fe 


^TiT 


n       m  m 

ksN  n  (s+pk)  in  cs+zi ) +<x  k  r  i  cs+z ) : 

)  =  k=i j=1        J j  =  i  J 


n  n  m  m 

|     j(s+pk)  +asN|     !(s+pv)  |     |cs+z.)+0k|   li 
k=i        K  k=i        K    j-i        J  j=i 


(8-34) 


The  characteristic  equation  of  this  compensation 
structure  was  found  to  be  in  the  general  form: 

l+aFgCs)  +  BkFg2(s)  =  0  (8-35) 

Also  when  steady  state  accuracy  considerations 
were  made  for  this  compensation  structure  the  results  were 
the  same  as  these  for  the  second  modified  structure.   Com- 
parison of  the  results  corresponding  to  these  two  cases 
reveals  that  they  are  equivalent  under  performance  con- 
siderations . 

f.   Sixth  Modified  Structure 

This  is  shown  in  figure  (8-6),  and  it  was  de- 
rived from  the  third  modified  structure  with  the  same 
reasoning  by  which  the  fifth  modified  structure  was  de- 
rived from  the  second  one. 
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It  was  also  found  that  the  third  and  sixth 
modified  structures  have  the  same  transfer  function, 
characteristic  equation,  error  constant  and  thus  under 
system  performance  considerations  they  are  equivalent. 

3.   Remarks 

Although  the  "second  singular  case"  can  be  used 
for  compensation  of  any  control  system,  i.e.,  it  may  have 
general  application,  its  "basic  structure"  has  certain 
disadvantages,  which  have  been  stated  in  the  previous 
analysis  and  which  may  not  be  acceptable  in  some  cases. 

When  the  characteristic  equation  of  the  basic 
structure  is  considered,  it  will  have  the  general  form: 

1+a.kF  Cs.)+6kF  2(s)  =  0 

On  the  other  hand  the  six  modified  structures, 
which  were  analyzed  above,  have  certain  advantages  over 
the  "basic  structure"  and  a  common  characteristic,  i.e., 
they  have  the  same  characteristic  equation,  which  has  the 
general  form: 

l+aFgCs)+BkFg2Cs)    =    0 

Also  their  performance  characteristics  can  be  mathematical- 
ly correlated. 

This  fact  aids  the  designer  in  making  his  decision 
on  the  problem  of  the  structure  without  additional  design 
effort  and  protects  from  arbitrary  choices. 
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In  the  following  example  all  these  modified 
singular  compensation  structures  were  used  to  compensate 
the  same  plant  in  order  that  the  derived  analytical 
results  could  be  compared.   The  compensated  system  per- 
formance was  also  studied  in  each  of  these  cases. 

4.   Example  (8-1) 

In  this  example  the  plant  of  example  (6-7)  was 
considered  for  which  the  transfer  function  G(s)  of  the 
system  before  any  singular  filter  has  been  introduced  is: 


20(s+2) 

G(s)  =  (8-36) 

s(s+l)  (s+5)  (s+20) 


Also  as  operating  point  from  the  dominant  root 
region  of  the  a-3  plane  corresponding  to  the  singular  root 
(?s=0.60,  u   ■1.23),  the  point  (a=920,  3=2500)  was  ar- 
bitrarily selected.   Then  each  of  the  six  modified  sin- 
gular compensation  structures  was  considered  separately 
and  the  performance  of  each  for  a  step  input,  R=l  Volt  and 
a  ramp  input  R=0.1  rad/sec,  was  studied  in  the  computer  as 
follows: 

a.   First  Modified  Structure 

This  is  shown  in  figure  (8-1) .   According 
to  equations  (8-7)  and  (8-36)  the  transfer  function  of 
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the  singular  cascade  compensator  G   (s)  was  found 
to  be : 


G_Cs)  = 


2500  (s+2) 


sc 


s(s+l)  (s+5)  (s+20)+920  (s+2) 


By  simulation  of  this  system  in  the  computer 
with  the  use  of  a  DSL  program  the  time  response  character- 
istics were  obtained  for  a  step  and  a  ramp  input.   They 
are  shown  in  figure  (8-1A)  and  (8-1B)  respectively. 

The  characteristics  of  the  transient  response 
corresponding  to  the  singular  root  (£  =0.60,  uj   =1.23) 
were  evaluated  and  they  were  found  to  be: 


t   =  3.19  sec 
P 


t  =  5.41  sec 
s 

Mp  =  9.48% 

oj  =  0.99  rad/sec 


Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot 

Transient  oscillatory  frequency: 

When  the  step  input  was  considered  then  it 
was  found  from  the  computer  print  output  that: 
Steady  state  output  value 
Time  to  peak  of  first  overshoot 
Settling  time 
Percent  overshoot 
Transient  oscillatory  frequency: 


Css  =  R  =  1  Volt 


*P 


=  3.2  sec 


ts  =  5.4  sec 


Mp  =  111 

wt  =  0.98  rad/sec 
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When  the  ramp  input  was  considered  then  it 
was  found  also  from  the  computer  print  output  that: 
Steady  state  output  value 
Steady  state  output  rate 
Settling  time 


C    =  R-0.002 
ss 

C    =  R  =  0.1 
ss 


ts   =  5.4  sec 

These  results  can  also  be  derived  approximately  from  the 
corresponding  figures  of  the  time  response  characteristics, 
figure  C8-1A.B) . 

b.   Second  Modified  Structure 

When  the  singular  compensator  G   Cs)  was 
located  at  the  feedback  path  as  shown  in  figure  (8-2),  the 
transient  response  characteristics  which  were  obtained  are 
shown  in  figure  (8-2A)  and  (8-2B)  for  a  step  and  a  ramp 
input  respectively. 

The  exact  transient  and  other  characteristic 
values  which  were  read  from  the  computer  print  output  have 
as  follows: 

(1)   Step  Input. 


Steady  state  output  value 

Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot  from  C 

ss 

Transient  oscillatory  frequency 


ss 


s 

Mp 

a). 


0.368    =   £  R 

p 

3 .  2    sec 

5 . 4    sec 

10.2% 

0.92    rad/sec 
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(2)   Ramp  Input . 

Steady  state  output  value      :   C   followed  a  ramp 
different  than  the  input  ramp. 

•  • 

Steady  state  output  rate       :   C    =  0.036  »  £L  R 

S  S  p 

Settling  time  :   t   =  5.4  sec 

c.   Third  Modified  Structure 

When  a  gain  adjustment  equals  with  2.717,  which 
corresponds  to  the  ratio  3/a  value,  i.e: 

8.  rsm  m   2>717 

a     920 

was  introduced  at  the  indicated  in  figure  C8-3)  location 
the  third  modified  singular  structure  was  obtained. 

The  transient  response  characteristics  which 
were  obtained  in  this  case  for  the  same  step  and  ramp 
input  as  in  the  previous  cases,  are  shown  in  figure  (8-3A) 
and  (8-3B)  respectively. 

The  exact  transient  and  other  characteristic 
values  which  were  read  from  the  computer  print  output  have 
as  follows: 
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(1)   Step  Input 


Steady  state  output  value 

Time  to  peak  of  first  overshoot 

Settling  time 

Percent  overshoot 

Transient  oscillatory  frequency 


ss 


ts 

Mp 

CO  j. 


1    =    R 

3. 2    sec 

5 . 4    sec 

10.21 

0.92    rad/sec 


(2)      Ramp    Input. 


Steady  state  output  value 
Steady  state  output  rate 
Settling  time 


C    =  R-0.002 
ss 

C    =  0.1  =  R 
ss 

tc   =5.4  sec 


d.   Fourth  Modified  Structure 

The  fourth  modified  structure  which  is  shown  in 
figure  [8-4)  requires  a  feedback  filter  G   (s) .   After  the 
transfer  function  G   (s)  was  evaluated  from  equations  (8-27) 
and  (8-36)  the  system  was  simulated  in  the  computer.   The 
transient  response  characteristics  for  a  step  R=l  Volt 
and  a  ramp  R=0.1  rad/sec  input,  which  were  obtained  in 
this  case  were  the  same  with  these  corresponding  to  the 
first  modified  structure.   They  are  shown  in  figure  (8-4A) 
and  (8-4B)  respectively. 
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e.  Fifth  Modified  Structure 

This  modified  structure  involves  a  different 
feedback  singular  filter  G  _(s)  whose  transfer  function 
was  evaluated  from  equations  (8-33)  and  (8-36).   The 
transient  response  characteristics  for  a  step  R=l  Volt 
and  a  ramp  R=0.1  rad/sec  input,  which  were  obtained 
when  this  structure  was  simulated  in  the  computer  are 
shown  in  figure  (8-5A)  and  (8-5B)  respectively,  which  are 
the  same  with  these  corresponding  to  the  second  modified 
structure. 

f.  Sixth  Modified  Structure 

This  is  shown  in  figure  (8-6)  and  was  implemented 
from  the  previous  one  by  introducing  a  gain  adjustment  as 
it  is  indicated  in  figure. 

Simulation  of  this  system  in  the  computer  for 
a  step  R=l  Volt  and  a  ramp  R=0.1  rad/sec  input,  gave  the 
same  transient  and  steady  state  characteristics  as  in  the 
case  corresponding  to  the  third  modified  structure.   They 
are  also  shown  in  figures  (8-6A)  and  (8-6B)  respectively 
for  each  input  which  was  considered. 
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5.   Comment 

The  above  presented  results  which  were  obtained 
from  the  computer  not  only  confirm  the  corresponding 
analytical  results  but  also  make  clear  the  flexibility 
which  the  designer  has  when  compensation  is  performed  by 
using  singular  line  theory. 


182 


IX.    DESIGN  OF  COMPENSATION  BY  USING  A 
SINGULAR  CASCADE  COMPENSATOR 

A.    INTRODUCTION 

Consider  the  plant  which  is  shown  in  figure  (4-1), 
and  that  this  plant  is  going  to  be  compensated  by  a  sin- 
gular cascade  compensator.   The  open  loop  transfer  function 
of  this  plant  is: 


(s  +  z, )  (s  +  z  )  .  .  .  (s  +  z  ) 

G(s)  =  kFgCs)  =  k  \      1   11 . K (9-1) 

s  (s+p^  (s+p2)  .  ..  (s+pn) 

It  was  assumed  that  the  above  transfer  function  is 
known  and  that  the  problem's  specifications  have  been 
interpreted  in  terms  of  the  desired  dominant  characteristic 
root  (?,  w  )  for  the  compensated  singular  system.   It  was 
also  assumed  that  tolerance  on  the  polar  coordinates  ?  and 
a>   values  for  this  root  have  been  established  from  the 
problem's  specifications. 

The  compensated  singular  system  is  shown  in  figure 
(6-6)  where   G   (s)  denotes  the  transfer  function  of  the 
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required  singular  cascade  compensator,  whose  standard  form 
is : 

(9-2) 


B(s  +  zJ  (s  +  z  )  .  .  .  Cs  +  zl 

g    (s)  =  — i— il — : — a: 

sc  v  J 


s    Cs+p^    Cs+p2)  .  .  .  Cs+pn)+a(s  +  zi)    (s  +  z2.  .  .  (s  +  zm) 


as  was  derived  in  Section  VI  C. 

Therefore  the  design  problem  in  general  will  be  limited 
to  the  problem  of  determining  the  required  parameter  pair 
values.   For  different  types  of  applications  this  problem 
can  be  stated  in  terms  of  the  parameter  plane  as:  What  must 
be  the  operating  point,  or  the  operating  line  segment,  or 
the  operating  region  of  the  parameter  plane  in  order  that 
the  performance  of  the  compensated  singular  system  be  as 
desired? 

The  required  design  procedure  in  all  these  cases  has 
been  developed  in  Section  VI.   In  this  section  the  design 
process  in  the  case  where  an  operating  point  is  desired  to 
be  defined  on  the  parameter  plane  and  from  this  the  pre- 
cise form  of  the  transfer  function  (9-2),  was  formulated 
in  a  sequence  of  basic  design  steps  which  must  be  followed 
in  a  case. 

B.    DESIGN  PROCESS 

The  required  design  process  in  general  when  an  ac- 
ceptable operating  point  of  the  parameter  plane  is  desired 
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to  be  found,  is  relatively  simple  and  can  usually  be 
completed  with  satisfactory  accuracy  for  most  applications 
with  graphical  and  mathematical  techniques,  by  the  fol- 
lowing four  basic  design  steps, 

1.  First  Step 

The  first  step  consists  of  drawing  the  Root  locus 
of  equation: 

kFg(s)  =  -1  (9-3) 

2 .  Second  Step 

The  second  step  consists  of  insertion  of  an  ini- 
tial cascade  compensator,  when  the  shape  of  the  above  Root 
locus  is  such  that  either  an  acceptable  operating  singular 
point  cannot  be  located  on  this  Root  locus  or  although  it 
can  be  located  it  is  not  acceptable  under  dominancy  consid- 
erations (See  VI,  F,2). 

When  an  initial  cascade  compensator  is  required 
then  this  requirement  will  actually  affect  only  the 
transfer  function  of  the  singular  cascade  compensator  and 
neither  the  complexity  of  the  rest  of  the  design  process 
nor  the  compensation  structure  (See  example  (6-6)  and 
figure  (6-28)). 

In  the  simple  case  where  such  an  initial  compen- 
sator will  not  be  required  then  this  step  will  consist  of 
determining  an  acceptable  operating  singular  point  on  the 
Root  locus  of  equation  (9-3), 
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3.  Third  Step 

The  third  step  consists  of  determining  the 
equation  of  the  singular  line  corresponding  to  the 
selected  operating  singular  point  (See  VI  D3a  equation 
(6-29)). 

4.  Fourth  Step 

This  step  which  is  the  last  one  consists  of  de- 
termining the  required  operating  point  of  the  parameter 
plane,  i.e.,  the  required  pair  of  a  and  6  parameter  values 
(See  examples  (6-6),  C6-7),  (6-8)). 
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EXTENSION  OF  THE  SINGULAR  LINE  THEORY 


TO  PARAMETER  SPACE 


A.    INTRODUCTION 

The  parameter  space  is  an  extension  of  the  parameter 
plane  when  there  are  three  or  more  parameters. 
Considering  the  characteristic  equation: 

n 

F(s)  =  H  a,  sk  =  0  (10-1) 

k=0   K 

where  n  is  a  real  integer  and  a,  is  a  real  valued  coef- 
ficient which  may  be  a  function  of  three  or  more  parameters, 
i.e: 

ak  -  gO,e,Y,<5,  ...  nth)  (10-2) 

where  n   denotes  the  n   parameter. 

The  point  s  can  be  defined  in  polar  coordinates  as: 


s  =  "^n  t  J  %  /T7?  C10-3) 


B.    SINGULAR  CONDITIONS 


If  for  a  given  pair  c  ,  u   values,  where  | t  | <1  and 
to  >0,    equation  (10-1)  has  an  infinite  number  of  real  valued 
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sets  of  a,3,Y,<5>  •••  n   parameter  values  which  satisfy 

it,  then  for  this  equation  singular  conditions  exist.   In 

this  case  the  point  (?  ,  to   )  is  a  singular  point.   This 
r       s'   ns  r 

equation  also  for  which  singular  conditions  exist  is 
called  singular. 

C.    SINGULAR  CONDITIONS  ANALYSIS  IN  THE  LINEAR  CASE 

A  simple  case  which  was  analyzed  is  the  linear  case, 
in  which  the  variable  parameters  appear  in  the  coefficients 
of  the  characteristic  polynomial  in  the  general  form: 

ak  *  bk  a+ck  6+dk  Y+ek  6+  *•   +  ••• 

+{Coefficient  of  the  n   parameters}, {n   parameter}+ 
+{Constant  term}, 

where  b-.,    c,,  .,.  are  real  constants. 

In  this  case  equation  (10-1)  can  be  rearranged  in  the 
form: 

f1Cs)a+f2Cs)    g+f3(s)    y+    ...+fn(s){n        parameter}* 

+fn+1(s)    =    0  (10-4) 

where  f , (s) ,  f~(s),  ...  are  polynomials  of  s.   By  analogy 
as  in  the  case  of  the  parameter  plane  (See  Section  IIB2) 
equation  C1Q-4)  can  be  written  as: 

r^6!  «+r2ej92  g+r  ej93  y+...+Tne?9n{nt    parameter}+rn+1e:,9n+1=0 
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or 


i-j  (cose  x+jsine  1)a  +  r2  Ccose2+j  sin02)  g+.  .  .  + 


+r    CcosO„+jsine„)  {n      parameter}+r  ^fcosQ       +jsine   in)  =  0 
n  ^  n  n  *  n+i  n+ 1  ^         n+ 1 


Equating    then   real    and   imaginary  parts    to    zero, 
yields : 


r.sine     a+r    sine     B+...+  r  sin  e  {  n     }  +  r         sine  =  0 

i         1         2  2  n  n  n+i  n+i 


and 


1~Vi 

r,cosen   a+r  cose„     g+...+  r  cose  {  n     }  +r-     cose   L     =  0 
1         l         2         2  n        n  n+i         n+i 


According    to    [1]    the   above    system  has   an    infinite 
number   of   solutions    if   and   only   if: 


r.sing  r    sineo      r.sine,  r   sine  r   ^    sine    ^ 

i  l    _      2  2_      3        v3_  n  n   .       n+i  n+i 


r1cosei        r2cose2      r3cose3  rncos9n        rn+icosen+i 


or 

tan9:  =  tane   =  tane„  = .  .  .  *  tane   =  tane  ^  (1Q-5) 

i       2        3  n       n+ 1         *-    J 


Equation  (10-5)  implies  that  for  the  characteristic 
equation  (10-1)  or  (10-4)  singular  conditions  exist  if 
and  only  if: 
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/flCs)      -/fn+1C5)    t    Kx    H 

AiS)     "  A+lCsJ    +  K2   n 

(10-6) 

,/VS)      »  Al*l.CsJ    +   K3   n 


Assuming  that  f   ,(s)  $    0  and  dividing  equation 
C10-4)  by  fn+1Cs)  yields: 

fl&  +   f2ls)         ft+   Vs)  +        +   fn^  th+1    fl 

ct  +   3+   Y+-'*+   n      +1  =  0 

fn+l<^  fn+l^  Ws>  Ws> 

(10-7) 


or 


F1(s)    a+   F2(s)    B+    F3(s)    Y+-..+    Fn(s)    nth+l    =    0 


(10-8) 


When  the  characteristic  equation  is  considered  in  the 
form  of  equation  (10-8)  then  the  necessary  and  sufficient 
conditions  for  singular  conditions  given  by  equation 
(10-6)  become: 
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/F2(s)  =  K2n  (10-9) 

/*S<^  =  K3n 


where  any  of  the  K,  ,  K2,...K  is  real  integers 


D.    SINGULAR  THEORY  AND  CONTROL  SYSTEMS 

From  equations  (10-8)  and  (10-9)  and  by  the  same  pro- 
cedure which  was  used  in  Section  III  for  the  case  of  the 
two  parameters,  it  was  found  that  in  this  case  of  the 
parameter  space  where  there  are  n  parameters  in  the 
characteristic  equation  of  a  control  system  the  general 
form  of  a  singular  characteristic  equation  is: 


m,  m9  m,  ,  m 

aCjF   1Cs)+    gC2F    Z(s)+Y    C3F    5(s)+...+    (ntn}CnFn      (s)+Cn+1=0 


whe 


(10-10) 


re  C, ,  C0,  C,,...,  C  ,  C  ^,  are  real  constants  and  m,  , 
1*      2'   3*    '   n'   n+1  1' 


m~,  m,,...,  m  are  n  different  real  integers.   (Otherwise 
when  there  are  two  or  more  of  these  identical  then  the 
number  of  the  actual  parameters  can  be  reduced  by  combining 
terms  of  equal  powers  of  F(s)). 
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For  equation  (10-10)  the  necessary  and  sufficient 
conditions,  for  singular  conditions  are  given  by  equation: 

/f(s)  =  Kn  CIO -11) 

where  K  is  a  real  integer. 

Equation  (10-11)  implies  that  all  the  points  of  the 
Root  locus  of  equation: 

KF(s)  =  -1  (10-12) 

which  do  not  lie  on  the  real  axis  are  singular  points  for 
the  characteristic  equation  (10-10),  where  K  in  this  case 
is  considered  to  be  the  Root  locus  variable  parameter  which 
takes  values  from  zero  to  plus  infinity  and  from  zero  to 
minus  infinity. 

Finally  any  control  system  with  n  variable  parameters 
appearing  in  its  characteristic  equation  linearly  depen- 
dent, is  singular,  if  and  only  if,  its  characteristic 
equation  can  be  written  in  the  form  of  equation  (10-10) 
which  in  the  simplest  case  takes  the  form: 

1+  aF(s)  +  BF2(s)  +  yF3(s)  +...nth  Fn(s)  =  0   (10-13) 

1.   Equation  of  a  Singular  n  Dimensional  Surface 

When  an  n  dimensional  parameter  space  was  con- 
sidered then  the  general  form  of  the  corresponding  singular 
characteristic  equation  of  a  control  system  is  given  by 
equation  (10-10).   Then  according  to  the  previous  analysis 
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when  the  value  of  a  singular  point  s   is  substituted  in 
equation  (10-10)  it  becomes: 


aC1  +  SC2   +  yC3  +...+  nL"  Cn  +  Cn+1  =  0       (10-14) 


X        .*  «•        ^ 


where  C, ,  C0,  C,,...,  C  ,  C  _,  are  real  constants. 
1*   2    3'    '   n'   n+1 

Equation  (10-14)  is  the  equation  of  the  singular  n 
dimensional  surface,  corresponding  to  the  singular  point 
s  ,  which  was  considered. 

In  order  to  determine  a  singular  point  of  this 
singular  surface  the  values  of  n-1  parameters  must  be 
defined  and  then  the  value  of  the  parameter  left  can  be 
evaluated  from  equation  (10-14).   Thus  one  among  an  infinite 
number  of  real  valued  sets  of  parameter  values  can  be  de- 
fined, which  in  other  words  defines  the  coordinates  of  a 
singular  point  in  the  parameter  space. 

E .    COMMENT 

From  the  previous  analysis  it  can  be  seen  that  any 
non-linear  case,  either  in  the  parameter  plane  on  in  the 
parameter  space  can  be  studied  by  defining  any  different 
parameter's  product  as  a  new  variable  parameter,  e.g.,  the 
non-linear  case  in  the  parameter  plane  can  be  considered 
as  a  linear  case  in  the  parameter  space  (three  dimensions 
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for  this  case)  by  defining  the  product  a3  equals  with  y 
which  will  be  the  third  variable  parameter. 

After  transforming  a  non-linear  characteristic 
equation  in  a  linear  characteristic  equation  as  it  was 
described  above  then  all  these  which  have  been  stated  for 
the  linear  case  can  be  applied  for  the  transformed 
characteristic  equation. 
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XI.    CONCLUSIONS  AND  RECOMMENDATIONS 


FOR  FURTHER  WORK 


This  thesis  has  developed  practical  methods  for  de- 
signing singular  compensators,  i.e.,  compensators  which 
produce  the  desired  singular  conditions  in  a  single  input- 
single  output  control  system. 

The  key  to  the  development  of  these  methods  was  the 
form  of  the  characteristic  equation  which  was  considered, 
i.  e: 

af1Cs)+ef2Cs)+f3(s)  =  0 

where:  f,  (s),  f2(s)  an^  ^Cs)  are  polynomials  of  s  and 
a, 6  are  real  variable  parameters. 

Necessary  and  sufficient  conditions  for  the  above 
characteristic  equation  to  be  singular  were  determined  and 
the  singular  line  theory  was  then  investigated  to  make 
possible  the-  development  of  the  required  standard  set  of 
rules  for  the  design  engineer. 

Application  of  singular  line  theory  in  control  system 
design  led  to  the  conclusion  that  it  can  be  used  in 
general  for  design  of  compensation  of  any  single  input 
single  output  linear  control  system  with  the  following 
basic  advantages. 
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Required  filters  for  any  desired  compensation 
structure  have  transfer  functions  which  are  known  as  a 
function  of  the  two  parameters  a   and  6.   The  parameter 
plane  diagram  required  for  analysis  and  synthesis  of  a 
design  problem  consists  of  straight  lines,  i.e.,  the 
singular  lines  which  can  be  drawn  by  hand,  without  needing 
a  computer.   Design  techniques  and  procedures  are  charac- 
terized by  simplicity,  flexibility  and  accuracy.   Finally, 
tolerances  on  the  parameter  a  and  3  values  can  be  pre- 
cisely established. 

When  considering  the  sensitivity  of  the  transient 

response  to  parameter  variations,  it  was  found  that  the 

control  system  is  less  sensitive  to  variations  in  a  and  3 

if  the  operating  point  on  the  singular  line  is  remote  from 

the  point  at  which  the  singular  line  is  tangent  to  the 

constant  t   and  a        curves, 
s      ns 

The  basic  design  tools  of  the  singular  line  theory 
design  method  are:   The  Root  locus  of  the  plant  and  the 
parameter  plane  diagram  of  the  compensated  singular 
system,  i.e.,  under  certain  conditions  or  in  other  words 
singular  conditions,  a  classical  design  method  [The  Root 
locus  method) ,  and  a  modern  design  method  [the  parameter 
plane  method)  can  work  together  in  such  a  way  that  each 
compensates  disadvantages  of  the  other. 
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For  example  stability  analysis  which  is  not  a  simple 
process  in  the  parameter  plane  is  performed  on  the  Root 
locus  diagram  in  the  s-plane  or  studying  the  effect  of 
variations  in  either  one  or  both  of  the  parameter  a  and  3 
values  on  the  transient  response  of  the  system,  which  is 
not  a  simple  process  on  the  Root  locus  diagram  is  per- 
formed in  the  parameter  plane,  etc. 

Conceptually  the  compatibility  of  the  Root  locus 
classical  design  method,  which  is  adequate  when  there  is 
a  single  variable  parameter  with  the  parameter  plane 
modern  design  method,  which  is  adequate  when  there  are 
two  variable  parameters,  when  singular  conditions  exist, 
is  based  on  the  fact  that  in  this  case  one  of  the  para- 
meters can  be  defined  through  the  equation  of  a  singular 
line,  as  a  function  of  the  other.   Thus  a  two  parameter 
problem  becomes  a  one  parameter  problem.   The  similarity 
also  of  the  Root  locus  of  the  plant  with  the  Root  locus 
of  the  compensated  singular  system,  independently  of  which 
singular  line  is  considered  generalized  the  validity  of 
the  above  concept  and  made  possible  the  development  of 
graphical  design  techniques  which  simplify  the  work  and 
markedly  increase  the  power  of  visualization  of  the  de- 
signer. 

It  was  also  realized  that  the  singular  line  theory 
could  be  applicable  in  self-adaptive  control  system,  since 
movement  of  the  operating  point  (a,  3)  along  a  dominant 
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line  segment  of  a  singular  line  or  within  a  dominant  root 
region  in  the  parameter  plane  implies  almost  fixed  system 
performance.   In  such  a  case  the  adaptive  controller  need 
only  to  drive  the  system  back  to  the  dominant  root  line 
segment  or  dominant  root  region  and  not  to  a  specific 
operating  point  in  the  parameter  plane. 

The  singular  line  theory  was  further  extended  in  the 
parameter  space  and  in  the  most  general  case  where  products 
of  the  variable  parameters  appear  in  the  characteristic 
equation. 

Further  study  is  needed  for  the  case  of  a  singular 
system  with  a  negative  error  constant.   This  study  can  be 
initialized  from  an  analysis  of  equation  (8-25) . 

Also  the  correlation,  if  any,  of  a  singular  feedback 
compensator  with  the  linear  state  variable  feedback  maybe 
interesting . 

Some  other  related  areas  for  future  work,  which  are 
recommended  are:  Applications  of  the  singular  line  theory 
in  self-adaptive  control  systems.   Singular  surface  theory 
and  control  systems,  singular  line  theory  and  multi-input, 
multi-output  linear  control  systems,  and  finally  any  of 
the  above  in  the  general  case  where  products  of  the  variable 
parameters  also  appear  in  the  characteristic  equation,  i.e., 
non-linear  cases. 
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APPENDIX  H 
SINGULAR  POINT  PROGRAM 


The  purpose  of  this  program  is  to  determine  the  roots  of  a 
characteristic  equation  which  are  singular  points  with  respect  to  the 
parameters  a   and  p. 

This  program  is  applicable  to  polynomials  whose  coefficients  are 
of  the  form  (b  a+cg+d),  where  a  and  e  are  variable  parameters  and  b,  c, 
and  d  are  real  constants. 

For  a  given  characteristic  equation  this  program  solves  the 
system  (2-4),  by  the  Cramer's  rule  method  for  a  desired  value  of  c, 
and  for  a  range  of  values  for  u>n,  which  has  been  determined  within 
the  program  and  can  be  changed  when  it  is  desired,  by  changing  the 
card  0035. 

The  four  columns  which  are  generated  at  the  print  output  are  the 
following:  The  first  is  a  list  of  the  values  for  wn,  which  were  used 
for  the  calculations  and  the  other  three  are  the  corresponding  nume- 
rators and  common  denominator  from  the  Cramer's  rule  solution  of  the 
system  (2-4). 

Any  value  of  un  in  the  first  column,  for  which  the  associated 
values  of  the  other  three  columns  are  zero's,  together  with  the  value 
of  t,   which  was  used  defines  a  singular  point  for  the  characteristic 
equation. 

The  input  data  for  the  program  are: 

a.  First  card: 

Value  of  ? and  N,  format  (F4.2,  12)  where  c  is  the  desired 
value  of  damping  ratio  and  N  is  the  order  of  the  characteristic  equation. 

b.  Second  card: 

Constant  coefficients  in  ascending  order  (8E  10.5) 

c.  Third  card: 

Alfa  coefficients  in  ascending  order  (8E  10.5) 

d.  Fourth  card: 

Beta  coefficients  in  ascending  order  (8E  10.5) 
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APPENDIX  I 
FIGURES" 


B %6d-Z >      G(s) 


FIGURE  (4-1) 
Block  diagram  of  a  unity  feedback  system 
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FIGURE  (4-2) 
First  compensation  structure 
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FIGURE  (4-3) 
Second  compensation  structure 
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FIGURE  (4-4) 
Third  compensation  structure 
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FIGURE  (4-5) 
Fourth  compensation  structure 


^-ir^ 


Cl 


(s) 


±o- 


H(s) 


FIGURE  (4-6) 
Fifth  compensation  structure 
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X-  SCALE -2.00E+0I     UNITS    INCH. 
Y-  SCALE"  6.00E*-0I     UNITS    INCH. 


FIGURE  (5-1) 
Example  ( 5-1 ): Parameter  plane  diagram, 
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X-SCALE  =  2 
Y-SCALE  =6 


Exampl e 


00  E+01  UNITS  INCH. 
00  E+01  UNITS  INCH. 
FIGURE  (5-2) 
(5-1 ): Parameter  plane  diagram  showing 
a  singular  line. 
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y- 


X-SCALE-5.00E-OI    UNITS    INCH. 
Y-SCALE'4.  00E*00  UNITS    INCH. 


FIGURE  (6-2) 
Example  ( 6-1 ): Parameter  plane  diagram 
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FIGURE  (6-3) 
Block  diagram  of  the  plant 


R  + 


Gc(s)=gFg(s) 


G(s)=KFg(s) 


H(s)=a 


FIGURE  (6-4) 
Block  diagram  of  the  compensated  singu- 
lar system  corresponding  to  the  "second 
singular  case" .scheme . 
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Singular  compensator  Gsr(s) 


1 


^<A 


rQ 


Gc(s)=Fg(s) 


H(s)=a 


G(s)=KFg(s) 


FIGURE  (6-5) 
Block  diagram  of  the  compensated  singular 
system  corresponding  to  the  modified  second 
singular  case, scheme. 


->       E 

GSC(S) 

G(s)=KFg(s) 
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f> 

FIGURE  (6-6) 
Simplified  block  diagram  of  the  modified 
second  singular  case, scheme. 
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Exampl e  (6-2) : 
sated  system. 


FIGURE  (6-7) 
Root  locus  of  the  uncompen 
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-003  -001 


X-SCALE-   I.OOE*02     UNITS    INCH. 
Y- SCALE*  I.OOE+-02     UNITS    INCH. 


FIGURE  (6-8) 
Example  (6-2):  Parameter  plane  di agram,K=20 
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FIGURE    (6-9) 
Example    (6-2):    Root    locus    of    the    compensated    singu- 
lar  system   corresponding    to   the   singular    root: 
s  =  -0.3884±j0.8899    or    (  Cs  =  0  .  40  ,%s  =  0  •  97096) 
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X-SCALE  »  I.OOE+02    UNITS     INCH. 
Y-SCALE  "  I.00E+02    UNITS     INCH. 


FIGURE  (6-10) 
Example  (6-2):  Parameter  plane  di agram  ,K=35 
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FIGURE    (5-1  1  ) 
Example      (6-3):    S-plane      diagram 
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FIGURE    (6-19) 
Exampl e(6-5  ) :    Root    locus    of   the    compensated    singular 
system   corresponding   to    the    singular    root    Us=0.30, 
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FIGURE  (6-24) 
Example  (6-5):  Root  locus  of  the  compensated  sin- 
gular system  corresponding  to  the  singular  root 
Us=0. 75,a>ns=0. 60667)  . 
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FIGURE  (6-27) 
Example  (6-6):  Block  diagram  of  the 
initially  compensated  plant. 
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FIGURE    (6-29) 
Example    (6-6):    Root    locus    of   the    uncompensated    system 
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FIGURE  (6-30) 
Example  (6-6):  Root  locus  of  the  plant  after 
it  was  compensated  by  an  initial  cascade  com- 
pensated with  transfer  f uncti on (s+2)/ (s+20 ) 
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FIGURE    (6-31  ) 
Example    (6-6):Root    locus    of   the    singular 
compensated    system    corresponding    to    the 
singular    root    ( 5    =0  .6,  u>      =1  .2326 )    or    (s: 
-0  .7395±j0  .9861  )s 
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X-SCALE=5.00E+01  UNITS  INCH. 
Y-SCALE=2.00E+02  UNITS  INCH. 

FIGURE  (6-49) 
Example  (6-6):  Parameter  plane  diagram  of  the 
compensated  singular  system. Si ngul ar  point 
(5S=0. 6, a)ns  =  l. 232586) 
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X-SCALE=4. OOE+02  UNITS  INCH. 
Y-SCALE=1.00E+03  UNITS  INCH. 

FIGURE  (6-50) 
Example  (6-6):  Parameter  plane  diagram  of  the 
compensated  singular  system. Si ngul ar  point 
Us  =  0.6,u>ns  =  1.2325) 
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FIGURE  (6-51  ) 
Example  (6-7):Root  locus  of  the  compensa 
ted  singular  system. 
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Exampl e (6-7) : Pa rameter  plane  diagram  and 
bounda  ry  lines. 
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FIGURE  (7-1) 
Example  (7-1):  Block  diagram 
initial  system  (Plant). 
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FIGURE  (7-2) 
Example  (7-2):  Block  diagram 
compensated  singular  system. 
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FIGURE  (8-1) 
First  modified  singular  compensation 
structure. 


R   i^   E 


<^ 


G(s)=KFg(s) 


FIGURE  (8-2) 
Second  modified  singular  compensation 
structure. 
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FIGURE  (8-3) 
Third  modified  singular  compensation 
structure . 
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FIGURE  (8-4) 
Fourth  modified  singular  compensation 
structure . 
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FIGURE  (8-5) 
Fifth  modified  singular  compensation 
structure . 


FIGURE  (8-6) 
Sixth  modified  singular  compensation 
structure . 
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